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Abstract 


Wave-theoretical  computer  codes  have  been  developed  to  model  pulse 
propagation  in  the  central  Arctic  Ocean.  Pulse  shapes  as  a  function  of  range 
and  depth  are  computed  from  the  Pulse  Fast  Field  Program  (PFFP)  and  the  pulse 
parabolic  equation  (PPE)  code.  Group-  and  phase-velocity  dispersion  and 
eigenfunctions  are  computed  from  the  PFFP  or  from  a  corresponding  normal-mode 
program.  Good  agreement  has  been  obtained  between  measured  and  computed  SOFAR 
signals.  The  effect  of  ice  roughness  on  Arctic  SOFAR  propagation  is  illus¬ 
trated  from  field  data  and  the  PFFP.  Hydroacoustic  signals  from  underwater 
explosions  that  have  propagated  over  the  Arctic  abyssal  plains  commonly 
display  marked  frequency  dispersion  in  pulses  that  are  bottom-interacting  and 
that  arrive  after  the  SOFAR  signal.  In  the  infrasonic  band  of  2  to  20  Hz,  the 
temporal  dispersion  for  each  pulse  that  has  interacted  with  the  flat  bottom  of 
the  plain  can  be  nearly  as  strong  as  that  observed  in  the  SOFAR  signal  for  the 
first  mode.  However,  the  bottom-interacting  pulses  correspond  to  a  coherent 
summation  of  many  higher-order  normal  modes  in  a  channel  bounded  above  by  the 
ocean  surface  and  below  by  the  zone  of  increasing  velocity  in  the  upper  400  m 
of  the  bottom  sediment.  Using  normal-mode  theory  and  the  Multiple  Scattering 
Pulse  Fast  Field  Program  (MSPFFP) ,  we  have  analyzed  the  dispersion  and  pulse 
shapes  and  have  derived  the  acoustic  properties  of  the  bottom  in  the  Pole, 
Barents,  and  Mendeleyev  Abyssal  Plains.  The  principal  properties  of  the 
bottom  controlling  the  propagation  are  compressional  velocity,  density,  and 
attentuation.  In  contrast,  the  ice  layer  has  a  negligible  effect  on  the 
dispersion  of  the  observed  waves.  The  effect  on  pulse  compression  of  this 
frequency  dispersion  of  the  bottom-interacting  signals  was  simulated 
numerically,  using  predistorted  waveforms  matched  to  the  dispersion  of  the 
SOFAR  channel  at  specified  ranges.  We  include  a  users  manual  (with  test 
computations)  for  our  computer  codes.  We  hope  that  this  manual  will  be  useful 
to  the  many  users  of  our  codes. 
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Introduction 

Our  research  in  Arctic  Marine  Acoustics  has  been  a  strong  interplay 
between  field  observations  and  computer  modeling  of  underwater  propagation. 
In  the  spring  of  1982  at  FRAM  IV/TRISTEN  we  completed  the  fourth  consecutive 
field  season  of  acoustic  measurements  in  the  Eurasian  Basin.  Project  TRISTEN 
for  methods  of  long-range  signaling  in  the  Arctic  channel  was  spawned  by  our 
research  in  long-range  propagation.  The  new  experiments  of  the  FRAM  series 
were  highly  successful,  and  they  provided  us  with  the  first  opportunity  to 
obtain  systematic  data  over  two  abyssal  plains  of  infrasonic  signals  (2  to 
20  Hz)  traveling  via  bottom-interacting  paths  to  long-ranges  in  the  Arctic 
channel.  The  frequency  dispersion  of  these  signals  was  of  practical  interest 
in  connection  with  TRISTEN,  but  our  main  objective  was  to  understand  the 
propagation  of  these  waves  in  the  Arctic  channel  and  to  derive  the  acoustic 
properties  of  the  upper  400  m  of  sediment  which  supports  the  propagation  of 
these  waves  to  long  ranges. 

The  new  experiments  also  provided  us  with  a  vastly  expanded  data  set  for 
a  systematic  investigation  of  environmental  factors  controlling  the  dispersion 
and  excitation  of  the  normal  modes  observed  in  Arctic  SOFAR  propagation.  We 
also  had,  for  the  first  time,  the  opportunity  to  investigate  long-range 
propagation  of  hydoracoustic  signals  from  earthquakes  at  the  Nansen  Ridge 
(Fig.  1). 

Simultaneously  with  the  field  measurements  and  data  analysis,  we  were 
improving  and  standardizing  our  full-wave  computer  models  for  transmission 
loss  and  pulse  propagation  to  make  them  as  automatic  and  fool-proof  as 
possible  so  that  we  could  conveniently  analyze  the  wave  propagation.  We  have 


completed  this  task  with  the  interim  standard  codes,  and  we  have  distributed 
the  most  recent  updates  of  the  combined  FFP  and  normal-mode  program  (and  the 
corresponding  pulse  code)  to  the  Navy  users  group. 

The  Arctic  SOFAR  signals  commonly  differ  markedly  in  charcter  from  those 
observed  in  the  deep  channel,  largely  because  of  the  predominance  of  very  low 
frequency  waves  in  the  Arctic  channel.  Waves  above  100  Hz  are  strongly 
scattered  by  the  rough  ice  boundaries,  and  thus  they  are  not  propagated  within 
the  waveguide  to  great  ranges.  Because  the  wavelengths  of  the  observed  waves 
are  commonly  comparable  with  the  dimensions  of  the  portion  of  the  duct 
controlling  the  propagation,  full  wave  theory  must  be  used  to  explain  the 
propagation,  particularly  for  waves  traveling  in  the  upper  thousand  meters  of 
the  water  column  and  for  those  waves  with  turning  points  in  the  upper  400  m  of 
bottom  sediment.  The  observations  of  marked  frequency  dispersion  of  explosion 
sounds  at  long  ranges  motivated  the  development  of  the  present  computer  codes. 

In  this  report  we  summarize  our  research  in  two  invited  papers  presented 
at  the  Acoustial  Society  of  America.  We  include  a  users  manual  (with  test 
computations)  for  our  computer  codes.  We  hope  that  this  manual  will  be  useful 
to  the  many  users  of  these  codes  within  the  Naval  community. 
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Abstract 


Wave-theoretical  computer  codes  have  been  developed  to  model  pulse 
propagation  in  the  central  Arctic  Ocean.  Pulse  shapes  as  a  function  of  range 
and  depth  are  computed  from  the  Pulse  Fast  Field  Program  (PFFP)  and  the  pulse 
parabolic  equation  (PPE)  code.  Group-  and  phase-velocity  dispersion  and 
eigenfunctions  are  computed  from  the  PFFP  or  from  a  corresponding  normal-mode 
program.  Good  agreement  has  been  obtained  between  measured  and  computed  SOFAR 
signals.  Computer  simulations  of  pulse  compression  show  in  detail  the 
possibility  of  achieving  significant  signal  gain  by  pulse  compression  in  the 
Arctic  channel.  This  is  accomplished  by  using  predistorted  waveforms  matched 
to  the  dispersion  of  the  channel  at  specified  ranges.  The  effect  of  ice 
roughness  on  Arctic  SOFAR  propagation  is  illustrated  from  field  data  and  the 
PFFP.  Recent  experiments  of  long-range  transmission  of  explosion  sound  over 
the  Pole  Abyssal  Plain  show  marked  frequency  dispersion  of  bottom-interacting 
signals  following  the  SOFAR  signal.  The  dispersion  of  each  pulse  interacting 
with  the  flat  ocean  bottom  of  the  Plain  at  a  depth  of  4200  m  can  be  nearly  as 
strong  as  that  observed  in  the  SOFAR  signal  for  the  first  mode,  but  the 
bottom- interacting  pulses  appear  to  correspond  to  very  high-order  normal 
modes.  Present  efforts  are  described  to  model  these  bottom-interacting 
pulses,  as  well  as  T-phase  propagation  from  earthquakes  at  the  Nansen  Ridge. 
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Introduction 


The  two  features  of  the  polar  environment  that  most  strongly  in¬ 
fluence  underwater  sound  are  the  permanent  ice  cover  and  the  velocity 
structure  in  the  water.  The  ice  modifies  propagation,  particularly  at 
high  frequencies,  by  scattering  waves  from  the  rough  ice  boundaries. 
Sound  velocity  is  increasing  with  depth  from  the  surface  to  the  bottom. 
In  sufficiently  deep  water,  sounds  are  transmitted  to  great  ranges  in 
this  sound  channel  by  upward  refraction  in  the  water  and  repeated  re¬ 
flection  from  the  ice  canopy.  The  surface  sound  channel  of  the  Arctic 
is  the  polar  extension  of  the  deep  sound  channel  or  SOFAR  channel  of 
the  nonpolar  oceans  (Ewing  and  Worzel,  1948) ,  but  the  Arctic  signals 
are  often  quite  different  in  character  from  those  observed  in  the  deep 
channel,  largely  because  of  the  predominance  of  low-frequency  waves 
(5  to  40  Hz)  in  the  Arctic  channel. 

The  experiments  of  SOFAR  propagation  in  the  central  Arctic  Ocean 
span  a  twenty-three  year  period.  This  work  was  initiated  by  NUSC  in 
1958-59.  Lamont  cooperated  in  the  field  with  NUSC  during  these  early 
experiments. 

These  early  experiments  revealed  the  unique  features  of  Arctic 
SOFAR  propagation.  Along  a  1000  km  path,  the  explosion  signals  were 
confined  to  a  frequency  band  from  about  7  to  35  Hz.  The  signals  showed 
marked  frequency  dispersion.  The  first  two  normal  modes  were  commonly 
observed.  These  waves  were  successfully  interpreted  by  normal-mode 
theory  for  the  Arctic  channel  (Kutschale,  1961). 

Frequency  dispersion  of  explosion  sounds  in  shallow  water  was  dis¬ 
covered  by  Ewing  (1936)  and  was  extensively  investigated  by  Worzel  and 
Ewing  (1948)  during  World  War  II.  These  observations  of  Worzel  and 
Ewing  in  shallow  water  were  explained  by  Pekeris  (1948)  who  had  devel¬ 
oped  the  normal-mode  theory  of  shallow-water  propagation.  There  is  a 
close  analogy  between  propagation  in  shallow  water  and  propagation  to 
long  ranges  in  deep  water  of  the  central  Arctic  Ocean.  In  the  central 
Arctic  Ocean,  the  dispersion  is  different  from  that  observed  in  shallow 
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water  because  the  Arctic  propagation  is  controlled  by  the  velocity 
structure  in  the  water/  but  the  methods  of  analyzing  the  Arctic  signals 
parallels  closely  the  methods  employed  by  these  early  workers. 


The  Waveguide 


The  central  Arctic  Ocean  is  covered  by  ice  the  year  round.  The 
pack  ice  averages  about  3  meters  thick.  The  ice  is  in  continual/  very 
slow  motion  under  the  influence  of  winds  and  currents.  The  pack  ice 
is  not  of  uniform  thickness,  but  it  is  broken  by  leads  of  open  water 
and  pressure  ridges.  The  irregularity  of  the  pack-ice  boundaries  is 
dependent  upon  location  and  time  of  year,  but  a  pattern  of  ice  rough¬ 
ness  useful  to  underwater  sound  has  been  found  for  broad  areas  of  the 
Arctic  Ocean  for  each  of  the  four  seasons  of  the  year  (Diachok,  1973) . 


This  seasonal  and  spatial  dependence  of  ice  roughness  is  of  great 
importance  for  long-range  sound  propagation  because  of  scattering  of 
waves  impinging  upon  the  ice.  We  illustrate  this  effect  in  Figures  1 
and  2.  In  Figure  1  we  show  a  computer  simulation  of  pulse  shapes  as  a 
function  of  range  with  no  ice  roughness  in  the  model.  In  Figure  2  we 
do  the  same  simulation  but  with  an  RMS  ice  roughness  of  20  m.  The 
amplitude  scale  is  the  same  as  in  the  previous  Figure  1,  but  we  note 
the  marked  decrease  in  the  amplitudes  of  the  waves,  particularly  the 
high-frequency  components  of  waves  corresponding  to  the  first  normal 
mode.  We  have  never  observed  such  rough  ice  in  the  Arctic,  but  this 
extreme  example  was  chosen  to  illustrate  a  dramatic  effect  on  the 
waves  at  moderately  short  ranges. 


In  contrast  to  the  variability  of  the  ice  cover,  the  sound-speed 
structure  with  depth  in  the  water  is  very  uniform  as  a  function  of 
location  and  season.  Sound  speed  is  an  increasing  function  of  depth 
from  the  surface  to  the  bottom.  Sounds  are  transmitted  to  long  ranges 
in  this  sound  channel  by  upward  refraction  in  the  water  and  repeated 
reflection  from  the  ice.  Figure  3  shows  a  typical  Arctic  profile.  The 
waveguide  consists  basically  of  three  parts:  the  upper  layer,  the 
lower  layer,  and  the  bottom  interface.  The  strong  gradient  to  a  depth 
of  300  m  is  controlled  by  the  rapid  increase  of  temperature  to  this 
depth,  while  the  deeper,  weaker  gradient  is  controlled  by  the  pressure 
effect.  The  profile  is  approximately  bilinear  with  the  sharp  decrease 
in  the  gradient  near  350  m  in  the  Amerasian  Basin  and  300  m  in  the 
Eurasian  Basin  (Crary,  1956) . 


Figure  4  shows  typical  ray  paths  computed  at  one-degree  intervals 
from  a  source  near  the  axis  of  the  Arctic  channel.  The  deep  penetrat¬ 
ing  RSR  waves  arrive  first,  and  the  waves  traveling  near  the  axis  last. 
From  this  diagram  we  see  that  the  onsets  of  the  signals  are  variable 
with  range  in  amplitude  and  mean  horizontal  velocity  in  response  to  the 
cycling  of  the  deep  penetrating  RSR  sounds,  but  that  waves  traveling 
near  the  surface  terminate  the  pulse  abruptly.  This  is  shown  in 
Figure  5  by  a  detailed  computer  simulation  of  pulse  shapes  as  a  func¬ 
tion  of  range.  Topographic  features  will  block  the  deep-penetrating 
RSR  waves.  The  highest  frequency  waves  correspond  to  the  deep-pene¬ 
trating  RSR  paths.  These  waves  make  fewer  contacts  with  the  ice  per 
unit  horizontal  distance  than  the  shallower  paths.  The  bottom  reflec- 


tions  commonly  correspond  only  to  the  lowest  frequency  waves  observed 
in  the  signal  because  of  attenuation  of  waves  at  both  boundaries  of 
the  waveguide.  The  principal  complications  for  making  detail  compari¬ 
sons  between  experiment  and  theory  are  introduced  by  wave  scattering 
from  the  ice  and  the  bottom. 

If  the  bottom  is  irregular  along  the  propagation  path,  we  generally 
neglect  the  bottom  interactions  in  the  computer  simulations.  This  as¬ 
sumption  is  justified  from  field  data  such  as  shown  in  Figure  6.  For 
this  series  of  experiments,  shots  of  identical  size  and  depth  were 
launched  from  T-3  and  recorded  aboard  ARLIS  II  employing  a  hydrophone  at 
a  fixed  depth.  Note  the  weakening  of  the  deep-penetrating  waves  as  the 
amount  of  shallow  water  along  the  path  increases.  This  experiment  shows 
qualitatively  that  the  onset  of  a  signal  is  largely  determined  by  waves 
which  pass  over  all  bottom  topography  without  suffering  bottom  reflec¬ 
tions.  Following  the  abrupt  termination  of  the  SOFAR  signal,  we  see 
weak  waves  scattered  from  the  bottom.  Propagation  of  these  incoherent 
scattered  waves  is  beyond  the  capabilities  of  our  present  computer  codes 
which  model  only  the  coherent  portion  of  the  pulses. 

At  low  frequencies  theories  for  the  coherent  component  of  the  scat¬ 
tered  waves  from  rough  surfaces  account  adequately  for  attenuation  of 
waves  reflected  from  the  ice.  Both  the  formulation  by  Kuperman  (1975) 
and  the  one  by  Mellen  and  Marsh  (1965)  are  used  in  the  computer  codes. 


Data  and  Some  Comparisons  with  Theory 

Figures  7  and  8  show  typical  Arctic  SOFAR  signals  at  a  range  of 
1300  km  and  720  km  from  the  shot.  Three  wave  groups  are  shown  clearly 
by  Figure  8.  The  first  group  corresponds  to  the  deep-penetrating  RSR 
waves.  This  group  lasts  about  four  seconds  on  this  oscillogram.  The 
first  two  pulses  of  this  group  are  well  separated  in  time,  but  toward 
the  end  of  this  group,  at  four  seconds  on  the  oscillogram,  it  is  diffi¬ 
cult  to  resolve  the  individual  arrivals.  Following  this  group  of  waves 
we  observe  a  strong  train  of  nearly  sinusoidal  waves  in  Figures  7  and  8 
showing  clear  frequency  dispersion,  with  low-frequency  waves  traveling 
faster  than  high-frequency  waves  in  this  wavetrain.  These  waves  corre¬ 
spond  to  the  first  normal  mode  propagated  near  the  axis  of  the  channel. 
Following  the  abrupt  termination  of  the  SOFAR  signal,  we  observe  weak 
bottom  reflections  rising  above  an  oscilliatory  wavetrain  corresponding 
to  waves  scattered  from  an  irregular  bottom. 

If  the  bottom  is  flat  over  the  whole  propagation  path,  we  observe 
many  strong  coherent  bottom-interacting  pulses  at  long  ranges.  These 
pulses  following  the  SOFAR  signal  can  show  masked  frequency  dispersion. 
This  is  shown  in  Figure  9  from  a  typical  signal  propagated  along  a 
316  km  path  over  the  Pole  Abyssal  Plain  which  is  at  a  depth  of  4200  m 
below  the  surface.  The  oscillogram  and  the  spectrogram  are  not  quite 
on  the  same  time  scale,  and  thus  we  have  indicated  corresponding  ar¬ 
rivals  on  the  two  displays  of  the  signal.  The  upper  trace  of  the  oscil 
logram  corresponds  to  a  slightly  lower  amplifier  gain  than  the  lower 
trace.  In  this  Figure  9  we  clearly  see  the  frequency  dispersion  of  the 
first  two  normal  modes  of  the  SOFAR  signal.  Following  the  abrupt  termi¬ 
nation  of  the  SOFAR  signal,  we  observe  sixteen  strong  bottom- interacting 
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pulses.  The  total  duration  of  the  SOFAR  signal  and  these  bottom¬ 
interacting  pulses  is  nearly  fifty  seconds.  The  bottom-interacting 
pulses  near  the  end  of  the  signal  are  nearly  sinusoidal. 
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In  Figure  10  we  show  in  detail  the  waveforms  and  spectra  of  two  of 
the  bottom  interacting  pulses  of  the  preceding  Figure  9.  The  signals 
are  normally  dispersed,  with  low-frequency  waves  traveling  faster  than 
high-frequency  waves.  We  are  presently  trying  to  explain  the  dispersion 
and  pulse  shapes  of  these  waves  from  a  model  which  extends  the  vertical 
velocity  profile  of  the  ocean  several  hundred  meters  into  the  bottom  of 
the  Plain.  We  believe  that  each  bottom-interacting  pulse  corresponds  to 
portions  of  dispersion  curves  for  very  high-order  normal  modes  derived 
from  such  a  model.  This  conclusion  is  based  upon  our  many  years  of  ex¬ 
perience  in  analyzing  dispersion  of  the  normal-modes  of  the  SOFAR  sig¬ 
nals. 


The  sound  spectrograms  of  Figures  11  and  12  show  the  frequency  dis¬ 
persion  of  two  typical  SOFAR  signals.  We  see  the  dispersion  clearly  of 
waves  corresponding  to  the  first  and  second  normal  modes.  Waves  of  the 
second  normal  mode  were  difficult  to  see  on  the  preceding  oscillograms 
because  the  amplitudes  are  considerably  weaker  than  waves  of  the  first 
normal  mode.  The  sharp  pulses  on  Figure  8  corresponding  to  the  deep- 
penetrating  RSR  waves  show  horizontal  bands  across  the  spectrogram  of 
Figure  12,  but  a  close  examination  of  many  spectrograms  shows  that  these 
pulses  correspond  to  a  superposition  of  many  higher-order  normal  modes 
with  group  velocities  equal  to  the  mean  horizontal  velocities  of  each  of 
the  pulses.  We  should  point  out  that  these  spectrograms  give  only  a 
qualitative  measure  of  the  relative  amplitudes  of  the  waves.  Our  pri- 
many  purpose  for  this  type  of  display  is  to  resolve  the  normal  modes  of 
the  signal. 

In  Figure  13  observed  dispersion  is  compared  with  computed  disper¬ 
sion  lor  the  first  two  modes.  The  data  were  taken  from  Figure  12.  We 
have  converted  the  group- velocity  scale  of  the  computations  to  a  time 
scale  to  make  a  direct  comparison  with  the  observed  spectrogram  possible. 
The  measured  and  computed  dispersion  curves  virtually  coincide  for  both 
modes.  For  these  computations  we  have  used  the  sound-speed  profile  of 
Figure  3  which  was  measured  during  the  transmission  experiment  in  the 
Eurasian  Basin. 

Figure  14  shows  oscillograms  from  three  shots  which  illustrate  the 
effect  of  source  depth  on  waves  of  the  first  mode.  We  note  that  waves 
corresponding  to  the  first  normal  mode  are  weak  for  the  366-meter  shot. 
This  is  due  to  the  pressure  distribution  with  depth  of  the  first  normal 
mode.  Each  oscillogram  corresponds  to  one  output  from  a  hydrophone  and 
one  from  a  vertical-component  geophone  on  the  ice  surface. 

Figure  15  shows  reasonably  good  agreement  between  theory  and  exper¬ 
iment  of  the  depth  variation  of  pressure  for  the  first  normal  mode  at 
the  two  frequencies  of  20  and  40  Hz. 


In  Figure  16  we  show  the  effect  of  ice  roughness  on  the  peak  ampli¬ 
tudes  of  the  waves  of  the  first  normal  mode  as  a  function  of  range. 

This  experiment  was  conducted  in  the  Amerasian  Basin  over  the  Canada 
Abyssal  Plain.  We  analyzed  the  waveforms  in  the  band  from  25  to  70  Hz 
because  the  effect  of  ice  roughness  on  propagation  is  far  stronger  in 
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this  band  than  at  lower  frequencies.  The  shots  were  dropped  from  a 
Navy  aircraft  on  a  straight  path  to  the  listening  station.  The  sound 
sources  were  1.8  lb  SUS  charges  fired  at  800  feet.  When  we  compare 
the  peak  intensity  as  a  function  of  range  for  waves  of  the  first  normal 
mode,  we  get  good  agreement  between  the  measurements  and  the  computer 
simulations  of  the  pulses  if  we  use  an  RMS  underice  roughness  of  3  m. 

It  was  gratifying  to  learn  subsequent  to  these  experiments  that  about 
the  same  value  of  RMS  roughness  had  been  derived  in  this  area  from  di¬ 
rect  measurements  of  ice  roughness. 


Pulse  Compression 


The  uniformity  of  dispersion  data  of  waves  corresponding  to  the 
first  normal  mode  in  the  frequency  band  8  to  30  Hz  suggests  the  possi¬ 
bility  of  achieving  significant  signal  gain  by  pulse  compression  in 
the  Arctic  channel  (DiNapoli,  Viccione,  and  Kutschale,  1978) .  This  is 
shown  from  computer  simulations  by  using  predistorted  waveforms 
matched  to  the  dispersion  of  the  channel  at  specified  ranges. 


The  pulse  compression  technique  was  originally  developed  for  radar 
systems  to  obtain  high  range  resolution  with  long  pulses.  In  that  ap¬ 
plication,  the  frequency  of  the  transmitted  pulse  is  varied  with  time 
to  produce  a  wideband  waveform.  A  matched  filter  is  then  implemented 
at  the  receiver  by  passing  the  waveform  through  a  filter  with  the  in¬ 
verse  frequency  versus  time  relationship.  All  frequency  components  of 
the  pulse  are  thus  properly  delayed  and  appear  simultaneously  at  the 
output  to  produce  a  pulse  of  shorter  duration  and  high  peak  amplitude 
(Walter,  1961). 


The  technique  was  first  used  in  underwater  acoustics  by  Parvalescu 
and  Clay  (1965)  to  examine  signal  transmission  stability  in  Tongue  of 
the  Ocean.  Their  experiments  were  performed  at  400  Hz  in  1.8  km  of 
water  at  a  range  of  36  km.  They  discovered  that  the  signal  enhancement 
was  realizable  but  sensitive  to  geometry  (source  and  receiver  depth 
variations).  The  Arctic  offers  a  far  better  channel  (one  strong  mode 
at  15  Hz  versus  over  a  hundred  modes  at  400  Hz)  for  the  application  of 
this  scheme. 


Figure  17  shows  computed  SOFAR  signals  in  the  range  interval  from 
300  to  500  km.  The  source  and  hydrophone  are  at  a  depth  of  100  m. 

Near  the  beginning  of  the  signals  the  second  mode  is  strongest,  followed 
by  the  strong  waves  of  the  first  normal  mode.  Note  the  clear  frequency 
dispersion  in  the  signals  which  was  shown  in  Figure  13  to  be  in  close 
agreement  with  field  data.  The  peak  amplitude  of  the  waves  decreases  by 
spherical  spreading  since  the  effect  of  ice  roughness  is  small  at  these 
frequencies  to  ranges  of  500  km.  We  have  assumed  in  these  simulations 
that  the  propagation  path  was  over  an  irregular  bottom  so  that  the 
bottom-interacting  signals  shown  previously  from  the  Pole  Abyssal  Plain 
could  be  neglected. 

The  computer  codes  are  highly  automated,  and  we  can  thus  investi¬ 
gate  pulse  compression  as  a  function  of  range  by  reversing  in  time 
pulses  such  as  those  shown  in  Figure  17  and  propagating  them  back  through 
the  channel.  Figure  18  shows  tne  evolution  of  pulse  compression  as  a 
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function  of  range.  The  input  signal  at  the  source  is  the  time  reverse 
of  the  400  km  signal  of  Figure  17.  The  ocean  filter  is  matched  at 
400  km.  This  is  shown  by  the  sharp  pulse  at  this  range.  We  note  that 
the  peak  amplitude  of  the  400  km  signal  is  greater  than  the  peak  am¬ 
plitude  of  the  300  km  signal.  Thus  the  desired  effect  has  been  achieved. 
At  500  km  the  signal  rapidly  decreases  in  amplitude  as  the  pulse  spreads 
apart.  At  300  km  the  dispersion  in  the  main  pulse  is  opposite  to  that 
characteristic  of  the  channel,  but  at  500  km  the  pulse  has  started  to 
spread  with  the  normal  dispersion  of  the  channel.  The  sidelobes  of  the 
compressed  pulse  correspond  to  crosscorrelations  between  the  modes. 

Note  the  frequency  dispersion  in  the  sidelobes.  This  dispersion  re¬ 
verses  on  either  side  of  the  central  pulse. 

If  we  had  included  all  the  bottom  interactions  in  the  simulations, 
we  would  expect  to  observe  many  more  sidelobes  corresponding  to  cross¬ 
correlations  between  the  normal  modes  of  the  water  column  and  those  af¬ 
fected  by  the  bottom.  The  effect  of  bottom- interacting  signals  on  pulse 
compression  is  currently  under  investigation. 

The  degradation  with  range  of  the  time  reversed  waveform  is  graph¬ 
ically  depicted  in  Figure  19  by  plotting  the  peak  amplitude  of  the  re¬ 
ceived  waveform  versus  range.  With  no  special  signal  shaping,  the  peak 
amplitude  would  decay  with  range  according  to  20  log  R.  This  is  due  to 
the  combined  effects  of  cyclindrical  spreading  and  also  a  stretching  of 
the  waveform  in  time.  The  initial  decay  of  the  peak  amplitude  with 
range  for  the  time-reversed  signal  is  seen  to  obey  a  10  log  R  dependence. 
This  occurs  because  as  the  range  increases  the  pulse  is  trying  to  com¬ 
press  itself  and  one  is  left  with  only  cyclindrical  spreading.  Beyond 
the  matched  range,  the  pulse  again  begins  to  spread  in  time  and  the  peak 
amplitude  obeys  the  20  log  R  decay  law. 
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Appendix  A 


Computer  Codes 

Pulse  shapes  as  a  function  of  range  and  depth  are  computed  from 
the  Pulse  Normal  Mode  Program  (PNMP)  (Kutschale,  1969  and  1970) ,  the 
Pulse  Fast  Field  Program  (PFFP)  (Kutschale  and  DiNapoli,  1977) ,  and 
the  Pulse  Parabolic  Equation  (PPE)  code  (Kutschale  and  Tappert,  1977) . 
Group  -  and  phase  -  velocity  dispersion  and  eigenfunctions  are  com¬ 
puted  from  the  PFFP  or  from  the  corresponding  normal  -  mode  program. 

Each  of  these  pulse  codes  uses  Fourier  synthesis  via  the  fast  Fourier 
transform  to  coherently  add  up  the  frequency  response  at  each  frequency 
point  of  the  broad-band  input  signal  to  compute  the  temporal  waveforms 
as  a  function  of  range. 

I  will  not  go  into  the  mathematical  details  of  each  of  the  models. 

I  have  given  a  list  of  references  which  provide  some  of  the  details. 

One  might  ask  why  have  we  developed  three  computer  codes.  The  answer 
is  that  each  code  enjoys  certain  advantages  over  the  others  depending 
upon  the  application. 

For  example,  the  PPE  will  accommodate  a  range  dependent  environ¬ 
ment,  and  we  can  investigate  the  effects  on  the  waveforms  of  small 
variations  in  the  observed  sound-speed  profiles  as  well  as  model  a 
range-dependent  ice  roughness.  However,  for  most  applications  we  can 
neglect  lateral  variations  in  the  sound-speed  profiles  and  the  ice 
roughness. 

The  PFFP  and  the  PNMP  are  each  based  upon  the  same  integral  repre¬ 
sentation  of  the  sound  field  in  a  horizontally  stratified  ocean.  Either 
of  these  codes  can  accommodate  a  stratified  medium  of  a  mixture  of  solid 
layers  and  liquid  layers.  Thus  we  can  investigate  the  effects  of  the 
ice  sheet  on  the  propagation,  or  we  can  calculate  waveforms  measured  by 
geophones  on  or  in  the  ice.  Furthermore,  these  codes  are  appropriate 
for  investigating  the  bottom- interacting  signals  since  no  paraxial  ap¬ 
proximation  is  made  to  the  full  wave  equation.  In  the  PFFP  the  integral 
representation  for  each  frequency  component  is  integrated  directly  over 
wavenumber,  while  in  the  PNMP  the  integral  representation  for  each  fre¬ 
quency  component  of  the  pulse  is  expanded  in  a  normal-mode  series  plus  a 
branch- line  integral.  The  contribution  of  the  branch- line  integral  is 
negligible  at  long  ranges.  If  we  obsdrve  only  a  few  modes  and  if  we 
want  waveforms  at  only  a  few  range  points,  the  PNMP  is  very  fast  com¬ 
pared  to  the  PPE  or  PFFP,  which  effectively  computes  all  modes  over  the 
specified  range  of  phase  velocities  and  frequencies. 
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In  the  PNMP  or  PPE  the  scattering  effect  of  the  ice  on  the  ampli- 
tudes  of  the  waves  in  the  pulses  is  modeled  by  using  a  modified  formula  -j 

of  Mellen  and  Marsh  (1965).  In  the  PFFP  we  use  either  this  approach  or,  ^ 

since  damping  in  each  layer  is  a  natural  feature  of  this  code,  we  use  a  £ 

very  convenient  formulation  of  scattering  by  Kuperman  (1975) . 


In  Figures  1A  and  2A  we  show  by  comparison  of  waveforms  the  effects 
on  the  signals  of  neglecting  some  of  the  modes.  For  the  PFFP  simulation 
we  have  included  the  full  contribution  to  the  signal  of  waves  with  phase 
velocities  between  the  surface  sound  velocity  and  1490  m/sec  in  the 


sound-speed  profile  of  Figure  3.  In  the  PNMP  simulation  we  have  com¬ 
puted  the  signal  corresponding  to  only  the  first  four  modes.  We  see 
from  this  comparison  that  the  ends  of  the  wavetrains  are  identical  but 
that  the  beginnings  differ  in  detail  because  in  order  to  model  a  sharp 
pulse,  corresponding  to  the  deep-penetrating  RSR  waves  shown  at  a 
range  of  50  km  in  Figure  1A  and  on  the  measured  signal  of  Figures  8 
and  12,  we  must  include  all  the  higher  modes.  When  this  is  done,  both 
the  PFFP  and  PNMP  give  identical  results. 
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Figure  1. 
Figure  2. 
Figure  3. 


Captions  of  Figures 


Computed  Arctic  SOFAR  signals.  No  ice  roughness. 

Computed  Arctic  SOFAR  signals.  RMS  ice  roughness  20  m. 

Sound  speed  as  a  function  of  depth  in  Eurasian  Basin  of  the 
Bastern  Arctic  Ocean. 

Figure  4.  Ray  paths  computed  at  a  one-degree  interval. 

Figure  5.  Computed  Arctic  SOFAR  signals  as  a  function  of  range. 

Figure  6.  Signals  showing  the  effect  of  bathymetry. 

Figure  7.  Waveform  of  typical  Arctic  SOFAR  signal  measured  in  Amer- 

asian  Basin.  Range  1300  km. 

Figure  8.  Waveform  of  typical  signal  measured  in  Eurasian  Basin.  Range 
720  km. 

Figure  9.  Waveform  and  corresponding  spectrogram  of  typical  signal  mea¬ 
sured  over  flat  bottom  of  the  Pole  Abyssal  Plain.  Water 
depth  4100  m.  Range  316  km. 

Figure  10.  Waveforms  and  corresponding  spectrogram  of  two  bottom  inter¬ 
acting  signals  of  Figure  9. 

Figure  11.  Sound  spectrogram  of  typical  SOFAR  signal  measured  in  Amer- 
asian  Basin.  Range  870  km. 

Figure  12.  Sound  spectrogram  of  typical  signal  measured  in  Eurasian 
Basin.  Range  720  km. 

Figure  13.  Comparison  of  measured  and  computed  group-velocity  dispersion 
for  first  two  normal  modes. 

Figure  14.  Waveforms  as  a  function  of  depth.  Range  about  600  km. 

Figure  15.  Comparison  of  measured  and  computed  vertical  pressure  disper¬ 
sion  for  waves  of  the  first  normal  mode. 

Figure  16.  Comparison  of  experiment  and  theory  of  peak  intensity  as  a 

function  of  range  of  waves  of  the  first  normal  mode.  RMS  ice 
roughness  3  m. 

Figure  17.  Computed  SOFAR  signals  as  a  function  of  range. 

Figure  18.  Waveforms  of  pulse  compression  as  a  function  of  range. 

Figure  19.  Peak  amplitude  of  compressed  pulse  as  a  function  of  range. 
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Figure  1A.  Waveforms  computed  by  PFFP  corresponding  to  inclusion  of 
all  normal  modes. 


Figure  2A.  Waveforms  computed  by  PNMP.  Only  first  four  normal  modes 
computed. 
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Abstract 


Hydroacoustic  signals  from  underwater  explosions  that  have  propagated 
over  the  Arctic  abyssal  plains  commonly  display  marked  frequency  dispersion  in 
pulses  that  are  bottom-interacting  and  that  arrive  after  the  SOFAR  signal.  In 
the  infrasonic  band  of  2  to  20  Hz,  the  temporal  dispersion  for  each  pulse  that 
has  interacted  with  the  flat  bottom  of  the  plain  can  be  nearly  as  strong  as 
that  observed  in  the  SOFAR  signal  for  the  first  mode.  However,  the  bottom¬ 
interacting  pulses  correspond  to  a  coherent  summation  of  many  higher-order 
normal  modes  in  a  channel  bounded  above  by  the  ocean  surface  and  below  by  the 
upper  400  m  of  the  bottom  sediments,  where  the  velocity  increases  with  depth. 
Using  normal-mode  theory  and  the  Multiple  Scattering  Pulse  Fast  Field  Program 
(MSPFFP) ,  we  have  analyzed  the  dispersion  and  pulse  shapes  and  have  derived 
the  acoustic  properties  of  the  bottom  in  the  Pole,  Barents,  and  Mendeleyev 
Abyssal  Plains.  The  principal  properties  of  the  bottom  controlling  the  propa¬ 
gation  are  compressional  velocity,  density,  and  attenuation.  In  contrast,  the 
ice  layer  has  a  negligible  effect  on  the  dispersion  of  the  observed  waves. 
The  effect  on  pulse  compression  of  this  frequency  dispersion  of  the  bottom¬ 
interacting  signals  was  simulated  numerically,  using  predistorted  waveforms 
matched  to  the  dispersion  of  the  SOFAR  channel  at  specified  ranges. 
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Long-range  transmissions  of  explosion  sounds  over  the  Arctic  abyssal 
plains  commonly  show  marked  frequency  dispersion  in  bottom-interacting  pulses 
that  arrive  after  the  SOFAR  signal1.  Two  types  of  frequency  dispersion  are 
observed.  First,  the  frequency  content  of  each  successive  bottom  bounce  is 
observed  to  decrease  in  the  arrival  sequence.  The  central  frequency  of  each 
bottom-interacting  pulse  is  dispersed  from  about  25  Hz  at  low-order  bottom 
bounces  to  about  8  Hz  at  high-order  bottom  bounces.  Second,  we  observe  the 
dispersion  in  frequency  of  individual  bottom-refracted  pulses.  These  bottom¬ 
interacting  pulses,  which  will  be  shown  to  have  part  of  their  travel  paths 
through  the  bottom  sediment,  are  markedly  distorted  in  the  frequency  band  from 
5-20  Hz. 

Previous  investigations  have  revealed  the  unique  characteristics  of  the 
Arctic  SOFAR  signal2.  Fig.  1  shows  a  typical  sound  velocity  profile  in  the 
Arctic  channel.  Fig.  2  shows  that  sound  waves  are  transmitted  long  distances 
in  this  channel  following  successive  RSR  paths. 

Shot  signals  in  the  Arctic  channel  clearly  show  low-order  normal  modes. 
The  Arctic  SOFAR  signals  are  commonly  nearly  sinusoidal.  This  effect  is 
caused  by  the  constructive  interference  of  RSR  rays  that  have  traveled  in  the 
upper  few  hundred  meters  of  water,  where  the  sound  velocity  increases  rapidly 
with  depth.  Usually,  at  long-ranges,  it  is  the  first  two  normal  modes  that 
are  observed  in  the  frequency  band  below  40  Hz.  However,  Fig.  3  shows  that  if 
the  channel  is  deep  enough  along  the  entire  propagation  path,  such  as  over  an 
abyssal  plain,  the  early  part  of  the  arriving  SOFAR  wave  train  displays  a 
discrete  and  impulsive  character  corresponding  to  the  arrivals  of  deep- 
penetrating  RSR  rays,  while  the  latter  part  of  the  wave  train  maintains  its 
nearly  sinusoidal  character. 
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The  dispersive  character  of  the  Arctic  SOFAR  signal  is  shown  in  Fig.  4, 
which  is  a  sound  spectrograa  of  the  signal  shown  in  Fig.  3.  The  dispersion  of 
waves  corresponding  to  the  first  and  second  normal  modes  is  clearly  visible  in 
the  frequency  band  below  40  Hz.  The  impulsive  character  of  the  deep* 
penetrating  RSR  waves  is  displayed  by  horizontal  bands  across  the  spectrogram, 
corresponding  to  a  superposition  of  many  high-order  modes  with  group  velo¬ 
cities  equal  to  the  mean  horizontal  velocities  of  each  of  the  RSR  pulses. 
Fig.  5  shows  a  comparison  between  the  observed  and  computed  dispersion  of  the 
first  two  normal  modes.  Agreement  between  observation  and  theory  is 
excellent.  The  group-velocity  dispersion  of  low-order  modes  in  the  Arctic 
channel  is  very  uniform  over  a  long  propagation  path,  but  the  arrivals  of  the 
deep-penetrating  RSR  waves  vary  with  range  in  both  amplitude  and  mean  hori¬ 
zontal  velocity  according  to  their  cycle  ranges.  Because  of  scattering  at  the 
rough  boundaries  of  the  ice,  only  low  frequencies  can  propagate  to  long 
distances  in  the  Arctic  channel. 

An  interesting  observation  in  experiments  of  long-range  acoustic  trans¬ 
mission  over  the  abyssal  plains  is  the  marked  frequency  dispersion  of  the 
bottom-interacting  pulses  that  arrived  after  the  SOFAR  signal.  This  is  shown 
in  Fig.  6  for  a  typical  signal  propagating  along  a  316  km  path  over  the  Pole 
Abyssal  Plain  in  a  water  depth  of  4200  m.  In  Fig.  6  the  temporal  dispersion 
of  the  SOFAR  signal  is  shown  clearly.  There  are  seventeen  bottom-interacting 
pulses  corresponding  to  the  bottom  bounces  of  order  11  to  order  27  observed  in 
Fig.  6  that  arrived  after  the  SOFAR  signal.  The  total  duration  of  the 
SOFAR  signal  and  these  bottom-interacting  pulses  is  nearly  fifty  seconds. 
The  bottom-interacting  pulses  arriving  just  after  the  SOFAR  signal  are  very 
impulsive  in  character,  but  the  bottom-interacting  signals  near  the  end  of  the 
wave  train  are  nearly  sinusoidal.  Fig.  7  shows  in  detail  the  waveforms  of  two 


of  the  bottom-interacting  pulses  from  Fig.  6.  The  pulses  are  normally 
dispersed  with  low-frequency  waves  traveling  faster  than  high-frequency  waves 
and  the  temporal  dispersion  of  these  bottom-interacting  pulses  is  as  strong  as 
that  of  the  first  mode  of  the  SOFAR  signal. 

The  observed  long-range  bottom-interacting  pulses  may  be  resolved  by  ray 
acoustics  into  the  multipath  arrivals  scattered  from  the  top  and  the  bottom 
boundaries  of  the  Arctic  channel.  Thus  they  correspond  to  the  summation  of 
many  high-order  normal  modes  as  do  the  deep-penetrating  RSR  waves.  However, 
the  frequency  dispersion  of  individual  bottom-interacting  pulses  is  nearly  as 
strong  as  that  obseved  in  the  SOFAR  signal  for  the  first  mode.  These  dis¬ 
persed  bottom-interacting  pulses  must  have  traveled  through  and  been  distorted 
in  the  bottom  sediment.  This  portion  of  the  waveguide  must  be  very  uniform 
and  have  a  dimension  comparable  to  that  of  the  channel  in  the  upper  few 
hundred  meters  of  water,  such  that  it  can  support  the  propagation  of  sound 
waves  to  long  distances  with  frequencies  comparable  to  that  of  the  SOFAR 
signal. 

The  frequency  dispersion  of  low-order  normal  modes  has  often  been  used  to 
derive  the  acoustic  properties  of  sediment  beneath  a  shallow  ocean  waveguide 
since  normal  mode  theory  was  developed  by  Pekeris3'4.  In  the  present  work, 
efforts  have  been  made  to  explain  the  propagation  in  the  Arctic  channel  of 
dispersed  bottom- interacting  pulses  and  to  derive  the  acoustic  properties  of 
the  bottom  sediment  from  the  dispersion  of  high-order  normal  modes.  The 
techniques  used  in  this  study  include  normal-mode  theory  and  the  Multiple 
Scattering  Pulse  Fast  Field  Program  (MSPFFP) .  Normal-mode  theory  is  used  to 
explain  the  frequency  dispersion  of  the  bottom-interacting  wave  trains. 
Fig.  8  shows  that,  although  the  dispersion  of  individual  high-order  modes 
cannot  be  resolved  in  the  bottom- interacting  wave  train,  the  trend  of  the 


frequency  reduction  for  bands  present  in  the  bottom-interacting  sequence  is 
found  to  correlate  well  with  dispersion  profiles  of  very  high-order  normal 
modes.  Mode  numbers  as  high  as  60  may  be  observed  in  bottom-interacting 
signals  with  frequencies  lower  than  25  Hz.  The  MSPFFP  is  used  to  derive  a 
geoacoustic  model  that  predicts  synthetic  waveforms  that  match  well  with  the 
dispersion  profiles  of  the  bottom  interacting  signals.  The  effect  of  the  ice 
cover  and  shear  waves  in  the  sediment  was  investigated,  and  it  was  found  that 
these  two  factors  can  be  neglected  in  the  model. 

In  general,  the  bottom-scattered  wave  trains  show  two  types  of  pulses: 
high-frequency  discrete  arrivals  and  low-frequency,  strongly  distorted 
arrivals.  The  high-frequency  discrete  arrivals  correspond  to  the  post  critical 
bottom-reflected  pulses,  while  the  low-frequency,  strongly  distorted  arrivals 
correspond  to  the  bottom-refracted  pulses.  Because  of  attenuation  in  the 
bottom-sediment,  the  bottom-refracted  pulses  contain  lower  frequencies  than 
the  bottom-reflected  pulses. 

The  high-frequency  very  impulsive  arrivals  after  the  SOFAR  signal  are  the 
bottom-reflected  signals.  Beyond  a  certain  number  of  bottom  bounces,  the 
arrival  of  the  low-frequency,  distorted  signals  verifies  the  arrival  of  the 
bottom-refracted  signals.  In  the  bottom-refracted  signals,  two  types  of  pulse 
distortion  can  be  found  on  the  spectrogram  of  Fig.  8.  First,  the  pulses  are 
normally  dispersed  corresponding  to  multipath  propgation  of  pulses  through 
sediment  with  a  positive  velocity  gradient.  Second,  after  the  arrival  of  the 
normally-dispersed  waves,  pulses  arrive  which  show  a  branch  of  inverse 
frequency  dispersion.  These  bottom-scattered  waves  are  still  bottom-refracted 
waves,  but  they  have  penetrated  into  a  deeper  sediment  layer  with  a  different 
velocity  gradient.  Therefore,  the  prototype  model  which  supports  the  propa¬ 
gation  of  the  observed  bottom-interacting  signals  consists  of  three  parts: 


the  ocean,  the  shallow  sediment  with  a  strong  velocity  gradient  and  the  deep 
sediment  with  a  weaker  velocity  gradient. 

The  Multiple  Scattering  Pulse  Fast  Field  Program  decomposes  the  Pulse 
Fast  Field  Program5'*  into  ray  path  type  contributions7.  Each  decomposed  term 
can  be  interpreted  as  the  desired  path  contribution  for  a  corresponding 
bottom-interacting  pulse.  The  Fast  Field  Program  (FFP)  algorithm  integrates 
directly  the  full  wave  solution.  Thus,  the  MSPFFP  technique  is  a  natural 
scheme  with  which  to  model  the  bottom-interacting  pulses,  which  correspond  to 
the  coherent  summation  of  many  modes  over  a  limited  time  interval.  Temporal 
waveforms  are  computed  by  Fourier  synthesis. 

We  started  our  analysis  by  using  a  simple  Pekeris  model  to  identify  the 
order  of  the  bottom-scattered  waves.  Then  we  applied  a  method  of  critical 
reflections  and  refractions  used  by  Katz  and  Ewing6  to  derive  a  preliminary 
geoacoustic  model.  This  geoacoustic  model  was  refined  by  matching  the  disper¬ 
sion  of  the  theoretical  bottom-refracted  waves  to  the  dispersion  of  the 
observed  bottom-refracted  waves. 

In  Fig.  9  the  observed  and  theoretical  reduced  travel-time  curves  for 
shots  recorded  from  two  seismic  refraction  lines  during  the  FRAM  II  Expedition 
are  shown.  The  triangle  mark  in  Fig.  9  represents  the  observed  reduced  travel 
time  at  the  onset  of  each  bottom  bounce.  In  general,  the  theoretical  travel 
time  curves  are  in  good  agreement  with  the  observed  ones.  This  is  not  to  say 
that  this  simple  model,  a  two- layered  Pekeris  waveguide,  is  good  enough  to 
explain  the  propagation  of  the  bottom- in ter acting  signals  in  the  Arctic  ocean; 
the  theoretical  travel  time  curves  are  the  travel  time  curves  of  the  bottom- 
reflected  signals  only.  The  bottom-interacting  pulses  of  order  more  than  17 
in  Fig.  8  are  bottom-refracted  pulses  with  part  of  their  travel  paths  through 
the  bottom  sediment.  However,  the  good  agreement  between  the  theoretical  and 
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observed  travel-tine  curves  may  verify  the  identification  of  the  order  of  the 
bottom-scattered  waves. 

For  shots  recorded  at  close  ranges,  the  number  of  bottom  bounces  is  easy 
to  count,  but  for  those  shots  recorded  at  long  ranges,  the  order  of  bottom 
bounces  is  not  as  easily  identified.  The  reduced  travel  time  versus  distance 
plot  offers  us  a  convenient  method  with  which  to  identify  the  orders  of  the 
arriving  pulses.  Bottom-scattered  waves  up  to  twelfth  order  were  observed  for 
100  km  shots  at  FRAM  II  and  bottom-scattered  waves  of  more  than  thirty-second 
order  were  observed  for  340  km  shots  at  Camp  I. 

On  the  basis  of  bottom-reflections  at  the  critical  angle,  the  sound 
velocity  in  sediments  near  the  ocean  floor  can  be  calculated  form  the  maximum 
number  of  bottom  reflections8. 

Assuming  that  the  curvature  of  the  reflected  ray  in  the  water  column  is 
negligible  and  that  the  maximum  number  of  bottom  reflections  on  a  given  record 
is  determined  by  the  critical  angle  of  reflection,  the  approximate  value  of 
the  sound  velocity  beneath  the  sea  floor  calculated  by  this  method  is  1.63 
km/sec.  The  sound  velocity  at  the  bottom  of  the  Arctic  Ocean,  4.2  km  in 
depth,  is  about  1.520  km/sec.  The  estimated  sound  velocity  beneath  the  sea 
floor  of  1.63  km/sec  does  not  firmly  suggest  that  a  first-order  discontinuity 
exists  at  the  bottom  interface.  A  thin  layer  with  a  strong  velocity  gradient 
is  more  likely  to  represent  the  real  situation. 

Assuming  a  constant  positive  velocity  gradient  in  the  shallow  sediment, 
the  velocity  gradient  and  the  depth  of  deepest  penetration  can  be  estimated 
from  the  cycle  range,  travel  time  and  phase  velocities.  The  phase  velocities 
are  estimated  from  the  theoretical  travel  time  curves  shown  in  Fig.  9.  These 
computations  agree  well  with  a  model  having  a  velocity  gradient  of  2.0  s”1. 
The  thickness  of  the  shallow  sediment  as  defined  by  this  velocity  gradient  is 


about  100  meters.  The  same  method  can  be  applied  to  derive  the  velocity 
gradient  in  the  deeper  sediment  layer.  Repeating  the  same  procedure,  a 
velocity  gradient  of  1.0  s~^  is  derived. 

From  this  estimate  of  the  velocity  gradients  in  the  sediment  and  from  the 
analysis  of  the  pulse  shapes  of  the  bottom-interacting  signals,  a  geoacoustic 
model  for  the  wave  propagation  in  the  Pole  Abyssal  Plain  may  be  derived.  Fig. 
10  shows  the  velocity  profile  of  this  acoustic  model,  and  Fig.  11  shows  the 
density  profile.  This  velocity  profile  extends  the  Arctic  SOFAR  channel 
several  hundred  meters  into  the  bottom  sediment. 

Fig.  12  shows  an  observed  wave  train  recorded  at  a  range  of  295  km.  The 
signal  was  recorded  by  a  hydrophone  50  meters  deep. 

In  Fig.  12,  only  seven  of  the  bottom-interacting  pulses  in  the  computed 
wave  train  are  shown.  This  is  to  indicate  the  ability  of  the  Multiple 
Scattering  Pulse  Fast  Field  Program  to  perform  the  computation  for  a  single 
pulse,  such  as  those  bottom  multiples  of  order  14,  16,  18  and  20,  or  for 
several  pulses  at  one  time,  such  as  the  bottom  multiples  of  order  22  to  24. 
One  can  start  by  modeling  the  lower  order  bottom  bounces  (such  as  M18  or  M20, 
which  travel  through  the  shallow  sediment  layer)  in  order  to  derive  the 
structure  of  the  shallow  sediment,  then  extend  the  structure  deeper  to  model 
the  high-order  bottom  bounces.  This  procedure  shows  the  advantage  of  the 
Multiple  Scattering  Pulse  Fast  Field  Program.  With  the  MSPFFP,  one  can 
calculate  one  or  more  bottom-interacting  pulses,  testing  and  modifying  the 
model,  rather  than  compute  the  whole  wave  train.  This  saves  a  great  deal  of 
computing  time,  especially  when  a  wide  range  of  model  parameters  has  to  be 
tested. 

In  Fig.  13  the  observed  SOFAR  signal  from  Fig.  12  and  a  computed  SOFAR 
signal  are  compared.  In  Fig.  14  the  dispersion  profiles  of  the  observed  and 


computed  SOFAR  signal  are  shown.  These  dispersion  profiles  were  calculated  by 
the  conventional  "graphical  peak  and  trough"  method  and  were  checked  by  the 
spectral  filtration  method.  The  agreement  in  the  dispersion  profile  for  the 
first  mode  of  the  SOFAR  signal  is  excellent.  We  used  the  same  procedure  to 
compare  the  dispersion  profiles  of  the  computed  and  observed  bottom¬ 
interacting  signals. 

The  next  three  figures,  15,  16  and  17,  show  the  comparison  of  the 
observed  and  computed  waveforms  for  the  bottom  bounces  of  orders  20,  22  and 
24.  These  bottom-interacting  signals  are  the  same  as  those  shown  in 
Fig.  12.  The  dispersion  profiles  of  the  bottom-interacting  signals  having 
similar  cycle  ranges  of  15  km,  13.5  km  and  12.5  km  are  shown  in  the  next  three 
figures,  18,  19  and  20.  If  the  distance  traveled  in  the  sediment  is  long 
enough  such  that  constructive  interference  of  rays  in  the  sediment  has 
occurred,  then  the  bottom-interacting  signals  having  the  same  cycle  range  (or 
the  same  phase  velocity)  should  have  nearly  the  same  dispersion  profiles.  As 
shown  in  these  figures,  the  dispersion  profiles  of  the  computed  waveforms  and 
those  of  the  observed  waveforms  with  about  the  same  cycle  ranges  are  in  good 
agreement. 

The  agreement  between  the  dispersion  profiles  of  shot  signals  at  differ¬ 
ent  ranges  suggests  that  a  very  uniform  sediment  waveguide  exists  below  the 
Pole  Abyssal  Plain  in  the  Arctic  Ocean.  The  constructive  interference  of  rays 
having  a  travel  path  through  the  sediment  waveguide  is  responsible  for  the 
observed  frequency  dispersion  of  the  individual  bottom-interacting  pulses. 
The  agreement  of  the  dispersion  profiles  for  the  computed  and  observed  bottom¬ 
interacting  pulses  suggests  that  the  geoacoustic  model  is  a  fairly  good  model 
for  the  sediment  waveguide  below  the  Pole  Abyssal  Plain.  In  this  geoacoustic 
model,  a  positive  acoustic  wave  velocity  gradient  decreases  with  depth  from 


about  2  s  at  the  sediment  surface  to  about  1.0  s”1  at  about  400-500  meters 
depth  which  is,  in  general,  consistent  with  observations  in  the  other  oceans. 

The  existence  of  a  low  velocity  layer  in  the  superficial  layer  of  the 
deep-sea  sediment  has  been  observed,  but  because  the  frequency  content  of  the 
bottom-interacting  signals  is  less  than  40  Hz,  the  existence  of  a  very  thin 
low-velocity  layer  cannot  be  resolved  from  our  data.  However,  in  the 
numerical  simulations  of  high-order  bottom-interacting  signals,  a  layer  of 
low-velocity  sediment  15  m  thick  is  used  because  of  better  agreement  in  the 
high-frequency  portion  of  the  24th  bottom  bounce  than  that  computed  without  a 
low-velocity  layer.  But,  in  general,  there  are  no  differences  in  the 
frequency  dispersion  of  these  two  waveforms. 

It  can  be  concluded  that  a  uniform  waveguide  exists  below  the  Pole 
Abyssal  Plain.  The  consistency  of  the  observed  signals  suggests  that  the 
bottom  sediment  waveguide  may  be  considered  a  part  of  the  Arctic  SOFAR 
channel.  The  velocity  profile  in  the  surface  channel  is  responsible  for  the 
frequency  dispersion  of  the  Arctic  SOFAR  signals,  and  the  positive  velocity 
gradient  in  the  bottom  sediment  is  responsible  for  the  dispersion  of  the 
bottom-refracted  pulses.  For  long-range  sound  propagation  over  the  Arctic 
abyssal  plains,  at  least  up  to  400  km  in  range,  the  bottom  sediment  structure 
several  hundred  meters  below  the  sea  floor  is  important. 

The  geoacoustic  model  derived  from  Pole  Abyssal  Plain  was  used  to  model 
the  bottom-interacting  signals  recorded  in  both  the  Mendeleyev  Abyssal  Plain 
and  the  Barents  Abyssal  Plain.  Good  agreement  in  wave  forms  of  the  observed 
and  computed  bottom-interacting  signals  suggests  that  similar  bottom  structure 
may  exist  below  these  abyssal  plains  of  the  Arctic  Ocean,  at  least  for  the 
upper  few  hundred  meters.  The  geoacoustic  model  derived  from  the  Pole  Abyssal 
Plain  is  also  similar  to  that  proposed  by  Hamilton9.  Although  the  depth- 


dependent  attenuation  coefficients  could  not  be  resolved  accurately  from  our 
data/  our  average  value  of  0.8  dB/km  at  10  Hz  is  in  the  range  of  the 
attenuation  coefficients  given  by  Hamilton. 

The  uniformity  of  group-velocity  dispersion  of  waves  corresponding  to  the 
first  normal  mode  in  the  frequency  band  8  to  30  Hz  suggested  the  possibility 
of  achieving  significant  signal  gain  by  pulse  compression  in  the  Arctic  SOFAR 
channel10.  Numerical  simulations  shown  in  Fig.  21  were  accomplished  by  propa¬ 
gating  time  reversed  predistorted  wave  forms  matched  to  the  dispersion  of  the 
channel  at  specified  ranges.  The  similarity  of  the  frequency  dispersion 
between  the  Arctic  SOFAR  signals  and  the  bottom-interacting  signals  over  the 
abyssal  plains  suggests  the  possibility  of  interference  with  a  pulse 
compression  scheme. 

A  computed  shot  signal,  including  all  the  bottom-interacting  pulses,  is 
shown  in  Fig.  22.  Fig.  23  shows  a  simulation  of  pulse  compression  by  cross- 
correlating  the  SOFAR  signal  and  the  total  response  including  all  the  bottom 
interactions.  Since  the  dispersion  of  the  bottom-refracted  pulses  is  similar 
to  that  of  the  SOFAR  signals,  pulse  compression  of  the  bottom- interacting 
response  by  the  SOFAR  channel  response  is  expected.  It  can  be  seen  in  Fig.  23 
that  after  a  travel  time  of  230  seconds,  the  bottom-refracted  pulses  are 
indeed  compressed,  especially  those  pulses  (M20,  M21)  which  deeply  penetrate 
the  shallow  sediment  waveguide  (with  a  strong  velocity  gradient)  but  do  not 
travel  into  the  deeper  sediment  layer  (with  a  weaker  velocity  gradient) .  Fig. 
24  shows  the  details  of  the  pulse  compression  of  the  bottom- interacting 
signals  of  order  18  to  24.  The  amplitude  scale  in  Fig.  24  is  three  times  that 
in  Fig.  23.  The  most  compressed  pulse  corresponds  to  the  bottom  bounce  of 
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The  effect  of  the  bottom-interacting  signals  on  pulse  compression  has 
been  shown  by  computer  simulation.  The  pulse  compression  of  the  SOFAR  signal 
is  observed,  but  many  more  sidelobes  occur.  The  sidelobes  generate  more  noise 
for  the  pulse  compression  scheme  to  deal  with.  The  dispersed  bottom-refracted 
pulses  are  indeed  compressed  during  the  simulation,  but  the  aaplitude  of  the 
compressed  bottom  refracted  pulses  is  much  smaller  than  that  of  the  compressed 
SOFAR  signal. 
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Fig.  1  Velocity  profile  in  Eurasian  Basin  calculated  frost  hydrographic 

data. 

Fig.  2  Ray  paths  confuted  at  one-degree  interval. 

Fig.  3  Waveform  of  typical  SOFAR  signal  measured  in  Eurasian  Basin. 

Fig.  4  Sound  spectrogram  of  typical  SOFAR  signal  measured  in  Eurasian 

Basin.  The  signal  is  shown  in  Fig.  3.  Ml  first  mode.  M2  second 
mode. 

Fig.  5  Comparison  of  measured  and  computed  group-velocity  dispersion  for 

the  first  two  normal  modes. 

Fig.  6  Measured  waveform  which  shows  strong  bottom-interacting 

signals.  Range  333.5  km. 

Fig.  7  Waveforms  of  two  bottom-interacting  signals  of  Fig.  6. 

Fig.  8  Waveform  and  corresponding  sound  spectrogram  of  typical  signal 

measured  over  flat  bottom  of  the  Pole  Abyssal  Plain.  Mean  water 
depth  4200  m.  Range  317  km.  Charge  size  55  pounds  TNT. 
Dispersion  profiles  for  the  17th  and  29th  mode  (M17,  M29)  are 
shown  on  the  spectrogram. 

Fig.  9  Plot  of  the  reduced  travel  time  curves  of  shot  records  from  two 

seismic  ref ration  lines  of  the  FRAM  II  Expedition.  The  shots  in 
the  0  -  80  km  range  were  recorded  at  FRAM  II  and  the  shots  in  the 
200  -  340  km  range  were  recorded  at  Camp  I.  The  continuous  lines 
are  theoretical  reduced  travel  time  curves  computed  from  a  two¬ 
layered  Pekeris  waveguide. 

Fig.  10  The  velocity  profile  of  the  geoacoustic  model. 


Fig.  11  The  density  profile  of  the  geoacoustic  model. 


Fig.  12  Measured  and  computed  waveforms  which  show  strong  bottom¬ 

interacting  signals.  Range  295  km.  Charge  size  880  pounds 
TNT.  The  computed  waveforms  were  calculated  by  the  Multiple 
Scattering  Pulse  Fast  Field  Program  (MSPFFP) . 

Fig.  13  Measured  and  computed  SOFAR  signal. 

Fig.  14  Dispersion  profile  corresponding  to  the  first  mode  of  the 

observed  and  computed  SOFAR  signal  in  Fig.  13. 

Fig.  15  Observed  and  computed  waveforms  of  20-th  bottom  bounce  (M20)  in 
Fig.  12. 


Fig.  1$  Observed  and  computed  waveforms  of  22-nd  bottom  bounce  (M22)  in 

Fig.  12. 

Fig.  17  Observed  and  computed  waveforms  of  24-th  bottom  bounce  (M24)  in 

Fig.  12. 

Fig.  18  Dispersion  profiles  for  bottom-interacting  signals  having  a  cycle 

range  about  15  km.  The  dispersion  profiles  of  waveforms  shown  in 
Fig.  15  are  included. 

Fig.  19  Dispersion  profiles  for  bottom- interacting  signals  having  a  cycle 
range  about  13.5  km.  The  dispersion  profiles  of  waveforms  shown 
in  Fig.  16  are  included. 

Fig.  20  Dispersion  profiles  for  bottom-interacting  signals  having  a  cycle 

range  about  12.5  km.  The  dispersion  profiles  of  waveforms  shown 
in  Fig.  17  are  included. 

Fig.  21  Computed  SOFAR  signals  at  ranges  from  300  to  900  km  and  pulse 

compression  of  SOFAR  signals.  Source  signal  for  pulse  compres¬ 
sion  is  the  time  reverse  of  second  SOFAR  signal  from  top. 

Fig.  22  Computed  waveform  at  295  km,  including  all  the  bottom-interacting 

signals.  Geoacoustic  model  of  the  bottom  is  shown  in  Figs.  10 
and  11.  The  water  depth  is  4200  m. 

Fig.  23  Pulse  compression  at  a  range  of  295  km.  All  the  bottom¬ 

interacting  signals  are  included. 

Fig.  24  The  details  of  the  pulse  compression  shown  in  Fig.  23  of  bottom¬ 

interacting  signals  of  order  18  to  24. 
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Introduction 


This  volume  is  being  published  to  present  and  describe  for  use 
the  current  computer  work  of  Dr.  H.  Kutschale. 

Dr.  Kutschale  is  presently  located  at  the  Lamont-Doherty 
Geological  Observatory  of  Columbia  University  in  Palisades>  New  York. 

The  work  described  here  is  a  continuation  of  many  years  of  work 
by  Dr.  Kutschale  in  the  field  of  underwater  acoustics  using  the  Fast 
Field  Program  (  FFP  ).  Investigators  in  the  Naval  community  who  have 
had  previous  exposure  to  Dr.  Kutschale 's  work  and  who  have  made  use  of 
his  work  will  find  the  current  work  to  be  easily  incorporated  into  their 
own  efforts  due  to  its  similarity  of  input  with  the  older  versions  of 
these  programs.  All  will  find  it  easily  understood  due  to  the  step  by 
step  style  of  documentation  in  each  of  the  user's  manuals.  In  most  cases 
the  user  groups  have  received  the  current  programs  along  with  many  of  the 
test  cases  described  in  this  volume. 

TLOSS11CC*  the  first  program  listed*  is  the  'basic'  program  and 
is  an  improvement  over  previous  FFP  programs  due  to  its  ability  to  handle 
sound  velocity  profiles  of  a  more  general  nature  than  before.  Equal 
velocity  layers  are  no  longer  required.  Continuously  varying  velocities 
are  handled  using  the  Airy  function.  Also*  double  precision  variables 
enable  the  program  to  remain  stable  when  the  input  is  far  from  the  norm* 
though  not  to  the  point  of  ridiculousness. 

TL0SS22CC  is  a  step  up  with  all  the  features  of  TLOSS11CC  along 
with  the  ability  to  print  out  the  results  of  the  transmission  loss  due 
to  one  particular  chosen  bottom  bounce. 

TLOSSIOHH  is  also  an  enhanced  version  of  TL0SS11CC.  This  program 
adds  the  ability  to  use  the  Normal  Mode  method  to  all  the  features  of 
TLOSS1 ICC. 

TL0SS17HH  starts  with  TLOSSIOHH  as  a  foundation  and  adds  the 
ability  to  calculate  the  branch  integral  with  the  Normal  Mode  method. 

TL0SS51HH  gives  the  user  all  of  the  power  of  TL0SS17HH  plus  the 
option  of  subdividing  the  input  sound  velocity  profile  in  order  to  do 
step  by  step  integration.  Also*  calculations  of  the  eigen  functions  are 
available  using  the  Normal  Mode  method  and  these  eigen  function  values 
along  with  the  integrand  values  are  put  out  to  disk  for  further  use  by 
other  programs. 

IMPULS11B  is  the  basic  impulse  program*  corresponding  in  function 
to  TL0SS11CC. 

IMPULS22C  is  the  impulse  program  that  corresponds  to  TL0S522CC 
and  has  the  same  BOUNCE  equations. 

IMPULS11G  is  the  impulse  program  that  corresponds  to  TLOSSIOHH. 

The  Normal  Mode  method  is  an  option  in  this  program. 

In  addition  to  these  eight  FFP  programs  Dr.  Kutschale  is  also 
currently  using  two  Parabolic  Equation  programs.  One  program  is  a 
transmission  loss  program*  PELOsSX*  and  one  program  is  an  impulse 
program*  PEPULSEX*  both  of  which  were  developed  by  Dr.  Kutschale  in 
conjunction  with  Dr.  F.  Tappert. 

These  ten  programs  are  put  forth  here  TO  ENABLE  A  PERSON  TO  USE 
THESE  PROGRAMS*  and  for  no  other  reason.  One  must  delve  deeper  than 
this  volume  in  order  to  revise  any  of  these  programs  to  meet  a  special 
need  not  as  yet  addressed.  What  is  included  for  each  program  is: 

1.  A  User's  Manual  describing  input  and  special  features 

2.  The  input  for  several  selected  test  cases 

3.  The  output  for  these  same  selected  test  cases 


TLQSS1 ICC  USER'S  MANUAL 


Introduction 

TLOSS11CC*  our  streamlined  transmission  loss  (  hereafter  abbreviated 
as  TL  )  program,  was  completed  in  October  of  1932. 

TL03S11CC  has  these  features  and  options  for  the  user: 

1.  The  FFP  computational  method 

3.  A  comprehensive  selection  of  modelling  parameters. 

These  include: 

a.  A  detailed  sound  velocity  profile 

b.  Choice  of  source  depth 

c.  Choice  of  detector  depth 

d.  Choice  of  frequency 

e.  Choice  of  wave  number  integrand  limits 

f.  Optional  beaming*  as  it  is  affected  by 
integrand  tapering  (  in  the  FFP  mode  ) 

g.  Optional  surface  roughness 

n.  Choice  of  either  fluid  or  vacuum  upper  half 
space 

i.  Possibility  of  the  source  and  detector  in 
the  fluid  upper  half  space 

j.  Optional  internal  roughness 

4.  Output.  This  includes: 

a.  A  printout  of  all  input 

b.  The  TL  plot 

c.  The  wave  integrand  plot  of  the  FFP;  either 
short  or  long. 

Note:  All  short  plots  except  the  wave  integrands 
are  repeated*  the  repeated  plot  being  a  smoothed 
copy  of  the  first  and  having  the  same  bounds  as 
the  first. 

5.  Additional  Optional  Tables  and  Plots.  These  include: 

a.  The  TL  plot  in  nautical  miles*  if  desired  and 
if  the  range  is  sufficiently  long.  The  original 
TL  plot  is  in  kilometers. 

6.  Numerical  output  of  the  TL  for  further  processing  or  plotting 
Formerly*  this  was  called  punched  output*  since  it  went  out 

to  the  card  punch. 


v.vra 


TL0SS1 ICC  USER'S  MANUAL 
Input  Variables  and  Formats 

Program  input  for  TLOSS11CC  includes  both  the  physical  model 
and  program  control.  This  section  of  the  manual  will  list 
and  explain  each  line  of  the  input. 


Entry 
MMAXF, 
(  213 


1/  one 
NCARD, 


line 

NREAO 


variable  names 
format 


MMAXF  is  a  control  variable  which  is  used  to  read 
in  the  correct  number  of  points  in  the  sound  velocity 
profile.  (  see  2.  -  next  Entry  ) 

NCARD  is  a  control  variable  that*  when  set  non-zero 
allows  the  program  to  put  out  TL  data  on  unit  7  in 
a  form  which  is  readable  as  input  bu  another  program. 

A  zero  value  will  prevent  any  output  on  unit  7.  This 
is  normally  used  to  make  graphic  plots.  <  Originally* 
on  the  IBM  360/91*  unit  7  was  the  punched  card  unit* 
hence  the  name  NCARD.  > 


NREAD  is  a  control  variable  that*  when  set  non-zero* 
directs  the  program  to  interpret  Entry  2  data  (  see  below  > 
as  thickness  input  instead  of  the  usual  depth  input. 


2.  Entry  2*  MMAXF  lines  (  max.  of  49  lines  ) 
((  V<M, N) <  N=1 #  7  ).  M=l,  MMAXF  ) 

(  7D10.  0  ) 


variable  name 
format 


V  is  an  array  that  contain 

For  each  point  in  the  profile* 

(  one  point  per  line  )>  the  pro 

a.  a  depth  value*  in  meters*  e 
in  which  case  the  program  ex 

b.  a  sound  (  pressure  )  veloci 

c.  a  sound  (  shear  >  veloci 
if  applicable 

d.  a  corresponding  attenuation 

e.  a  corresponding  attenuation 
if  applicable 

f.  a  density  value*  in  grams  p 

g.  a  roughness  value*  if  appli 

Note:  No  depth  value  can 

depth  value.  The  program  will 


s  the  sound  velocity  profile, 
and  on  each  line* 
gram  expects: 

xcept  when  NREAD  is  set  non-zero* 
peccs  a  thickness  value,  in  meters 
ty  value*  in  meters  per  second 
ty  value*  in  meters  per  second. 

value  for  pressure 
value  for  shear. 


er  c.  c.  *  or  specific  density 
cable 

be  less  than  a  previous 
stop  if  this  is  so. 


3.  Entry  3*  one  line 

KSF.  LRF,  MM,  LS.  KKK.  I  BEAM,  NW,  FREQ*  PO.  C.  CMAX 
(  715,  4D10.  0  ) 

KSF  is  a  pointer  to  the  source  depth  value. 

It  paints  to  one  of  the  points  in  the  sound  velocity 
profile*  the  one  which  is  the  source  depth  value. 

LRF  is  a  pointer  to  the  detector  (  or  receiver  > 


variable  names 
format 


depth  val 


is  a  poin 
ue.  It's 


function  is  similar  to  KSF. 


MM  is  a  TL  plot  pruning  factor*  with  a  minimum 
acceptable  value  of  2.  MM  is  applicable  to  FFP  only. 

If  an  MM  value  of  N  is  read  in*  then  the  first  1/N 
of  the  plot  will  be  plotted.  This  is  useful  either  to 
get  a  more  detailed  look  at  the  output  or  to  prune 
away  that  part  of  the  plot  affected  by  range  aliasing. 

LS  is  the  number  of  points  used  in  the  wave 
integration.  It  also  affects  the  range  of  the  TL  plot, 
for  the  FFP  method. 

KKK  is  a  control  variable  that  indicates: 

a.  If  zero,  this  is  the  last  case. 

b.  If  negative,  the  next  case  starts  with  Entry  1* 

(  new  profile  ). 

c.  If  positive,  the  next  case  starts  with  Entry  3. 

IBEAM  is  a  control  variable  that*  when  set  non¬ 
zero*  requires  the  program  to  taper  the  wave  integrand 
at  each  end. 


NW  is  a  flag,  that  whan  sat  non-zero.  diracts  tha  program 
to  put  out  a  printar  plot  of  tha  wave  intagrand. 

FREQ  is  the  frequency  for  this  casa.  in  Hartz. 

PO  is  tha  sourca  laval.  For  TL  PO  ■  1.0 

C  is  the  lower  bound  sound  velocity  in  tha  wave 
integration,  in  maters  par  second. 

CMAX  is  tha  upper  bound  sound  velocity  in  tha  wave 
integration,  in  maters  par  second.  If  a  zero  is  entered 
hare,  tha  wave  number  integration  will  start  from  zero. 

Entry  4.  one  line 


Entry  4.  one  line 

REFLTO.  ALPHAO,  ATTENO.  RHOO.  H,  INAUT.  LONG 
(  5D10.  0.  215  ) 


variable  name 
format 


REFLTO  is  tha  RMS  surface  roughness  value,  in  meters. 
If  this  value  is  non-zero,  tha  program  explicitly  sets 
tha  next  three  values  to  zero. 

ALPHAO  is.  if  present,  the  sound  velocity  in  the 
fluid  upper  half  space,  in  meters  per  second.  This  is 
the  variable  that  determines  whether  the  present  case  has 
either  a  fluid  or  vacuum  upper  half  space. 

ATTENO  is.  if  present,  the  attenuation  in  the 
fluid  upper  half  space. 

RHOO  is.  if  present,  the  density  of  the  fluid 
upper  half  space,  in  grams  per  c. c. 

H  is.  if  present,  the  distance  between  the  source 
and  detector  in  the  fluid  upper  half  space.  To  use  this 
option: 

a.  The  first  two  sound  velocity  depths  must  be  zero. 

b.  There  must  be  some  variation  in  the  pressure 
sound  velocity  profile. 

c.  There  must  not  be  any  shear  sound  velocity. 

d.  ALPHAO  and  H  must  be  present. 

INAUT  is  a  control  variable  that,  when  set  non¬ 
zero.  allows  the  plotting  of  the  TL  plot  in  nautical 
miles  if  the  range  is  greater  than  100  kilometers. 


.lometers. 


LONG  is  a  control  variable  that,  when  set  non¬ 
zero.  switches  the  wave  integrand  plot  from  short  to 
long.  This  is  used  when  it  is  necessary  to  examine 
each  point  of  the  wave  integrand. 
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TLQSSl ICC  USER'S  MANUAL 
C.  Operating  Procedures 

The  primary  computer  at  Lamont-Doherty  Geological  Observatory 
is  a  VAX  11/780.  which  runs  under  the  VMS  operating  system. 

Listed  below  are  four  distinctly  different  operations  that 
are  used  to  process  TLOSS11CC. 

1.  Compiling-Linking 

a.  Compiling  (  PORT  command  > 

Compiling  the  TLOSS11CC  program  on  the  VAX  11/780 
has  always  been  done  with  the  following  switches 
on: 

1.  /CHECK  Forbids  array  pointer  overflow. 

2.  /N0I4  Integer  variables  are  set  to  2  bytes. 

3.  /N0F77  Where  the  compiler  needs  to  choose 
between  Fortran  77  and  Fortran  4+.  the  latter 
is  chosen. 

4.  /LIST  A  compiler  listing  is  made  and 
stored  on  the  disk  as  TLOSSl ICC.  LIS. 

b.  Linking  (  LINK  command  > 

Linking  is  straight-forward.  No  extra  routines 
are  linked  with  TLQSS11CC. 

2.  Running 

TLoSSl ICC  is  normally  run  through  a  batch  stream* 
in  particular  the  SLOW  batch  stream*  since  this  is  the 
cheapest  route*  especially  at  night.  However,  there  is 
another  batch  stream  with  higher  priority  (  and  higher 
cost  )  that  can  be  used  for  quicker  results,  the 
SYS4BATCH  batch  stream.  In  order  to  use  either  of  these 
two  batch  streams: 


a. 


The  executable  -file  <  TL0SS11CC.EXE  )  must  be  ready 
to  use.  The  LINK  operation  satisfies  this  requirement. 


b.  The  data  file  must  also  be  ready  to 
normally  prepared  by  the  user  using  the 
<  and  the  TL0SS1 ICC  USER'S  MANUAL,  if 
at  a  terminal. 


use.  This 
editor, 
necessary  > 


is 


c.  Submit  the  command  file,  TL0SS11CC.COM*  to  the 
batch  stream  (  or  batch  queue  >  using  a  command 
such  as  one  of  these: 

SUBMIT/NAME-LEC010CT/QUEUE-SLQW/AFTER-20: 00  TLQSSl ICC 
SLOW  queue.  Job  starts  at  8:00  P.  M. 

SUBMIT/NAME-LEC010CT/QUEUE-SLQW  TLQSSl ICC 
SLOW  queue.  Job  starts  now 

SUBMIT/NAME-LEC010CT  TL0SS1 ICC 

SYS4BATCH  queue  (  the  'regular'  queue  ). 

Job  starts  now 

d.  The  following  command  file*  also  prepared  using 
the  editor,  is  a  typical  file  submitted*  using  one 
of  the  previous  submit  commands. 

4 ASSIGN  LEC010CT  SYS4PRINT  <1) 

4 ASSIGN  LEC010CT  SYS40UTPUT  (2) 

4 ASSIGN  DRO: LEC010CT.  PLT  F0R007  (3) 

4ASSIGN  DRO: LEC010CT. OUT  F0R006  (4) 

4ASSIGN  DRO: STANDARD.  DAE  F0R005  (S) 

4TYPE  TL0SS1 ICC  < A) 

4C0PY  DRO:  STANDARD.  DAE  DRO:  LEC 01  OCT.  DAT  <7> 

4RUN  TL0SS1 ICC  <8> 

This  run  is  called  the  LEC010CT  run.  'LEC'  stands  for 
LOSS  ELEVEN  C.  '010CT'  indicates  that  this  is  the 
first  run  in  October. 
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log 

The 

The 

The 

The 

The 

The 


first  two  lines  in  this  command  file  put  the 
file  on  the  disk  instead  of  the  printer, 
third  line  puts  unit  7  output  into  file  LEC010CT.  PLT. 
fourth  line  puts  unit  6  output  into  file  LEC010CT.  OUT. 
fifth  line  draws  unit  9  input  from  STANDARD.  DAE. 
sixth  line  copies  TLOSS11CC.  LIS  to  LEC010CT.  LIS. 
seventh  line  copies  STANDARD.  DAE  to  LECOIOCT.  DAT. 
eighth  line  runs  program  TLOSS11CC. 


A  command  file  such  as  this  yields  a  file  group; 
each  file  having  the  same  name  but  a  different 
extension.  These  files  are  then  treated 
collectively.  They  are  printed  together  and  stored 
on  tape  together.  Such  a  group  is*  for  example: 


LECOIOCT.  DAT 
LECOIOCT.  LIS 
LECOIOCT.  LOO 
LECOIOCT.  OUT 
LECOIOCT.  PLT 


a  copy  of  the  data  file 
the  program  listing  file 
the  log  file 
the  output  file 
the  numerical  output  file 


The 

run 


user  normally 
by  typing  the 


gets  a  complete 
command: 


printout  of  the 


PRINT  LECOIOCT.*;* 


on  the  terminal. 


3. 


Copying  Runs  to  Tape  (  Archiving  > 

Periodically,  it  is  necessary  to  clear  the  old  runs 
off  the  disk  and  put  them  on  tape.  This  is  a  very  simple 
operation  and  proceeds  as  follows: 


a.  Initialize  a  tape.  Give  it  a  name  tnat  is  somehow 
indicative  of  the  present  month.  For  example: 


INITIALIZE/DENSITY-1600  MT1:  OCTNOV 


This  tape  would  have  runs  made  in  October  and  November. 
This  step  is  not  necessary  if  the  user  is  adding  runs 
to  a  tape  that  is  already  partly  filled.  In  such  a 
case,  be  careful  not  to  initialize! 


Mount  the  tape  with  the  MOUNT  command. 
MOUNT  MT1 :  OCTNOV 


This  is  not  the  same  as  physically  mounting  the  tape 
on  the  tape  drive,  which  is  assumed  to  have  taken 
place  previously. 


c.  Copy  all  the  run  files  from  disk  to  tape  with  thi 
COPY  command.  For  example  in  this  case: 


COPY  *OCT.  *;  *  MT1 : 
COPY  *NOV.  *;  *  MT1 : 


would  copy  all  the  October  and  all  the  November  runs  onto 
the  tape  mounted  on  tape  drive  1. 


d.  For  future  reference, 
names  of  the  files  on  the 
as  follows: 


it  is 
tape. 


good  to  get  a  list  of 
This  is  accomplished 


DIRECTORY/PRINT  MT1 : 


command  does  not 
user  should  write 


The  printout  resulting  from  this 
have  the  tape  name  on  it.  so  the 
this  on  immediately. 

If.  for  some  reason,  this  printout  is  lost,  th 
can  be  remounted  and  the  directory  can  be  taken 
again.  But  this  is  no  reason  to  be  careless  with 
tne  original  directory  printout. 


e  tape 


e.  Dismount  the  tape  with  the  dismount  command: 

DISMOUNT  MT1 : 

f .  Delete  the  files  on  the  disk  that  have  just  been 
copied  to  the  tape.  In  this  case: 

DELETE  *OCT.  *;  * 

DELETE  *NOV.  *;  * 

4.  Distributing  the  Program  by  Tape 

The  program  and  test  cases  are  distributed  using 
tapes  with  the  same  format  as  the  archival  tapes.  Tnis 
is  the  Files-11  format-  peculiar  to  Digital  Equipment 
Corporation.  The  tapes  from  this  system  have  a  density 
of  1600  bpi-  and  a  blocksize  of  204$  bytes. 


If  any  further  questions  arise-  feel  free  to  call: 

Dr.  H.  Kutschale  (914)  3S9-2900  x382 

Mr.  L.  Carroll  (914)  359-2900  x221 
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06  1 
0.  00 
50.  00 
100.  0 
100.  0 
500.  0 
500.  0 
2 

0.  0 
00. 


TEST  CASE  DATA 


1500.  0 
1500.  0 
1500.  0 
3000.  0 
3000.  0 
3000.  0 
02  1024 
1500.  0 
00.  0 


1400.  0 
1400.  0 
1400.  0 
00  000 
.  0000150 


•  0000150 
. 0000150 

•  0000150 

. 00007700 
.  00007700 
.  00007700 
1  020.  0 


1.  00 


DATA  FOR  TLQSS1 ICC 
1.  00 
1.  00 
1.  00 

.  0001600  2.  50 
0001600  2.  50 
.0001600  2.  50 
•  0  10000.  100< 


0000.  0 


1000.  0  0000  0 
1  00  1 


A. 


Introduction 


TL0SS22CC.  one  of  our  streamlined  transmission  loss  < 
abbreviated  as  TL  )  programs,  uas  completed  in  March  of 

TL0SS22CC  has  these  features  and  options  for  the  user: 

1.  The  FFP  computational  method 

a.  The  complete,  unabridged  version 

b.  An  abbreviated  version  to  calculate  the 
one  particular  number  of  bounces 

3.  A  comprehensive  selection  of  modelling  paramete 
These  include: 

a.  A  detailed  sound  velocity  profile 

b.  Choice  of  source  depth 

c.  Choice  of  detector  depth 

d.  Choice  of  frequency 

e.  Choice  of  uave  number  integrand  limits 

f.  Optional  beaming,  as  it  is  affected  by 
integrand  tapering  (  in  the  FFP  mode  ) 

g.  Optional  surface  roughness 

n.  Choice  of  either  fluid  or  vacuum  upper 
spac  e 

i.  Possibility  of  the  source  and  detector 
the  fluid  upper  half  space 

j.  Optional  internal  roughness 

4.  Output.  This  includes: 

a.  A  printout  of  all  input 

b.  The  TL  plot 

c.  The  uave  integrand  plot  of  the  FFP;  eit 
short  or  long. 

Note:  All  short  plots  except  the  uave  integ 
are  repeated,  the  repeated  plot  being  a  smo 
copy  of  the  first  and  having  the  same  bound 
the  first. 

5.  Additional  Optional  Tables  and  Plots.  These  in 

a.  The  TL  plot  in  nautical  miles,  if  desir 
if  the  range  is  sufficiently  long.  The 
TL  plot  is  in  kilometers. 

6.  Numerical  output  of  the  TL  for  further  processi 
Formerly,  this  uas  called  punched  output,  since 
to  the  card  punch. 
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Input  Variables  and  Formats 


Program  input  for  TL0SS22CC  includes  both  the  physical  model 
and  program  control.  This  section  of  the  manual  will  list 
and  explain  each  line  of  the  input. 


Entry 
MMAXF, 
{  213 


1,  one 
NCARD, 


line 

NREAO 


variable  names 
format 


MMAXF  is  a  control  variable  which  is  used  to  read 
in  the  correct  number  of  points  in  the  sound  velocity 
profile.  (  see  2.  -  next  Entry  ) 


Entry  ) 


NCARD  is  a  control  variable  that,  when  set  non-zero 
allows  the  program  to  put  out  TL  data  on  unit  7  in 
a  form  which  is  readable  as  input  bu  another  program. 

A  zero  value  will  prevent  any  output  on  unit  7.  This 
is  normally  used  to  make  graphic  plots.  (  Originally/ 
on  the  IBM  360/91/  unit  7  was  the  punched  card  unit/ 
hence  the  name  NCARD.  > 


NREAD  is  a  control  variable  that/  when  set  non-zerOi 
directs  the  program  to  interpret  Entry  2  data  (  see  below  > 
as  thickness  input  instead  of  the  usual  depth  input. 

Entry  2/  MMAXF  lines  (  max.  of  49  lines  ) 

((  V(M,  N) ,  N*l,7  )/  M*1,MMAX*  )  variab 

(  7D10.  0  )  format 


variable  name 
format 


V  is  an  array  that  contain 

For  each  point  in  the  profile/ 

(  one  point  per  line  )z  the  pro 

a.  a  depth  value/  in  meters/  e 
in  which  case  the  program  ex 

b.  a  sound  (  pressure  )  veloci 

c.  a  sound  (  shear  >  veloci 
if  applicable 

d.  a  corresponding  attenuation 

e.  a  corresponding  attenuation 
if  applicable 

f.  a  density  value/  in  grams  p 

g.  a  roughness  value/  if  appli 

Note:  No  depth  value  can 

depth  value.  The  program  will 


s  the  sound  velocity  profile, 
and  on  each  line/ 
gram  expects: 

xcept  when  NREAD  is  set  non-zero# 
pects  a  thickness  value>  in  meters 
ty  value/  in  meters  per  second 
ty  value/  ip  meters  per  second# 

value  for  pressure 
value  for  shear. 

er  c.  c.  #  or  specific  density 
cable 

be  less  than  a  previous 
stop  if  this  is  so. 


3.  Entry  3.  one  line 

KSF,  LRF,  MM,  LS,  KKK,  IBEAM#  NW#  FREQ,  PO,  C,  CMAX 
(  715,  4D10.  0  ) 

KSF  is  a  pointer  to  the  source  depth  value. 

It  points  to  one  of  the  points  in  the  sound  velocity 
profile,  the  one  which  is  the  source  depth  value. 


variable  names 
format 


LRF  is 
depth  value 


s  a  pi 

e.  I< 


function 


similar  to  KSF. 


MM  is  a  TL  plot  pruning  factor,  with  a  minimum 
acceptable  value  of  2.  MM  is  applicable  to  FFP  only. 

If  an  MM  value  of  N  is  read  in,  then  the  first  1/N 
of  the  plot  will  be  plotted.  This  is  useful  either  to 
get  a  more  detailed  look  at  the  output  or  to  prune 
away  that  part  of  the  plot  affected  by  range  aliasing. 

LS  is  the  number  of  points  used  in  the  wave 
integration.  It  also  affects  the  range  of  the  TL  plot, 
for  the  FFP  method. 

KKK  is  a  control  variable  that  indicates: 

a.  If  zero,  this  is  the  last  case. 

b.  If  negative,  the  next  case  starts  with  Entry  1, 

<  new  profile  ). 

c.  If  positive,  the  next  case  starts  with  Entry  3. 

IBEAM  is  a  control  variable  that,  whan  set  non¬ 
zero,  requires  the  program  to  taper  the  wave  integrand 
at  each  end. 
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NW  is  a  flag*  that  whan  sat  non-zero.  directs  tha  program 
to  put  out  a  printar  plot  of  tha  wave  integrand. 

FREQ  is  tha  frequency  for  this  case*  in  Hertz. 

P 0  is  tha  sourca  laval.  For  TL  PO  ■  1.0 

C  is  tha  louiar  bound  sound  valocity  in  tha  uiava 
integration,  in  metars  par  sacond. 

CMAX  is  the  upper  bound  sound  valocity  in  tha  wave 
integration,  in  maters  par  sacond.  If  a  zero  is  entered 
here,  tha  wave  number  integration  will  start  from  zero. 

Entry  4.  one  line 

REFLTo.  ALPHAO.  ATTENO.  RHOO.  H.  INAUT.  LONG.  BOUNCE  variable  name 
(  5010.  0.  315  )  format 

REFLTO  is  tha  RMS  surface  roughness  value,  in  maters. 

If  this  value  is  non-zero,  tha  program  explicitly  sets 
the  next  three  values  to  zero. 

ALPHAO  is.  if  present,  tha  sound  velocity  in  tha 
fluid  upper  half  spaca.  in  maters  par  sacond.  This  is 
the  variable  that  datarminas  whether  tha  present  case  has 
either  a  fluid  or  vacuum  upper  half  spaca.  ! 

ATTENO  is.  if  present,  tha  attenuation  in  tha  I 

fluid  upper  half  spaca. 

RHOO  is.  if  present,  tha  density  of  tha  fluid  I 

upper  half  spaca.  in  grams  par  c.  c.  ,1 

H  is.  if  present,  tha  distance  between  tha  sourca 
and  detector  in  tha  fluid  upper  half  spaca.  To  use  this 
option: 

a.  Tha  first  two  sound  valocity  depths  must  be  zero. 

b.  Thera  must  ba  some  variation  in  tha  pressure 
sound  velocity  profile. 

c.  Thera  must  not  ba  any  shear  sound  valocity. 

d.  ALPHAO  and  H  must  ba  present. 

INAUT  is  a  control  variable  that,  whan  set  non¬ 
zero.  allows  tha  plotting  of  tha  TL  plot  in  nautical 
miles  if  the  range  is  greater  than  100  kilometers. 

LONG  is  a  control  variable  that,  when  set  non¬ 
zero.  switches  the  wave  integrand  plot  from  short  to 
long.  This  is  used  when  it  is  necessary  to  examine  i 

each  point  of  the  wave  integrand. 

BOUNCE  is  a  control  variable  that,  when  set  non-zero, 
directs  the  program  to  calculate  only  the  effects  of  a  certain 
number  of  bounces.  BOUNCE  being  that  number.  This  is  somewhat 
similar  to  the  'H  variable'  output  and  should  not  be  used  in 
conjunction  with  that  feature. 
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C.  Operating  Procedures 

The  primary  computer  at  Lamont-Doherty  Geological  Observatory 
is  a  VAX  11/780*  which  runs  under  the  VMS  operating  system. 

Listed  below  are  four  distinctly  different  operations  that 
are  used  to  process  TL0SS22CC. 

1.  Compiling-Linking 

a.  Compiling  (  FORT  command  ) 

Compiling  the  TL0SS22CC  program  on  the  VAX  11/780 
has  always  been  done  with  the  following  switches 
on: 

1.  /CHECK  Forbids  array  pointer  overflow. 

2.  /NO I 4  Integer  variables  are  set  to  2  bytes. 

3.  /N0F77  Where  the  compiler  needs  to  choose  - 
between  Fortran  77  and  Fortran  4+*  the  latter 

is  chosen. 

4.  /LIST  A  compiler  listing  is  made  and 
stored  on  the  disk  as  TLQSS22CC.  LIS. 


b.  Linking  (  LINK  command  ) 

Linking  is  straight-f orward.  No  extra  routines 
are  linked  with  TL0SS22CC. 


2.  Running 

TL0SS22CC  is  normally  run  through  a  batch  stream* 
in  particular  the  SLOW  batch  stream*  since  this  is  the 
cheapest  route*  especially  at  night.  However*  there  is 
another  batch  stream  with  higher  priority  (  and  higher 
cost  )  that  can  be  used  for  quicker  results*  the 
SYS*BATCH  batch  stream.  In  order  to  use  either  of  these 
two  batch  streams: 

a.  The  executable  file  (  TLQSS22CC.EXE  >  must  be  ready 
to  use.  The  LINK  operation  satisfies  this  requirement. 

b.  The  data  file  must  also  be  ready  to  use.  This  is 
normally  prepared  by  the  user  using  the  editor* 

(  and  the  TL0SS22CC  USER'S  MANUAL,  if  necessary  ) 
at  a  terminal. 

c.  Submit  the  command  file*  TLOSS22CC.COM*  to  the 
batch  stream  (  or  batch  queue  >  using  a  command 
such  as  one  of  these: 

SUBM I T/NAME-LTT2SMAR /QUEUE-SLOW/ AFTER-20: 00  TL0SS22CC 
SLOW  queue*  Job  starts  at  8: 00  P.  M. 

SUBM I T/N AME-LTT2 SM AR /QUEUE-SLOW  TL0SS22C C 
SLOW  queue*  Job  starts  now 

SUBMIT/NAME-LTT25MAR  TL0SS22CC 

SYS4BATCH  queue  (  the  'regular'  queue  )* 

Job  starts  now 


d.  The  following  command  file*  also  prepared  using 
the  editor*  is  a  typical  file  submitted*  using  one 
of  the  previous  submit  commands. 

* ASSIGN  LTT29MAR  SYS4PRINT  (1) 

4 ASSIGN  LTT25MAR  SYS40UTPUT  (2> 

•ASSIGN  DRO: LTT25MAR.  PLT  F0R007  (3) 

•ASSIGN  DRO:  LTT25MAR.  OUT  F0R006  <4> 

•ASSIGN  DRO: SHALLOW.  DAE  F0R005  (5) 

•TYPE  TL0SS22CC  (6) 

•COPY  DRO: SHALLOW.  DAE  DRO: LTT29MAR.  DAT  (7) 
•RUN  TL0SS22CC  (8) 

This  run  is  called  the  LTT25MAR  run.  'LTT'  stands  for 
LOSS  TWENTY  TWO  C.  '25MAR '  indicates  that  this  is  the 
twenty-fifth  run  in  March. 


The  first  two  lints  in  this  command  filt  put  tht 
log  filt  on  tht  disk  instead  of  tht  printtr. 

Tht  third  lint  puts  unit  7  output  into  filt  LTT29MAR.  PLT. 
Tht  fourth  lint  puts  unit  6  output  into  filt  LTT29MAR.  OUT. 
Tht  fifth  lint  draws  unit  9  input  from  SHALLOW.  DAE. 

Tht  sixth  lint  copits  TL0SS22CC.  LIS  to  LTT25MAR.  LIS. 

Tht  stvtnth  lint  copits  SHALLOW.  DAE  to  LTT23MAR.DAT. 

Tht  tighth  lint  runs  program  TL0SS22CC. 

A  command  filt  such  as  this  yields  a  filt  group; 
tach  filt  having  tht  samt  name  but  a  difftrtnt 
extension.  These  files  art  then  treated 
collectively.  They  art  printed  together  and  stored 
on  tape  together.  Such  a  group  is#  for  example: 

LTT25HAR.DAT  a  copy  of  the  data  file 

LTT29MAR.  LIS  the  program  listing  file 

LTT23MAR.LOG  the  log  file 

LTT23MAR.  OUT  the  output  file 

LTT29MAR.  PLT  the  numerical  output  file 

The  user  normally  gets  a  complete  printout  of  the 
run  by  typing  the  command: 

PRINT  LTT29MAR.  *;* 

on  the  terminal. 

Copying  Runs  to  Tape  (  Archiving  ) 

Periodically!  it  is  necessary  to  clear  the  old  runs 
off  tht  disk  and  put  them  on  tape.  This  is  a  very  simple 
operation  and  proceeds  as  follows: 

a.  Initialize  a  tape.  Give  it  a  name  that  is  somehow 
indicative  of  the  present  month.  For  example: 

INITIALIZE/DENSITY-1600  MT1:  MAR APR 

This  tape  would  have  runs  made  in  March  and  April. 

This  step  is  not  necessary  if  the  user  is  adding  runs 
to  a  tape  that  is  already  partly  filled.  In  such  a 
case#  be  cartful  not  to  initialize! 

b.  Mount  the  tape  with  the  MOUNT  command. 

MOUNT  MT1 :  MAR APR 

This  is  not  the  same  as  physically  mounting  the  tape 
on  the  tape  drive#  which  is  assumed  to  have  taken 
place  previously. 

c.  Copy  all  the  run  files  from  disk  to  tape  with  the 
COPY  command.  For  example  in  this  case: 

COPY  *MAR.  *;*  MT1 : 

COPY  *APR.  *;*  MT1: 

would  copy  all  the  March  and  all  the  April  runs  onto 
the  tape  mounted  on  tape  drive  1. 

d.  For  future  reference#  it  is  good  to  get  a  list  of 
names  of  the  files  on  the  tape.  This  is  accomplished 
as  follows: 

DIRECTORY/PRINT  MT1 : 

The  printout  resulting  from  this  command  does  not 
have  the  tape  name  on  it#  so  the  user  should  write 
this  on  immediately. 

If#  for  some  reason#  this  printout  is  lost#  the  tape 
can  be  remounted  and  the  directory  can  be  taken 
again.  But  this  is  no  reason  to  be  careless  with 
the  original  directory  printout. 
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e.  Dismount  ths  tap*  with  th*  dismount  command: 

DISMOUNT  MT1  : 

f.  Deist*  th*  files  on  th*  oisk  that  have  just  b**n 
copied  to  th*  tap*.  In  this  ease: 

DELETE  *MAR.  *;  * 

DELETE  *APR.  *;  * 

4.  Distributing  th*  Program  by  Tap* 

Th*  program  and  test  eases  ar*  distributed  using 
tapes  with  tn*  sam*  format  as  th*  archival  tapes.  This 
is  th*  Files-11  format,  peculiar  to  Digital  Equipment 
Corporation.  Th*  tapes  from  this  system  have  a  density 
of  1600  bpi.  and  a  blocksiz*  of  2048  bytes. 


If  any  further  questions  arise,  feel  free  to  call: 

Dr.  H.  Kutschal*  (914)  399-2900  *382 

Mr.  L.  Carroll  (914)  399-2900  *221 
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TEST  CASE  DATA 


05  1 

0 

DATA  FOR 

TL0SS22CC 

00.  0 

1500.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

50. 

1500.  0 

0.  000 

0.  00001 

0. 00000 

1.  00 

75. 

1500.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

4200. 

1300.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

4200. 

1600.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

2 

3 

2  1024 

-1 

0 

001  020.  0  10000.  0 

1499.  00 

1600.  01 

0.  0 

00.  0 

0.  00 

00.  00 

0.  000 

0 

1 

1 

03  1 

0 

DATA  FOR 

TL0SS22CC 

00.  0 

1300.  0 

0.  000 

0. 00001 

0.  00000 

1.  00 

30. 

1300.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

73. 

1300.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

4200. 

1300.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

4200. 

1600.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

2 

3 

2  1024 

-1 

0 

001  020.  0  10000. 0 

1499.  00 

1600.  01 

0.  0 

00.  0 

0.  00 

00.  00 

0.  000 

0 

0 

2 

03  1 

0 

DATA  FOR 

TL0SS22CC 

00.  0 

1300.  0 

0.  000 

0. 00001 

0.  00000 

1.  00 

SO. 

1300.  0 

0.  000 

0.  00001 

0.  00000 

l.  00 

75. 

1500.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

4200. 

1500.  0 

0.  000 

0. 00001 

0. 00000 

1.  00 

4200. 

1600.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

2 

3 

2  1024 

-1 

0 

001  020.  0  10000. 0 

1499.  00 

1600.  01 

0.  0 

00.  0 

0.  00 

00.  00 

0.  000 

0 

1 

3 

03  1 

0 

DATA  FOR 

TLQSS22CC 

00.  0 

1300  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

50. 

1500.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

75. 

1500.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

4200. 

1500.  0 

0.  000 

0.  00001 

0.  00000 

1.  00 

4200. 

1600. 0 

0.  000 

0.  00001 

0.  00000 

1.  00 

2 

3 

2  1024 

00 

0 

001  020.  0  10000.  0 

1499.  00 

1600.  01 

0.  0 

00.  0 

0.  00 

00.  00 

0.  000 

0 

0 

4 
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TEST  CASE  DATA 


04  1 

0 

00.  0 

1500.  0 

0.  OCO 

00.  0 

1500.  0 

0.  000 

4125. 

1500.  0 

0.  000 

4125. 

1600.  0 

0.  000 

2 

2  2  1024 

00 

0.  0 

1500.  0 

0.  OOOl 

DATA  FOR  TLQSS22CC 

O.OOOOl  0.00000  1.00 

0.00001  0.00000  1.00 

0.00001  0.00000  1.00 

0.00001  0.00000  1.00 

0  001  020.  0  10000.  0 
01.00  25.00  0 


1499.  00 
1  0 


1600.  01 


: j 


\1 
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RANGE  (KM) 
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06  t 
0.  00 
50.  00 
100.  0 
100.  0 
125.  0 
125.  0 
2 

O.  0 


TEST  CASE  DATA 


1500.  0 
1500.  0 
1500.  0 
3000.  0 
3000. 0 
3000.  0 
02  1024 
1500.  0 


1400.  0 
1400.  0 
1400.  0 
00  000 
. 0000000 


6-APR-1983  15:31 


0000000 
0000000 
0000000 
00006286 
00006286 
00006286 
1  020.  0 


DATA  FOR  TL0SS22CC 
1.  00 


1.  00 


1.  00 
1.  00 

0001281  2.  50 
0001281  2.  50 
0001281  2.  50 

i  10000. 


0000.  0 


1000.  0 
0 


8000.  0 
0 
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TEST  CASE  DATA 


08  1 

X 

000.  0 

1536.  5 

.  0000100 

1.  0 

152.  4 

1539.  2 

.  0000100 

1.  0 

253.  9 

1530.  0 

.  0000100 

1.  0 

M 

406.  3 

1501.  1 

.  0000100 

1.  0 

■ 

863.  5 

1475.  0 

. 0000100 

1.  0 

1015.  9 

1471  9 

. 0000100 

1.  0 

5587. 9 

1549.  6 

. 0000100 

1.  0 

5587. 9 

1549.  6 

.  0000100 

1.  9 

\', 

3 

5  2  512 

01  00 

1  25. 

10000.  0 

%  t 

00.  0 

0000.  0 

0.  00000 

0.  0 

000.  0 

3 

5  2  512 

00  00 

1  25. 

10000.  0 

00.  0 

1536.  5 

0.  00001 

1.  0 

000.  0 

1470.  0 
0 

1470.  0 
0 


1549.  4 
1549.  4 
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Introduction 

TLOSSIOHH*  our  transmission  loss  (  hereafter  abbreviated  as  TL  ) 
program,  was  completed  in  June  of  1982. 

TL0SS10HH  has  these  features  and  options  for  the  user: 

1.  A  choice  of  computational  methods; 

a.  The  FFP  method 

b.  The  Normal  Mode  method 

2.  Optional  horizontal  or  vertical  displacement  in  addition 
to  the  usual  pressure  values  for  TL. 

3.  A  comprehensive  selection  of  modelling  parameters. 

These  include: 

a.  A  detailed  sound  velocity  profile 

b.  Choice  of  source  depth 

c.  Choice  of  detector  depth 

d.  Choice  of  frequency 

e.  Choice  of  wave  number  integrand  limits 

f.  Optional  beaming*  as  it  is  affected  by 
integrand  tapering  (  in  the  FFP  mode  ) 

g.  Optional  surface  roughness 

h.  Choice  of  either  fluid  or  vacuum  upper  half 
space 

i.  Possibility  of  the  source  and  detector  in 
the  fluid  upper  half  space 

j.  Optional  internal  roughness 

Ic.  Choice  of  initial  range*  range  step,  and 
range  count  (  for  the  Normal  Mode  method  ) 

4.  Output.  This  includes: 

a.  A  printout  of  all  input 

b.  The  TL  plot 

c.  The  wave  integrand  plot  of  the  FFP;  either 
short  or  long. 

d.  A  table  of  the  denominator  values  for  the 
Normal  Mode  method 

e.  A  table  of  roots  and  their  corresponding 
phase  velocities*  slopes*  and  excitation 
functions  for  the  Normal  Mode  method 

Note:  All  short  plots  except  the  wave  integrands 
are  repeated*  the  repeated  plot  being  a  smoothed 
copy  of  the  first  and  having  the  same  bounds  as 
the  first. 

5.  Additional  Optional  Tables  and  Plots.  These  include: 

a.  The  TL  plot  in  nautical  miles*  if  desired  and 
if  the  range  is  sufficiently  long.  The  original 
TL  plot  is  in  kilometers. 

b.  The  denominator  plot  for  the  Normal  Mode  method. 

6.  Numerical  output  of  the  TL  for  further  processing  or  plotting 
Formerly*  this  was  called  punched  output,  since  it  went  out 

to  the  card  punch. 

7.  A  profile  check  for  multiple  ducts. 

This  check  is  made  for  the  Normal  Mode  method. 

If  multiple  ducts  are  present*  special  care  is  taken 
to  locate  those  'hidden'  roots  which  are  only  found 
by  reducing  the  step  size  of  the  denominator  function. 
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TLOSS 1 OHH  USER'S  MANUAL 
Input  Variables  and  Formats 

Program  input  for  TLGSSIOHH  includes  both  the  physical  model 
and  program  control.  This  section  of  the  manual  will  list 
and  explain  each  line  of  the  input. 

1.  Entry  1>  one  line 

MMAXF*  NCARD  variable  names 

(  213  )  format 


MMAXF  is  a  control  variable  which  is  used  to  read 
in  the  correct  number  of  points  in  the  sound  velocity 
profile.  (  see  2.  -  next  Entry  > 

NCARD  is  a  control  variable  that*  when  set  non-zero 
allows  the  program  to  put  out  TL  data  on  unit  7  in 
a  form  which  is  readable  as  input  bu  another  program. 

A  zero  value  will  prevent  any  output  on  unit  7.  This 
is  normally  used  to  make  graphic  plots.  (  Originally* 
on  the  IBM  360/91*  unit  7  was  the  punched  card  unit* 
hence  the  name  NCARD.  ) 

2.  Entry  2*  MMAXF  lines  (  max.  of  49  lines  ) 

((  V(M*N>.  N-l*  7  >.  M-l.  MMAXF  >  variable  name 

(  7D10.  0  )  format 

V  is  an  array  that  contains  the  sound  velocity  profile. 
For  each  point  in  the  profile*  and  on  each  line* 

(  one  point  per  line  )*  the  program  expects: 

a.  a  depth  value*  in  meters 

b.  a  sound  <  pressure  )  velocity  value*  in  meters  per  second 

c.  a  sound  (  shear  >  velocity  value*  in  meters  per  second* 

if  applicable 

d.  a  corresponding  attenuation  value  for  pressure 

e.  a  corresponding  attenuation  value  for  shear, 

if  applicable 

f.  a  density  value*  in  grams  per  c. c. .  or  specific  density 

g.  a  roughness  value,  if  applicable 

Note:  No  depth  value  can  be  less  than  a  previous 

depth  value.  The  program  will  stop  if  this  is  so. 

3.  Entry  3.  one  line 

KSF.  LRF,  MM.  LS.  KKK.  IBEAM*  IMAX,  FREQ.  PO.  C.  CMAX 


<  713.  4D10.  0  ) 

KSF  is  a  pointer  to  the  source  depth  value. 

It  points  to  one  of  the  points  in  the  sound  velocity 
profile*  the  one  which  is  the  source  depth  value. 

LRF  is  a  pointer  to  the  detector  (  or  receiver  > 


variable  name! 
format 


depth  valu 


s  a  poin 
e.  It's 


function  is  similar  to  KSF. 


MM  is  a  TL  plot  pruning  factor,  with  a  minimum 
acceptable  value  of  2.  MM  is  applicable  to  FFP  only. 

If  an  MM  value  of  N  is  read  in*  then  the  first  1/N 
of  the  plot  will  be  plotted.  This  is  useful  either  to 
get  a  more  detailed  look  at  the  output  or  to  prune 
away  that  part  of  the  plot  affected  by  range  aliasing. 

LS  is  the  number  of  points  used  in  the  wave 
integration.  It  also  affects  the  range  of  the  TL  plot, 
for  the  FFP  method. 

KKK  is  a  control  variable  that  indicates: 

a.  If  zero*  this  is  the  last  case. 

b.  If  negative*  the  next  case  starts  with  Entry  1* 

(  new  profile  ). 

c.  If  positive*  the  next  case  starts  with  Entry  3. 

IBEAM  is  a  control  variable  that*  when  set  non¬ 
zero*  requires  the  program  to  taper  the  wave  integrand 
at  each  end. 


IMAX  is  the  rang*  count  or  th*  number  of  point*  in 
th*  TL  plot  for  a  Normal  Mod*  case. 

FREQ  is  th*  frequency  for  this  case*  in  Hertz. 

PO  is  th*  source  level.  For  TL  PO  *  1.0 

C  is  th*  lower  bound  sound  velocity  in  th*  wav* 
integrationi  in  meters  per  second.  For  the  Normal  Mod*  m*thod< 
if  a  zero  is  entered/  a  default  value  will  be  found 
from  th*  sound  velocity  profile. 

CMAX  is  the  upper  bound  sound  velocity  in  th*  wav* 
integration#  in  meters  per  second.  If  a  zero  is  entered 
her*/  the  wave  number  integration  will  start  from  zero 
for  th*  FFP  method/  or,  for  th*  Normal  Mod*  method/  a 
default  value  will  be  found  from  th*  sound  velocity 
profile. 

Entry  4/  on*  line 

REFLTOz  ALPHAO.  ATTENO/  RHOO,  H/  1NAUT,  LONG.  KHPV,  IDPL 

variable  names 

(  3D10.  0/  415  )  format 

REFLTO  is  the  RMS  surface  roughness  value/  in  meters. 

If  this  value  is  nan-zero/  th*  program  explicitly  sets 
th*  next  three  values  to  zero. 

ALPHAO  is/  if  present/  the  sound  velocity  in  th* 
fluid  upper  half  space/  in  meters  per  second.  This  is 
th*  variable  that  determines  whether  th*  present  case  has 
either  a  fluid  or  vacuum  upper  half  space. 

ATTENO  is/  if  present/  the  attenuation  in  th* 
fluid  upper  half  space. 

RHOO  is/  if  present/  the  density  of  the  fluid 
upper  half  space/  in  grams  per  c.  c. 

H  is#  if  present/  th*  distance  between  the  source 
and  detector  in  th*  fluid  upper  half  space.  To  us*  this 
option: 

a.  Th*  first  two  sound  velocity  depths  must  be  zero. 

b.  There  must  be  some  variation  in  the  pressure 
sound  velocity  profile. 

c.  There  must  not  be  any  shear  sound  velocity. 

d.  ALPHAO  and  H  must  be  present. 

INAUT  is  a  control  variable  that/  when  set  non¬ 
zero/  allows  th*  plotting  of  th*  TL  plot  in  nautical 
miles  if  th*  rana*  is  greater  than  100  kilometers. 

This  is  for  th*  FFP  method  only. 

LONQ  is  a  control  variable  that/  when  set  non¬ 
zero/  switches  the  wave  integrand  plot  from  short  to 
long.  This  is  used  when  it  is  necessary  to  examine 
each  point  of  th*  wave  integrand.  FFP  only. 

KHPV  is  the  control  variable  that  determines 
whether  output  will  be: 

a.  pressure/  if  zero. 

b.  vertical  displacement/  if  positive. 

c.  horizontal  displacement/  if  negative. 

IDPL  is  th*  control  variable  that/  when  set/  allows  th* 
denominator  function  to  be  plotted.  Normal  Mod*  method  only. 


Entry  5.  ona  line 
RA,  OLRA 
<  2D10.  0  ) 


variabla  nama 
format 


RA  is  tha  initial  rang*  of  tha  TL  plot#  in  matars. 

OLRA  is  tha  ranga  stap  of  tha  TL  plot#  in  matars. 

Thasa  two  variablas  ara  raad  in  for  tha  Normal  Moda 
mathod.  Sinca  tha  TL  plot  is  in  kilomatars#  it  is  vary 
important  to  ba  cartful  with  tha  units  of  thasa  two  variablas. 

;a: 

Tha  program  datarminas  tha  computational  mathod# 

(  whathar  FFP  or  Normal  Moda  )>  by  axaminina  all  of 
tha  attanuation  valuas  and  tha  surfaca  rougnnass  valua. 

If  thay  ara  all  taro#  tha  Normal  Moda  mathod  is  salactad. 
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Operating  Procedures 


The  primary  computer  at  Lamont-Doherty  Geological  Observatory 
is  a  VAX  11/780.  which  runs  under  the  VMS  operating  system. 

Listed  below  are  four  distinctly  different  operations  that 
are  used  to  process  TL0SS10HH. 


1. 


Compiling-Linking 
a.  Compiling 
Comp i 1 ing 
has  always 
on: 


(  FORT  command  ) 

the  TL0SS10HH  program  on  the  VAX  11/780 
been  done  with  tne  following  switches 


1.  /CHECK  Forbids  array  pointer  overflow. 

2.  /NQI4  Integer  variables  are  set  to  2  bytes. 

3.  /N0F77  Where  the  compiler  needs  to  choose 
between  Fortran  77  and  Fortran  4+.  the  latter 
is  chosen. 

4.  /LIST  A  compiler  listing  is  made  and 
stored  on  the  disk  as  TLOSSIOHH.  LIS. 


b.  Linking  (  LINK  command  ) 

Linking  is  straight-f orward.  No  extra  routines 
are  linked  with  TLOSSIOHH. 


2.  Running 

TLOSSIOHH  is  normally  run  through  a  batch  stream, 
in  particular  the  SLOW  batch  stream,  since  this  is  the 
cheapest  route,  especially  at  night.  However,  there  is 
another  batch  stream  with  higher  priority  (  and  higher 
cost  )  that  can  be  used  for  quicker  results,  the 
SYS4BATCH  batch  stream.  In  order  to  use  either  of  these 
two  batch  streams: 


a.  The  executable  file  (  TLQSS10HH.EXE  )  must  be  ready 
to  use.  The  LINK  operation  satisfies  this  requirement. 


b.  The  data  file  must  also  be  ready  to  use.  This 
normally  prepared  by  the  user  using  the  editor. 

(  and  the  TLOSSIOHH  USER'S  MANUAL,  if  necessary  ) 
at  a  terminal. 


is 


c.  Submit  the  command  file.  TL0SS10HH.COM.  to  the 
batch  stream  (  or  batch  queue  >  using  a  command 
such  as  one  of  these: 

SUBMIT/NAME-LTH22JUN/QUEUE-SL0W/AFTER-20: 00  TLOSSIOHH 
SLOW  queue.  Job  starts  at  8: 00  P.  M. 

SUBM I T /NAME-LTH22 JUN/ QUEUE-SLOW  TLOSS 1 OHH 
SLOW  queue.  Job  starts  now 

SUBMIT/NAME-LTH22JUN  TLOSSIOHH 

SYStBATCH  queue  (  the  'regular'  queue  >. 

Job  starts  now 


d.  The  following  command  file,  also  prepared  using 
the  editor,  is  a  typical  file  submitted,  using  one 
of  the  previous  submit  commands. 

•ASSIGN  LTH22JUN  SYStPRINT  (1) 

•ASSIGN  LTH22JUN  SYS»OUTPUT  (2) 

•ASSIGN  DRO: LTH22JUN.  PLT  F0R007  <3> 

•ASS I ON  DRO: LTH22JUN.  OUT  F0R006  (4) 

•ASSIGN  DRO: STANDARD.  DAF  F0R009  (9) 

•TYPE  TLOSSIOHH  (6) 

•COPY  DRO: STANDARD.  DAF  DRO: LTH22JUN.  DAT  <7> 
•RUN  TLOSSIOHH  (8) 

This  run  is  called  the  LTH22JUN  run.  'LTH'  stands  for 
LOSS  TEN  H.  '22 JUN '  indicates  that  this  is  the 

22nd  run  in  June. 


• 


'.  .V 


The  first  two  linos  in  this  command  file  put  the 
log  file  on  the  disk  instead  of  th#  printer. 

Th#  third  lin#  puts  unit  7  output  into  file  LTH22JUN.  PLT. 
Th#  fourth  lin#  puts  unit  6  output  into  fil#  LTH22JUN.  OUT. 
Th#  fifth  lin#  draws  unit  9  input  from  STANDARD.  OAF. 

Th#  sixth  lin#  copies  TLOSSIOHH.  LIS  to  LTH22JUN.  LIS. 

Th#  seventh  lin#  copies  STANDARD.  DAF  to  LTH22JUN.  DAT. 

Th#  eighth  lin#  runs  program  TLOSSIOHH. 

A  command  fil#  such  as  this  yields  a  fil#  group; 
each  fil#  having  the  same  name  but  a  different 
extension.  Th#s#  files  are  then  treated 
collectively.  They  are  printed  together  and  stored 
on  tape  together.  Such  a  group  is*  for  example: 


LTH22JUN.  DAT 
LTH22JUN.  LIS 
LTH22JUN.  LOG 
LTH22JUN.  OUT 
LTH22JUN.  PLT 


a  copy  of  the  data  file 
the  program  listing  file 
the  log  file 
the  output  file 
the  numerical  output  file 


The  user  normally  gets  a  complete  printout  of  the 
run  by  typing  the  command: 

PRINT  LTH22JUN.  *; * 

on  the  terminal. 

Copying  Runs  to  Tape  (  Archiving  ) 

Periodically*  it  is  necessary  to  clear  the  old  runs 
off  the  disk  and  put  them  on  tape.  This  is  a  very  simple 
operation  and  proceeds  as  follows: 

a.  Initialize  a  tape.  Give  it  a  name  that  is  somehow 
indicative  of  the  present  month.  For  example: 

INITIALIZE/DENS ITY*1600  MT1:  JUNJUL 

This  tape  would  have  runs  made  in  June  and  July. 

This  step  is  not  necessary  if  the  user  is  adding  runs 
to  a  tape  that  is  already  partly  filled.  In  such  a 
case*  be  careful  not  to  initialize! 

b.  Mount  the  tape  with  the  MOUNT  command. 

MOUNT  MT1 :  JUNJUL 

This  is  not  the  same  as  physically  mounting  the  tape 
an  the  tape  drive*  which  is  assumed  to  have  taken 
place  previously. 

c.  Copy  all  the  run  files  from  disk  to  tape  with  the 
COPY  command.  For  example  in  this  cast: 

COPY  *JUN.  *;*  MTl: 

COPY  *JUL.  *;  *  MTl: 

would  copy  all  the  June  and  all  the  July  runs  onto 
the  tape  mounted  on  tape  drive  1. 

d.  For  future  reference*  it  is  good  to  get  a  list  of 
names  of  the  files  on  the  tape.  This  is  accomplished 
as  follows: 

DIRECTORY/PRINT  MTl: 

The  printout  resulting  from  this  command  does  not 
have  the  tape  name  on  it*  so  the  user  should  write 
this  on  immediately. 

If*  for  some  reason*  this  printout  is  lost*  the  tape 
can  be  remounted  and  the  directory  can  be  taken 
again.  But  this  is  no  reason  to  be  careless  with 
tne  original  directory  printout. 
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e.  Dismount  the  taps  with  the  dismount  command: 

DISMOUNT  MTi: 

f.  Dalata  th*  files  on  the  disk  that  have  just  bean 
copied  to  the  tape.  In  this  case: 

DELETE  *JUN.  *«  * 

DELETE  *JUL.  *i  * 

Distributing  the  Program  by  Tape 

The  program  and  cast  cases  are  distributed  using 
tapes  with  the  same  format  as  the  archival  tapes.  This 
is  the  Files-11  format*  peculiar  to  Digital  Equipment 
Corporation.  The  tapes  from  this  sustem  have  a  density 
of  1600  bpi»  and  a  blocksize  of  2048  bytes. 


bytes. 


If  any  further  questions  arise*  feel  free  to  call: 

Dr.  H.  Kutschale  <914>  359-2900  x382 


Mr.  L.  Carroll 


(914)  359-2900  x221 


wv 


LTH010CT.  DAT;  2 


6-APR-1983  15:  31 


TEST  CASE  DATA 


0.  00 

1500.  0 

. 0000150 

50.  00 

1500.  0 

. 0000150 

100.  0 

1500.  0 

. 0000150 

100.  0 

3000.  0 

1400.  0 

. 00007700 

500.  0 

3000.  0 

1400.  0 

. 00007700 

500.  0 

3000.  0 

1400.  0 

. 00007700 

2 

3  02  1024 

-1  000 

1  020 

0.  0 

00. 

1500.  0 

CO.  0 

. 0000150 

1.  00  ( 

0.  00 

1500. 0 

. 0000150 

50.  00 

1500.  0 

. 0000150 

100.  0 

1500.  0 

. 0000150 

100.  0 

3000.  0 

1400.  0 

. 00007700 

500.  0 

3000.  0 

1400.  0 

. 00007700 

500.  0 

3000.  0 

1400.  0 

.  00007700 

2 

3  02  1024 

-1  000 

1  020 

0.  0 

00. 

1500.  0 

00.  0 

. 0000150 

1.  00  ( 

0.  00 

1500.  0 

. 0000150 

50.  00 

1500.  0 

. 0000150 

100.  0 

1500.  0 

. 0000150 

100.  0 

3000.  0 

1400.  0 

. 00007700 

500.  0 

3000.  0 

1400.  0 

. 00007700 

500.  0 

3000. 0 

1400.  0 

.  00007700 

2 

3  02  1024 

00  000 

1  020 

0.  0 

00. 

1500.  0 

00.  0 

. 0000150 

1.  00  i 

DATA  FOR  TLOSSIOHH 
1.  00 
1.  00 
1.  00 

. 0001600  2.  50 
.  0001600  2.  50 
0001600  2.  50 

>.  o  10000.  1000.  0 

0000. 0  01 

DATA  FOR  TLOSSIOHH 
1.  00 
1.  00 
1.  00 

.  0001600  2.  50 
.  0001600  2.  50 
. 0001600  2.  SO 
0  10000.  1000.  0 

0000. 0  0  1 

DATA  FOR  TLOSSIOHH 
1.  00 
1.  00 
1.  00 

.  0001600  2.  50 
.  0001600  2.  50 
.  0001600  2.  50 
.  0  10000.  1000.  0 

0000.  0  0  1 


8000.  0 
1 


8000.  0 
1 


8000.  0 
1 
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A.  Introduction 


TL0SS17HH.  our  nor*  complete  transmission  loss  (  hereafter 
abbreviated  as  TL  >  program*  was  completed  in  Juna  of  1982. 

It  is  an  extension  of  TLOSSIOHH*  tha  intarim  standard  TL 
program. 

TL0SS17HH  has  thasa  faaturas  and  options  for  tha  user: 

1.  A  choica  of  computational  methods* 

a.  Tha  FFP  mathod 

b.  Tha  Normal  Moda  mathod 

If  tha  Normal  Moda  mathod  is  chosen*  tha  usar  also 
has  tha  option  of  including  aithar  of  two  branch  cut 
intagrals. 

2.  Optional  horizontal  or  vortical  d isp lacamant  in  addition 
to  tha  usual  prassura  valuas  for  TL. 

3.  A  comprahansiva  solaction  of  modalling  paramatars. 

Thasa  includa: 

a.  A  datailad  sound  velocity  profila 

b.  Choica  of  sourco  dapth 

c.  Choica  of  datoctor  dapth 

d.  Choica  of  frequency 

a.  Choica  of  wava  numbar  intagrand  limits 

f.  Optional  beaming*  as  it  is  affactad  by 
intagrand  taparing  (  in  tha  FFP  moda  > 

g.  Optional  surfaca  roughnass 

n.  Choica  of  aithar  fluid  or  vacuum  uppar  half 
space 

i.  Possibility  of  tha  sourca  and  datactor  in 
tha  fluid  uppar  half  spaca 

j.  Optional  internal  roughnass 


Optional  intarnal  roughnass 
Choica  of  initial  range*  ranga  stap>  and 
ranga  count  (  for  tha  Normal 


oda  mathod  ) 


4.  Output.  This  includas: 

a.  A  printout  of  all  input 

b.  Tha  TL  plot 

c.  Tha  wava  intagrand  plot  of  tha  FFP;  aithar 
short  or  long. 

d.  A  tabla  of  tha  danominator  valuas  for  tha 
Normal  Moda  mathod 

a.  A  tabla  of  roots  and  thair  corraspond ing 
phasa  velocities*  slopes.  and  excitation 
functions  for  tha  Normal  Moda  mathod 

Nota:  All  short  plots  axcapt  tha  wava  intagrands 

ara  repeated*  tha  rapaatad  plot  baing  a  smoothad 

copy  of  tha  first  and  having  tha  sama  bounds  as 

tha  first. 

5.  Additional  Optional  Tablas  and  Plots.  Thasa  includa: 

a.  Tha  TL  plot  in  nautical  miles*  if  dasirad  and 
if  tha  ranga  is  sufficiantly  long.  Tha  original 
TL  plot  is  in  kilomatars. 

b.  Tha  danominator  plot  for  tha  Normal  Moda  mathod. 

c.  Modal  TL  plots*  Branch  TL  plots*  and  Total  TL 
plots*  if  a  branch  cut  option  is  chosan. 

6.  Numarical  output  of  tha  TL  for  furthar  procassing  or  plotting. 
Formarly*  this  was  callad  punchad  output*  since  it  want  out 

to  tha  card  punch. 

7.  A  profila  check  for  multiple  ducts. 

This  check  is  made  for  tha  Normal  Moda  mathod. 

If  multiple  ducts  ara  present*  special  care  is  taken 
to  locate  those  'hidden'  roots  which  ara  only  found 
by  reducing  the  step  size  of  tha  danominator  function. 
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Input  Variable*  and  Formats 

Program  input  for  TL0SS17HH  includes  both  tha  physical  modal 
and  program  control.  This  section  of  tha  manual  will  list 
and  explain  each  line  of  tha  input. 


1.  Entry  1*  one  line 
MMAXF,  NCARD 
(  213  > 


variable  names 
format 


MMAXF  is  a  control  variable  which  is  used  to  read 
in  the  correct  number  of  points  in  the  sound  velocity 
profile.  (  see  2.  -  next  Entry  > 

NCARD  is  a  control  variable  that*  when  set  non-zero 
allows  the  program  to  put  out  TL  data  on  unit  7  in 
a  form  which  is  readable  as  input  by  another  program. 

A  zero  value  will  prevent  any  output  on  unit  7.  This 
is  normally  used  to  make  graphic  plots.  (  Originally* 
on  the  IBM  360/91*  unit  7  was  the  punched  card  unit* 
hence  the  name  NCARD.  ) 

2.  Entry  2<  MMAXF  lines  (  max.  of  49  lines  ) 

((  V(M,  N) *  N-1.7  ).  M-l. MMAXF  )  variable  name 

(  7D10.  0  >  format 

V  is  an  array  that  contains  the  sound  velocity  profile. 
For  each  point  in  the  profile*  and  on  each  line* 

(  one  point  per  line  >*  the  program  expects: 

a.  a  depth  value*  in  meters 

b.  a  sound  (  pressure  )  velocity  value*  in  meters  per  second 

c.  a  sound  (  shear  )  velocity  value*  in  meters  per  second, 

if  applicable 

d.  a  corresponding  attenuation  value  for  pressure 

e.  a  corresponding  attenuation  value  for  shear* 

if  applicable 

f.  a  density  value*  in  grams  per  c. c. *  or  specific  density 

g.  a  roughness  value*  if  applicable 

Note:  No  depth  value  can  be  less  then  a  previous 

depth  value.  The  program  will  stop  if  this  is  so. 

3.  Entry  3*  one  line 

KSF.  CRF,  MM.  LS.  KKK.  IBEAM*  IMAX.  FREQ*  PO.  C,  CMAX 


(  715,  4D10.  0  ) 


variable  names 
format 


KSF  is  a  pointer  to  the  source  depth  value. 

It  points  to  one  of  the  points  in  the  sound  velocity 
profile**  the  one  which  is  the  source  depth  value. 

LRF  is  a  pointer  to  the  detector  (  or  receiver  > 
depth  value.  It's  function  is  similar  to  KSF. 

MM  is  a  TL  plot  pruning  factor,  with  a  minimum 
acceptable  value  of  2.  MM  Is  applicable  to  FFP  only. 

If  an  MM  value  of  N  is  read  in*  then  the  first  1/N 
of  the  plot  will  be  plotted.  This  is  useful  either  to 
get  a  more  detailed  look  at  the  output  or  to  prune 
away  that  part  of  the  plot  affected  by  range  aliasing. 

LS  is  the  number  of  points  used  in  the  wave 
integration.  It  also  affects  the  range  of  the  TL  plot* 
for  the  FFP  method. 

KKK  is  a  control  variable  that  indicates: 

a.  If  zero,  this  is  the  last  case. 

b.  If  negative*  the  next  case  starts  with  Entry  1* 

<  new  profile  >. 

c.  If  positive,  the  next  case  starts  with  Entry  3. 

IBEAM  is  a  control  variable  that*  when  set  non¬ 
zero,  requires  the  program  to  taper  the  wave  integrand 
at  each  end. 


IW 


IMAX  is  the  rang*  count  or  the  number  of  points  in 
the  TL  plot  for  a  Normal  Node  case. 

FREQ  is  the  frequency  for  this  case/  in  Hertz. 

PO  is  the  source  level.  For  TL  PO  ■  1.0 

C  is  the  lower  bound  sound  velocity  in  the  wave 
integration/  in  meters  per  second.  For  the  Normal  Mode  method 
if  a  zero  is  entered/  a  default  value  will  be  found 
from  the  sound  velocity  profile. 

CMAX  is  the  upper  bound  sound  velocity  in  the  wave 
integration/  in  meters  per  second.  If  a  zero  is  entered 
here/  the  wave  number  integration  will  start  from  zero 
for  the  FFP  method/  or/  for  the  Normal  Mode  method/  a 
default  value  will  be  found  from  the  sound  velocity 
profile. 

Entry  4/  one  line 

REFLTO,  ALPHAOz  ATTENO.  RHOO,  Hz  INAUT,  LONG,  KHPV,  IBR, 

IDPL  variable  names 

<  5D10.  0,  515  )  format 

REFLTO  is  the  RMS  surface  roughness  value,  in  meters. 

If  this  value  is  non-zero,  the  program  explicitly  sets 
the  next  three  values  to  zero. 

ALPHAO  is,  if  present,  the  sound  velocity  in  the 
fluid  upper  half  space,  in  meters  per  second.  This  is 
the  variable  that  determines  whether  the  present  case  has 
either  a  fluid  or  vacuum  upper  half  space. 

ATTENO  is,  if  present,  the  attenuation  in  the 
fluid  upper  half  space. 

RHOO  is,  if  present,  the  density  of  the  fluid 
upper  half  space,  in  grams  per  c.  c. 

H  is.  if  present,  the  distance  between  the  source 
and  detector  in  the  fluid  upper  half  space.  To  use  this 
option: 

a.  The  first  two  sound  velocity  depths  must  be  zero. 

b.  There  must  be  some  variation  in  the  pressure 
sound  velocity  profile. 

c.  There  must  not  be  any  shear  sound  velocity. 

d.  ALPHAO  and  H  must  be  present. 

INAUT  is  a  control  variable  that,  when  set  non¬ 
zero,  allows  the  plotting  of  the  TL  plot  in  nautical 
miles  if  the  range  is  greater  than  100  kilometers. 

This  is  for  the  FFP  method  only. 

LONG  is  a  control  variable  that,  when  set  non¬ 
zero,  switches  the  wave  integrand  plot  from  short  to 
long.  This  is  used  when  it  is  necessary  to  examine 
each  point  of  the  wave  integrand.  FFP  only. 

KHPV  is  the  control  variable  that  determines 
whether  output  will  be: 

a.  pressure,  if  zero. 

b.  vertical  displacement,  if  positive. 

c.  horizontal  displacement,  if  negative. 

IBR  is  the  control  variable  that  determines  if 
a  branch  cut  integral  will  be  included  with  the  Normal 
Mode  method. 

a.  If  zero,  no. 

b.  If  positive,  yes  <  exact  ) 

c.  If  negative,  yes  <  FFP  ) 

IDPL  is  the  control  variable  that,  when  set,  allows  the 
denominator  function  to  be  plotted.  Normal  Mode  method  only. 
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5.  Entry  5.  one  line 

RA.  DLRA  variable  names 

(  2D10.  0  )  format 

RA  is  the  initial  range  of  the  TL  plot#  in  meters. 

DLRA  is  the  range  step  of  the  TL  plot#  in  meters. 

These  two  variables  are  read  in  for  the  Normal  Mode 
method#  when  IBR  is  zero  or  positive#  only.  Since  the 
TL  plot  is  in  kilometers#  it  is  very  important  to  be 
careful  with  the  units  of  these  two  variables. 

Note: 

The  program  determines  the  computational  method# 

(  whether  FFP  or  Normal  Mode  >#  by  examining  all  of 
the  attenuation  values  and  the  surface  roughness  value. 

If  they  are  all  zero#  the  Normal  Mode  method  is  selected. 
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Operating  Procedures 

The  primary  computer  at  Lamont-Doherty  Geological  Observatory 
is  a  VAX  11/780.  which  runs  under  the  VMS  operating  system. 

Listed  below  are  four  distinctly  different  operations  that 
are  used  to  process  TL0SS17HH. 


Comp ilina -Linking 

a.  Compiling  <  FORT  c 

Compiling  the  TL0SS17HH  progr 
has  always  been  done  with  the 


command  ) 

ram  on  the  VAX  11/780 
e  following  switches 


1.  /CHECK  Forbid 

2.  /NO I 4  Intege 

3.  /N0F77  Where 
between  Fortran  77 
is  chosen. 

4.  /LIST  A  comp 
stored  on  the  disk 


Forbids  array  pointer  overflow. 
Integer  variables  are  set  to  2  bytes. 
Where  the  compiler  needs  to  choose 
ran  77  and  Fortran  4+»  the  latter 

A  compiler  listing  is  made  and 
e  disk  as  TL0SS17RH.  LIS. 


Linking  (  LINK 
Linking  is  straight-f orward. 
are  linked  with  TLOSS 17HH. 


S17HH. 


command  > 

No  extra  routines 


Running 

TLflSS 


TL0SS17HH  is  normally  run  through  a  batch  stream, 
in  particular  the  SLOW  batch  stream,  since  this  is  the 
cheapest  route,  especially  at  night.  However,  there  is 
another  batch  stream  with  higher  priority  (  and  higher 
cost  >  that  can  be  used  for  quicker  results,  the 
SYS4BATCH  batch  stream.  In  order  to  use  either  of  these 
two  batch  streams: 

a.  The  executable  file  (  TL0SS17HH.EXE  >  must  be  ready 

to  use.  The  LINK  operation  satisfies  this  requirement. 

b.  The  data  file  must  also  be  ready  to  use.  This  is 
normally  prepared  by  the  user  using  the  editor. 

(  and  the  TL0SS17HH  USER 'S  MANUAL,  if  necessary  ) 
at  a  terminal. 

c.  Submit  the  command  file.  TLOSS17HH.COM.  to  the 
batch  stream  (  or  batch  queue  >  using  a  command 
such  as  one  of  these: 

SUBMIT/NAME-LSH22JUN/QUEUE-SL0W/AFTER-20: 00  TLQSS17HH 
SLOW  queue.  Job  starts  at  8:00  P.  M. 

SUBM IT/ NAME-LSH22 JUN / QUEUE-SLOW  TLOSS 1 7HH 
SLOW  queue.  Job  starts  now 

SUBM I T /NAME-LSH22 JUN  TLOSS 17HH 

SYS4BATCH  queue  (  the  'regular'  queue  ). 

Job  starts  now 

d.  The  following  command  file,  also  prepared  using 
the  editor,  is  a  typical  file  submitted,  using  one 
of  the  previous  submit  commands. 

♦ASSIQN  LSH22JUN  SYS4PRINT  (1) 

♦ASSIGN  LSH22JUN  SYS40UTPUT  (2) 

♦ASSIGN  DRO: LSH22JUN.  PLT  F0R007  (3) 

♦ASSIGN  DRO: LSH22JUN.  OUT  F0R006  (4) 

♦ASSIGN  DRO: SIMPLE19.  DAF  F0R0O5  <3> 

♦TYPE  TLOSS 17HH  (6) 

♦COPY  DRO: SIMPLE19.  DAF  DRO: LSH22JUN.  DAT  (7) 

♦RUN  TLOSS 17HH  <8> 

This  run  is  called  the  LSH22JUN  run.  'LSH'  stands  for 
LOSS  SEVENTEEN  H.  '22 JUN '  indicates  that  this  is  the 

22nd  run  in  June. 
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log 

The 

The 

The 

The 

The 
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first  two  lines  in  this  command  file  put  the 
file  on  the  disk  instead  of  the  printer, 
third  line  puts  unit  7  output  into  file  LSH22JUN.  PLT. 
fourth  line  puts  unit  6  output  into  file  LSH22JUN.  OUT. 
fifth  line  draws  unit  S  input  from  SIMPLE19.  DAF. 
sixth  line  copies  TL0SS17HH.  LIS  to  LSH22JUN.  LIS. 
seventh  line  copies  SIMPLE19.  DAF  to  LSH22JUN.DAT. 
eighth  line  runs  program  TL0SS17HH. 


A  command  file  such  as  this  yields  a  file  groupi 
each  file  having  the  same  name  but  a  different 
extension.  These  files  are  then  treate-d 
collectively.  They  are  printed  together  and  stored 
on  tape  together.  Such  a  group  is.  for  example: 


LSH22JUN.  OAT 
LSH22JUN.  LIS 
LSH22JUN.  LOO 
LSH22JUN.  OUT 
LSH22JUN.  PLT 


a  copy  of  the  data  file 
the  program  listing  file 
the  log  file 
the  output  file 
the  numerical  output  file 


The  user  normally  gets  a  complete  printout  of  the 
run  by  typing  the  command: 


PRINT  LSH22JUN.  *» * 


on  the  terminal. 


Copying  Runs  to  Tape  (  Archiving  ) 

Periodical ly<  it  is  necessary  to  clear  the  old  runs 
off  the  disk  and  put  them  on  tape.  This  is  a  very  simple 
operation  and  proceeds  as  follows: 


a.  Initialize  a  tape.  Oive  it  a  name  that  is  somehow 
indicative  of  the  present  month.  For  example: 


INITIALIZE/DENSITY-1600  MT1:  JUNJUL 


This  tape  would  have  runs  made  in  June  and  July. 

This  step  is  not  necessary  if  the  user  is  adding  runs 
to  a  tape  that  is  already  partly  filled.  In  such  a 

-  .  -  - -  -  —  —  *£ze! 


case,  be  careful  not  to  initial 
b.  Mount  the  tape  with  the  MOUNT  command. 
MOUNT  MT1 :  JUNJUL 


This  is  not  the  same  as  physically  mounting  the  tape 
on  the  tape  drive,  which  is  assumed  to  have  taken 
place  previously. 


c.  Copy  all  the  run  files  from  disk  to  tape  with  the 
COPY  command.  For  example  in  this  case: 


COPY  *JUN.  +>  *  MT1 : 
COPY  *JUL.  *«  *  MT1 : 


would  copy  all  the  June  and  all  the  July  runs  onto 
the  tape  mounted  on  tape  drive  1. 
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d.  For  future  reference,  it  is  good  to  get  a  list  of 
names  of  the  files  on  the  tape.  This  is  accomplished 
as  follows: 


DIRECTORY/PRINT  MT1 : 


The  printout  resulting  from  this  command  does  not 
have  the  tape  name  on  it.  so  the  user  should  write 
this  on  immediately. 

If.  for  some  reason,  this  printout  is  lost,  the  tape 
can  be  remounted  and  the  directory  can  be  taken 
again.  But  this  is  no  reason  to  be  careless  with 
tne  original  directory  printout. 


Dismount  the  tap*  with  the  dismount  command: 
DISMOUNT  MT1 : 


f .  Deist*  the  files  on  the  disk  that  have  just  been 
copied  to  the  tap*.  In  this  case: 

DELETE  *JUN.  *;  * 

DELETE  *JUL.  ** * 

4.  Distributing  the  Program  by  Tape 

The  program  and  test  cases  are  distributed  using 
tapes  with  the  same  format  as  the  archival  tapes.  This 
is  the  Files-11  format*  peculiar  to  Digital  Equipment 
Corporation.  The  tapes  from  this  system  have  a  density 
of  1600  bpi<  and  a  blocksiz*  of  2048  bytes. 


If  any  further  questions  arise*  feel  free  to  call: 

Dr.  H.  Kutschal*  (914)  359-2900  x382 

Mr.  L.  Carroll  (914)  359-2900  x221 
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TEST  CASE  DATA 
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DATA  FOR  TL0SS17HH 
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TEST  CASE  DATA 
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TL0SS51HH  USER'S  MANUAL 


Introduction 

TL0SS51HH.  our  most  complete  transmission  loss  (  hereafter 
abbreviated  as  TL  >  program.  utas  completed  in  February  of  1983. 

It  is  an  extension  of  TL0SS17HH.  which  is  an  extension  of  TLQSSIOHH. 


the  interim  standard  TL  program 
TL0SS91HH  has  these  features  a 


and  options  for  the  user: 

1.  A  choice  of  computational  methods; 

a.  The  FFP  method 

b.  The  Normal  Mode  method 

If  the  Normal  Mode  method  is  chosen.  the  user  also 
has  the  option  of  including  either  of  two  branch  cut 
integrals. 

2.  Optional  horizontal  or  vertical  displacement  in  addition 
to  the  usual  pressure  values  far  TL. 

3. 


A  comprehensive  selection  of  modelling  parameters. 
These  include: 

a.  A  detailed  sound  velocity  profile 
Choice  of  source  depth 
Choice  of  detector  depth 
Choice  of  frequency 

Choice  of  wave  number  integrand  limits 
Optional  beaming,  as  it  is  affected  by 
integrand  tapering  (  in  the  FFP  mode  > 
Optional  surface  roughness 

Choice  of  either  fluid  or  vacuum  upper  half 
space 

Possibility  of  the  source  and  detector  in 
the  fluid  upper  half  space 
Optional  internal  roughness 
Choice  of  initial  range,  range  step,  and 
range  count  (  for  the  Normal  Mode  method  > 


b. 

c. 

d. 

e. 

f. 


l. 


Output. 

a. 

b. 

c. 

d. 


5. 


This  includes: 

A  printout  of  all  input 
The  TL  plot 

The  wave  integrand  plot  of  the  FFP;  either 
short  or  long. 

A  table  of  the  denominator  values  for  the 
Normal  Mode  method 

e.  A  table  of  roots  and  their  corresponding 
phase  velocities,  slopes,  and  excitation 
functions  for  the  Normal  Mode  method 

Note:  All  short  plots  except  the  wave  integrands 
are  repeated,  the  repeated  plot  being  a  smoothed 
copy  of  the  first  and  having  the  same  bounds  as 
the  first. 

Additional  Optional  Tables  and  Plots.  These  include: 
a.  The  TL  plot  in  nautical  miles,  if  desired  and 
if  the  range  is  sufficiently  long.  The  original 
TL  plot  is  in  kilometers. 

The  denominator  plot  for  the  Normal  Mode  method. 

Modal  TL  plots.  Branch  TL  plots,  and  Total  TL 
plots,  if  a  branch  cut  option  is  chosen. 

A  long  printer-plot  of  the  sound  velocity  profile  after 
the  profile  has  been  broken  into  equal  increment  depths 
and  also  a  table  of  this  expanded  sound  velocity  profile 
Eigen  value  printer-plots  far  each  mode,  if  there  are  mode’ 

Numerical  output  of  the  TL  for  further  processing  or  plotting. 

put  out  in  a  separate  file,  if  there  are  ai 


b. 

c. 

d. 


since  it  went  out 


7. 


Also,  eigen  values  are  , 

Formerly,  this  was  called  punched  output 
to  the  card  punch. 

A  profile  check  for  multiple  ducts. 

This  check  is  made  for  the  Normal  Mode  method. 

If  multiple  ducts  are  present,  special  care  is  taken 
to  locate  those  'hidden'  roots  which  are  only  found 
by  reducing  the  step  size  of  the  denominator  function. 


TL0SS51HH  USER'S  MANUAL 


Input  Variables  and  Formats 

Program  input  for  TLOSS31HH  includes  both  the  physical  model 
and  program  control.  This  section  of  the  manual  will  list 
and  explain  each  line  of  the  input. 

1.  Entry  1>  one  line 

MMAXF.  NCARD  variable  names 

(  213  >  format 

MMAXF  is  a  control  variable  which  is  used  to  read 
in  the  correct  number  of  points  in  the  sound  velocity 
profile.  <  see  2.  -  next  Entry  ) 

NCARD  is  a  control  variable  that*  when  set  non-zero 
allows  the  program  to  put  out  TL  data  on  unit  7  in 
a  form  which  is  readable  as  input  bu  another  program. 

A  zero  value  will  prevent  any  output  on  unit  7.  This 
is  normally  used  to  make  graphic  plots.  (  Originally, 
on  the  IBM  360/91.  unit  7  was  the  punched  card  unit, 
hence  the  name  NCARD.  ) 

2.  Entry  2.  MMAXF  lines  <  max.  of  49  lines  > 

((  V(M.  N).  N-1,7  ),  M*1 ,  MMAXF  )  variable  name 

<  7D10.  0  )  format 

V  is  an  array  that  contains  the  sound  velocity  profile. 
For  each  point  in  the  profile,  and  on  each  line. 

(  one  point  per  line  >.  the  program  expects: 

a.  a  depth  value,  in  meters 

b.  a  sound  (  pressure  )  velocity  value,  in  meters  per  second 

c.  a  sound  (  shear  )  velocity  value,  in  meters  per  second, 

if  applicable 

d.  a  corresponding  attenuation  value  for  pressure 

e.  a  corresponding  attenuation  value  for  shear, 

if  applicable 

f.  a  density  value,  in  grams  per  c. c. .  or  specific  density 

g.  a  roughness  value,  if  applicable 

Note:  No  depth  value  can  be  less  than  a  previous 

depth  value.  The  program  will  stop  if  this  is  so. 

3.  Entry  3.  one  line 

KSF.  LRF.  MM.  LS.  KKK,  IBEAM.  IMAX,  FREQ.  PO.  C,  CMAX 

variable  names 

(  715.  4D10.  0  )  format 

KSF  is  a  pointer  to  the  source  depth  value. 

It  points  to  one  of  the  points  in  the  sound  velocity 
profile;  the  one  which  is  the  source  depth  value. 

LRF  is  a  pointer  to  the  detector  (  or  receiver  ) 
depth  value.  It's  function  is  similar  to  KSF. 

MM  is  a  TL  plot  pruning  factor,  with  a  minimum 
acceptable  value  of  2.  MM  is  applicable  to  FFP  only. 

If  an  MM  value  of  N  is  read  in.  then  the  first  1/N 
of  the  plot  will  be  plotted.  This  is  useful  either  to 
get  a  more  detailed  look  at  the  output  or  to  prune 
away  that  part  of  the  plot  affected  by  range  aliasing. 

LS  is  the  number  of  points  used  in  the  wave 
integration.  It  also  affects  the  range  of  the  TL  plot, 
for  the  FFP  method. 

KKK  is  a  control  variable  that  indicates: 

a.  If  zero,  this  is  the  last  case. 

b.  If  negative,  the  next  case  starts  with  Entry  1. 

<  new  profile  ). 

c.  If  positive,  the  next  case  starts  with  Entry  3. 

IBEAM  is  a  control  variable  that,  when  set  non¬ 
zero.  requires  the  program  to  taper  the  wave  integrand 
at  each  end. 
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IMAX  is  the  range  count  or  the  number  of  points  in 
the  TL  plot  for  a  Normal  Mode  case. 

FREQ  is  the  frequency  for  this  case,  in  Hertz. 

PO  is  the  source  level.  For  TL  PO  »  1.0 

C  is  the  lower  bound  sound  velocity  in  the  wave 
integration,  in  meters  per  second.  For  the  Normal  Mode  method, 
if  a  zero  is  entered,  a  default  value  will  be  found 
from  the  sound  velocity  profile. 

CMAX  is  the  upper  bound  sound  velocity  in  the  wave 
integration,  in  meters  per  second.  If  a  zero  is  entered 
here,  the  wave  number  integration  will  start  from  zero 
for  the  FFP  method,  or.  for  the  Normal  Mode  method,  a 
default  value  will  be  found  from  the  sound  velocity 
profile. 

4.  Entry  4.  one  line 

REFLTO.  ALPHAO,  ATTENO.  RHOO,  H.  INAUT.  LONG.  KHPV.  IBR. 

IDPL  variable  names 

(  5010.  0.  515  )  format 

REFLTO  is  the  RMS  surface  roughness  value,  in  me-ters. 

If  this  value  is  non-zero,  the  program  explicitly  sets 
the  next  three  values  to  zero. 

ALPHAO  is.  if  present,  the  sound  velocity  in  the 
fluid  upper  half  space,  in  meters  per  second.  This  is 
the  variable  that  determines  whether  the  present  case  has 
either  a  fluid  or  vacuum  upper  half  space. 

ATTENO  is.  if  present,  the  attenuation  in  the 
fluid  upper  half  space. 

RHOO  is.  if  present,  the  density  of  the  fluid 
upper  half  space,  in  grams  per  c. c. 

H  is.  if  present,  the  distance  between  the  source 
and  detector  in  the  fluid  upper  half  space.  To  use  this 
option: 

a.  The  first  two  sound  velocity  depths  must  be  zero. 

b.  There  must  be  some  variation  in  the  pressure 
sound  velocity  profile. 

c.  There  must  not  be  any  shear  sound  velocity. 

d.  ALPHAO  and  H  must  be  present. 

INAUT  is  a  control  variable  that,  when  set  non¬ 
zero.  allows  the  plotting  of  the  TL  plot  in  nautical 
miles  if  the  range  is  greater  than  100  kilometers. 

This  is  for  the  FFP  method  only. 

LONG  is  a  control  variable  that,  when  set  non¬ 
zero.  switches  the  wave  integrand  plot  from  short  to 
long.  This  is  used  when  it  ts  necessary  to  examine 
each  point  of  the  wave  integrand.  FFP  only. 

KHPV  is  the  control  variable  that  determines 
whether  output  will  be: 

a.  pressure,  if  zero. 

b.  vertical  displacement,  if  positive. 

c.  horizontal  displacement,  if  negative. 

IBR  is  the  control  variable  that  determines  if 
a  branch  cut  integral  will  be  included  with  the  Normal 
Mode  method. 

a.  If  zero.  no. 

b.  If  positive,  yes  <  exact  ) 

c.  If  negative,  yes  (  FFP  ) 

IDPL  is  the  control  variable  that,  when  set.  allows  the 
denominator  function  to  be  plotted.  Normal  Mode  method  only. 
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5.  Entry  5.  ont  line 

RA.  DLRA.  THLAYR,  IEXP  variable  names 

(  3010.  0.  19  )  format 

RA  is  the  initial  range  of  the  TL  plot,  in  meters. 

DLRA  is  the  range  step  of  the  TL  plot#  in  meters. 

These  two  variables  are  read  in  for  the  Normal  Mode 
method,  when  IBR  is  zero  or  positive,  only.  Since  the 
TL  plot  is  in  kilometers,  it  is  very  important  to  be 
careful  with  the  units  of  these  two  variables. 

THLAYR  is  the  desired  layer  thickness  for  the  sound 
velocity  profile  which  is  used  if  IEXP  is  set  non-zero. 

New  interface  depths  are  chosen  by  interpolation  every 
THLAYR  meters  and  the  original  data  points  are  retained. 
This  is  valid  for  FFP  and  Modes. 

IEXP  is  a  control  flag  which,  when  set  non-zero, 
directs  the  program  to  split  up  the  input  sound  velocity 
profile  and  calculate  TL  on  this  new  profile.  If  this  is  a 
normal  mode  run.  and  if  there  are  modes  then  plot  out  on 
the  printer  eigen  value  plots  for  each  mode.  Also,  eigen 
values  vs.  depth  are  put  out  on  Unit  S  to  a  file  for 
further  processing  or  plotting. 


Note: 

The  program  determines  the  computational  method. 

(  whether  FFr  or  Normal  Mode  ).  by  examining  all  of 
the  attenuation  values  and  the  surface  roughness  value. 

If  they  are  all  zero,  the  Normal  Mode  method  is  selected. 
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C. 


TLQSS51HH  USER'S  MANUAL 
Operating  Procedures 


The  primary  computer  at  Lamont-Doherty  Geological  Observatory 
is  a  VAX  ll/780<  which  runs  under  the  VMS  operating  system. 

Listed  below  are  four  distinctly  different  operations  that 
are  used  to  process  TL0SS31HH. 


1.  Compiling-Linking 

a.  Compiling  <  FORT  command  ) 

Compiling  the  TL0SS91HH  program  on  the  VAX  11/780 
has  always  been  done  with  the  following  switches 
on: 

1.  /CHECK  Forbids  array  pointer  overflow. 

2.  /N0I4  Integer  variables  are  set  to  2  bytes. 

3.  /NQF77  Where  the  compiler  needs  to  choose 
between  Fortran  77  and  Fortran  4+.  the  latter 
is  chosen. 

4.  /LIST  A  compiler  listing  is  made  and 
stored  on  the  disk  as  TL0SS51HH.  LIS. 


b.  Linking  (  LINK  command  ) 

Linking  is  straight-f orward.  No  extra  routines 
are  linked  with  TLQSS51HH. 

2.  Running 

TL0SS31HH  is  normally  run  through  a  batch  stream* 
in  particular  the  SLOW  batch  stream*  since  this  is  the 
cheapest  route*  especially  at  night.  However*  there  is 
another  batch  stream  with  higher  priority  (  and  higher 
cost  >  that  can  be  used  for  quicker  results*  the 
SYS4BATCH  batch  stream.  In  order  to  use  either  of  these 
two  batch  streams: 

a.  The  executable  file  (  TL0SS51HH.EXE  >  must  be  ready 
to  use.  The  LINK  operation  satisfies  this  requirement. 

b.  The  data  file  must  also  be  ready  to  use.  This  is 
normally  prepared  by  the  user  using  the  editor* 

<  and  the  TL0SS31HH  USER'S  MANUAL,  if  necessary  > 
at  a  terminal. 

c.  Submit  the  command  file*  TL0SSS1HH.COM*  to  the 
batch  stream  (  or  batch  queue  )  using  a  command 
such  as  one  of  these: 

SUBM I T /N AME»LFQ22FEB / QUEUE “SLOW/ AFTER*20 : 00  TLQSS31HH 
SLOW  queue*  Job  starts  at  8: 00  P.  M. 

SUBMIT/NAME •LFQ22FEB/QUEUE*SL0W  TL0SSS1HH 
SLOW  queue*  Job  starts  now 


SU8MIT/NAME«LF022FEB  TL0SSS1HH 

SYS48ATCH  queue  (  the  'regular'  queue  )* 

Job  starts  now 

d.  The  following  command  file*  also  prepared  using 
the  editor*  is  a  typical  file  submitted*  using  one 
of  the  previous  submit  commands. 
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4 ASSIGN  LF022FEB  SYS4PRINT  (1) 

4 ASSIGN  LF022FEB  SYS40UTPUT  (2> 

4 ASSIGN  DRO: LF022FEB.  PLT  F0RO07  (3) 

4 ASSIGN  DRO: LF022FEB.  OUT  FOR 006  (4) 

4 ASSIGN  DRO: SIMPLE19.  DAF  F0R005  (S) 

tTYPP  Tl  f  A) 

4C0PY  DRO: SIMPLE!?.  DAF  DRO: LF022FEB.  DAT  (7) 
4RUN  TL0SS91HH  (8) 

This  run  is  called  the  LF022FEB  run.  'LFO'  stands  for 
LOSS  FIFTY  ONE  H.  '22FEB '  indicates  that  this  is  the 
22nd  run  in  February. 


The  first  two  lines  in  this  command  file  put  the 
log  file  on  the  disk  instead  of  the  printer. 

The  third  line  puts  unit  7  output  into  file  LF022FEB.  PLT. 
The  fourth  line  puts  unit  6  output  into  file  LFQ22FEB.  OUT. 
The  fifth  line  draws  unit  5  input  from  SIMPLE19.  DAF. 

The  sixth  line  copies  TLOSS51HH.  LIS  to  LF022FEB.  LIS. 

The  seventh  line  copies  SIMPLE19.  DAF  to  LF022FEB.DAT. 

The  eighth  line  runs  program  TLQSS31HH. 

A  command  file  such  as  this  yields  a  file  groups 
each  file  having  the  same  name  but  a  different 
extension.  These  files  are  then  treated 
collectively.  They  are  printed  together  and  stored 
on  tape  together.  Such  a  group  is.  for  example: 

LF022FEB.DAT  a  copy  of  the  data  file 

LF022FEB.  LIS  the  program  listing  file 

LF022FEB.LOG  the  log  file 

LF022FEB.  OUT  the  output  file 

LF022FEB.  PLT  the  numerical  output  file 

The  user  normally  gets  a  complete  printout  of  the 
run  by  typing  the  command: 

PRINT  LF022FEB.  *»  ■* 

on  the  terminal. 

Copying  Runs  to  Tape  <  Archiving  > 

Periodically,  it  is  necessary  to  clear  the  old  runs 
off  the  disk  and  put  them  on  tape.  This  is  a  very  simple 
operation  and  proceeds  as  follows: 

a.  Initialize  a  tape.  Give  it  a  name  that  is  somehow 
indicative  of  the  present  month.  For  example: 

INITIALIZE/DENSITY*lAOO  MT1 :  FEBMAR 

This  tape  would  have  runs  made  in  February  and  March. 

This  step  is  not  necessary  if  the  user  is  adding  runs 
to  a  tape  that  is  already  partly  filled.  In  such  a 
case,  be  careful  not  to  initialize! 

b.  Mount  the  tape  with  the  MOUNT  command. 

MOUNT  MT1 :  FEBMAR 

This  is  not  the  same  as  physically  mounting  the  tape 
on  the  tape  drive,  which  is  assumed  to  have  taken 
place  previously. 

c.  Copy  all  the  run  files  from  disk  to  tape  with  the 
COPY  command.  For  example  in  this  case: 

COPY  *FEB.  *.  *  MT1 : 

COPY  *MAR.  *i  *  MT1 : 

would  copy  all  the  February  and  all  the  March  runs  onto 
the  tape  mounted  on  tape  drive  1. 

d.  For  future  reference,  it  is  good  to  get  a  list  of 
names  of  the  files  on  the  tape.  This  is  accomplished 
as  follows: 

DIRECTORY/PRINT  MTl : 

The  printout  resulting  from  this  command  does  not 
have  the  tape  name  on  it.  so  the  user  should  write 
this  on  immediately. 

If.  for  some  reason,  this  printout  is  lost,  the  tape 
can  be  remounted  and  the  directory  can  be  taken 
again.  But  this  is  no  reason  to  be  careless  with 
the  original  directory  printout. 


•.  Dismount  the  tap#  with  the  dismount  commend: 
DISMOUNT  MT1: 

f .  Delete  the  files  on  the  disk  that  have  just  been 
copied  to  the  tape.  In  this  case: 

DELETE  *FEB.  *:  * 

DELETE  *MAR.  *»  * 


Distributing  the  Program  by  Tape 

The  program  and  test  cases  are  distributed  using 
tapes  witn  the  same  format  as  the  archival  tapes.  This 
is  the  Files-11  format*  peculiar  to  Digital  Equipment 
Corporation.  The  tapes  from  this  system  have  a  density 
of  1600  bpii  and  a  blocksize  of  204s  bytes. 


bytes. 


If  any  further  questions  arise*  feel  free  to  call 
Dr.  H.  Kutschale  (914)  359-2900  x3 


Mr.  L.  Carroll 


(914)  359-2900  x382 
(914)  359-2900  x221 
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LF002MAR.  DAT;  2 


06 


TEST  CASE  DATA 


0.  00 

1500.  0 

. 0000000 

50.  00 

1500.  0 

.  0000000 

100.  0 

1500.  0 

.  0000000 

100.  0 

3000.  0 

1400.  0 

.  00006286 

125.  0 

3000.  0 

1400.  0 

. 00006286 

125.  0 

3000.  0 

1400.  0 

.  00006286 

2 

3  02  1024 

00  000 

1  020 

0.  0 

1500.  0 

. 0000000 

1.  00  i 

00. 

00.  0 

00.  000 

1 

6- APR- 1993  15:30 
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.  0001281 
.  0001281 
.  0001281 
.  0  10000 
0000.  0 


TL0SS51HH 
1.  00 
1.  00 
1.  00 
2.  50 
2.  50 
2.  SO 

1000.  0 
0  0 


8000.  0 
00  O 


LF003MAR.  DAT;  3 


6-APR-1933  15:30 
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04  1  0 

00.  0 
25.  0 
50.  0 
50.  0 
2 

0.  0 

0200.  0 


TEST  CASE  DATA 


1500.0  0.000 

1525.0  0.000 

1550.  0  0.  000 

1600.  0  0.  000 

2  2  1024  00 

00.  0  0.  00 

200.  0  0.  00 


DATA  FOR  TL0SS51HH 

0.  00000  0.  00000  1.  00 

0.00000  0.00000  1.00 

0.00000  0.00000  1.00 

0.  00000  0.  00000  2.  00 

0  100  110.0  10C00.  0  1 500.. 001 

00.  00  0.  000  0  0  0 

1  RANGE  &  DLRANGE  IN  METERS, 


1599.  999 
-1 

OUTPUT  IN 


LC«d 


LF004MAR.  OAT;  2 


6-AFR-1983  i5:30 


TEST  CASE  DATA 


08  1  0 
000.  0 

1336.  3 

. 0000100 

1.  0 

132.  4 

1539.  2 

.  0000100 

1.  0 

233.  9 

1330.  0 

. 0000100 

1.  0 

406.  3 

1501.  1 

. 0000100 

1.  0 

863.  3 

1475.  0 

.  0000100 

1.  0 

1013.  9 

1471.  9 

.  0000100 

1.  0 

3387.  9 

1349.  6 

.  0000100 

1.  0 

3387.  9 

1349.  6 

.  0000100 

1.  9 

3  3 

2  312 

-1  00 

1  23. 

10000.  0 

00.  0 

0000.  0 

0.  00000 

0.  0 

000.  0  0 

00.  0 

0000.  0 

1470.  0 
0 


1349.  4 


08  1  0 

000.  0 

1336.  3 

. 0000000 

1.  0 

132.  4 

1339-.  2 

.  0000000 

1.  0 

233.  9 

1 530.  0 

.  0000000 

1.  0 

406.  3 

1301.  1 

.  0000000 

1.  0 

863.  3 

1475.  0 

. 0000000 

1.  0 

1015.  9 

1471.  9 

.  0000000 

1.  0 

3587.  9 

1349.  6 

.  0000000 

1.  0 

3587.  9 

1549.  6 

.  0000000 

1.  9 

3  5 

2  256 

00 

00 

1  25. 

10000.  0 

00.  0 

0000.  0 

0.  00000 

0.  0 

000.  0  < 

00.  0 

0000.  0 

1470.  0 
0 


1549.  4 


LF003MAR.  DATi  2  6-APR-1983  15: 30 

TEST  CASE  DATA 


OS  1  0 

000.  0 

1336.  3 

. 0000100 

1.  0 

132.  4 

1339.  2 

. 0000100 

1.  0 

233.  9 

1330.  0 

. 0000100 

1.  0 

406.  3 

1301.  1 

. 0000100 

1.  0 

863.  3 

1473.  0 

.  0000100 

1.  0 

1013.  9 

1471.  9 

. 0000100 

1.  0 

3387.  9 

1349.  6 

. 0000100 

1.  0 

3387.  9 

1349.  6 

. 0000100 

1.  9 

3  3  2  312  -1  CO  1  23.  10000.  0 

00.  0  0000.  0  0.  00000  0.  0  000.  0  0 

00.  O  0000.  0  030.  00  1 


08  1  0 
000.  0 

1336.  3 

. 0000000 

1.  0 

132.  4 

1339.  2 

.  0000000 

1.  0 

233.  9 

1330. 0 

. 0000000 

1.  0 

406.  3 

1301.  1 

.  0000000 

1.  0 

863.  3 

1473.  0 

.  0000000 

1.  0 

1013.  9 

1471.  9 

.  0000000 

1.  0 

3387.  9 

1349.  6 

.  0000000 

1.  0 

3387.  9 

1349.  6 

.  0000000 

1.  9 

3  3 

2  256 

00  00 

1  23. 

10000.  0 

00.  0 

0000.  0 

0.  00000 

0.  0 

000.  0 

00.  0 

0000.  0 

050.  00 

1 

INTEGRAND 


INTEGRAND  PLOT 


900. 


.640  .660  .680 


WAVE  SLOWNESS 

1  25.0  HZ  L51H  SM  0 


, 

•J'’- 

77T\<;\  .-.V.KT^ 

-’■*  -*  *'  '  *  ^ 

Tpr  Tr  '*T"r  T»5 

■77TV7-'77.~t 
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LF006MAR. 

DAT;  3 

6-APR-1983  15:30 

TEST 

CASE  DATA 

?:v 

06  1 

DATA  FOR 

TLQSS51HH 

0.  00 

1500.  0 

.  0000000 

1. 

00 

.v' 

50.  00 

1500.  0 

.  0000000 

1. 

00 

■V-1 

100.  0 

1500.  0 

.  0000000 

1. 

00 

100.  0 

2800.  0 

1200.  0 

.  00006286 

. 0001281 

2. 

50 

.■  1:  ■ 

125.  0 

3000.  0 

1400.  0 

. 00006286 

.  0001281 

2. 

50 

125.  0 

3000.  0 

1400.  0 

. 00006286 

.  0001281 

2. 

50 

r : 

2 

3  02  1024 

00  000 

1  020.  0  10000. 

1000.  0 

0.  0 

1500.  0 

.  0000000 

1.  00 

0000.  0 

0  0 

00. 

00.  0 

00.  000 

1 
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IMPULS1 IB  USER'S  MANUAL 


Introduction 

IMPULS11B/  our  streamlined  FFP  impulse  program*  was  completed  in 
October  of  1982.  It  is  an  extension  of  IMPULSIOBB.  the  former  interim 
standard  impulse  program. 

IMPULS11B  has  these  features  and  options  for  the  user: 

1.  The  FFP  method  of  computation 


2. 


3. 


5. 


f 


A  comprehensive  selection  of  modelling  parameters. 
These  include: 

A  detailed  sound  velocity  profile 
Choice  of  source  depth 
Choice  of  detector  depth 

Choice  of  frequency  band  and  source  frequency 
Choice  of  wave  number  integrand  limits 
Optional  beaming*  as  it  is  affected  by 
integrand  tapering 
Optional  surface  roughness 

Choice  of  either  fluid  or  vacuum  upper  half 
space 

Possibility  of  the  source  and  detector 
the  fluid  upper  half  space 
Optional  internal  roughness 
Choice  of  initial  range*  range  spacing* 
range  count 


f! 


i 


m 


and 


Output. 

a. 

b. 


This  includes: 

A  printout  of  all  input 
The  pulse  signal  plots 


Additional  Optional  Tables  and  Plots.  These  include: 


The  wave  integrand  plots  of  the  FFP  for  certain 
frequencies 


Numerical  output  of  the  program  for  further  processing  or  | 
Formerly,  this  was  called  punched  output,  since  it  went  out 
the  card  punch. 


punch 


lotting^ 


IMPULS11B  USER'S  MANUAL 


Input  Variable*  and  Formats 

Program  input  for  IMPULS11B  includes  both  the  physical  model 
and  program  control.  This  section  of  the  manual  will  list 
and  explain  each  line  of  the  input. 


1.  Entru  1.  one  line 

MMAXF,  NOUTPT,  NREVRS,  NOISE.  NCARD,  NPLOT .  NREAD 
(  713  ) 


variable  names 
format 


MMAXF  is  a  control  variable  which  is  used  to  read 
in  the  correct  number  of  points  in  the  sound  velocity 
profile.  (  see  Entry  2.  -  next  Entry  ) 

NOUTPT  is  the  number  of  output  signal  plots. 

NREVRS  is  both  a  flag  and  a  pointer.  If  non-zero,  it 
signals  that  there  will  be  pulse  reversal,  and  it  points 
to  the  output  signal  which  is  reversed  and  used  as  a 
source  function. 

NOISE  is  a  control  flag  that,  when  set  to  zero. 

(  be  careful!  >  zeroes  out  the  output  signal  which  is 
outside  the  UMIN  to  UMAX  interval.  This  eliminates  the 
noise  outside  the  specified  time  interval,  which  can 
be  severe  at  long  ranges. 

NCARD  is  a  control  variable  that,  when  set  non-zero, 
allows  the  program  to  put  out  signal  data  on  unit  7  in  a  form  which 
is  readable  as  input  by  another  program.  A  zero  value  will  prevent 
anu  output  on  unit  7.  This  is  normally  used  to  make  graphic  plots. 

(  Originally,  on  the  IBM  360/91.  unit  7  was  the  punched  card  unit, 
hence  the  name  NCARD.  ) 

NPLOT  is  a  control  variable  that,  when  set  non-zero, 
allows  the  program  to  put  out  printer  plots  of  the 
output  signals. 

NREAD  is  a  control  variable  that,  when  set  non-zero, 
directs  the  program  to  interpret  Entry  2  data  (  see  below  > 
as  thickness  input  instead  of  the  usual  depth  input. 


Entry  2.  MMAXF  lines  <  max.  of 
((  V(M,  N) »  N-1,7  ).  M-l,  MMAXF  ) 
(  7D10.  0  ) 


variable  name 
format 


V  is  an  array  that  contains  the  sound  velocity  profile. 

For  each  point  in  the  profile,  and  on  each  line. 

<  one  point  per  line  ).  the  program  expects: 

a.  a  depth  value,  in  meters,  except  when  NREAD  is  set  non-zero, 
in  which  case  the  program  expects  a  thickness  value,  in  maters 

b.  a  sound  (  pressure  >  velocity  value,  in  meters  per  second 

c.  a  sound  (  shear  )  velocity  value,  in  meters  per  second, 

if  applicable 

d.  a  corresponding  attenuation  value  for  pressure 

e.  a  corresponding  attenuation  value  for  shear, 

if  applicable 

f.  a  density  value,  in  grams  per  c. c. .  or  specific  density 

g.  a  roughness  value,  if  applicable 

Note:  No  depth  value  can  be  less  than  a  previous 

depth  value. 

Entry  3.  one  line 

KSF.  LRF,  IRES.  LS.  KKK.  IBEAM.  NW,  RSPACE.  PO.  C.  CO 

variable  names 

<  715.  4D10.  0  )  format 

KSF  is  a  pointer  to  the  source  depth  value. 

It  points  to  one  of  the  points  in  the  sound  velocity 
profilei  the  one  which  is  the  source  depth  value. 

LRF  is  a  pointer  to  the  detector  (  or  receiver  ) 
depth  value.  It's  function  is  similar  to  KSF. 

IRES  is  a  control  variable  that,  when  set  non-zero, 
requires  the  program  to  taper  the  response  function 
at  each  end. 


.•  V.  V.  .\ 


wwinw 


LS  is  the  number  of  points  used  in  the  wave  integration. 

KKK  is  a  control  variable  that  indicates: 

a.  If  0  or  2*  this  is  the  last  case. 

b.  If  -1  or  -2*  the  next  case  starts  with  Entry  1« 

<  new  profile  ). 

c.  If  >1  or  +3*  the  next  case  starts  with  Entru  3. 

The  second  value  in  each  case  (  2.  -2,  +3  )»  will  direct 
the  program  to  take  each  signal  and  compress  it  back  to 
its  source.  A  plot  of  this  compressed  pulse  follows 
each  output  signal. 

IBEAM  is  a  control  variable  that*  when  set  non- 
zero*  requires  the  program  to  taper  the  wave  integrand 
at  each  end. 

NW  is  a  flag  that*  when  set  non-zero*  directs  the 
program  to  print  out  selected  integrand  plots. 

RSPACE  is  the  range  spacing  between  output  signals* 
in  kilometers.  However*  the  program  can  only  approximate 
this  value. 

PO  is  the  source  level.  For  this  program  we  normally 
use  a  value  of  10000. 

C  is  the  lower  bound  sound  velocity  in  the  wave 
integration*  in  meters  per  second. 

CO  is  the  upper  bound  sound  velocity  in  the  wave 
integration*  in  meters  per  second. 

Entry  4.  one  line 

FREQo7  NX,  INTF,  FREQM,  DT,  RO,  FREQS.  UMAX.  UMIN 

variable  name 

(  F10.  0*  215*  6F10.  0  >  format 

FREQO  is  the  lower  bound  of  the  frequency  integration, 
in  Hertz. 

NX  is  the  number  of  FFT  points  and  the  number  of 
points  in  each  output  signal. 

INTF  is  a  variable  that  specifies  how  often*  within  the 
frequency  integration,  the  wave  integrand  will  be  plotted  on  the 
printer*  (  every  INTF'th  frequency  ).  Enter  NX+1  for  no  wave 
integrand  plots. 

FREQM  is  the  upper  bound  of  the  frequency  integration* 
in  Hertz. 

OT  is  the  time  difference  between  the  successive  points  in 
the  output  signals*  in  secs. 

RO  is  the  initial  range  or  the  range  of  the  first 
output  signal*  in  kilometers. 

FREQS  is  the  sound  source  frequency*  in  Hertz. 

UMAX  is  a  sound  velocity  (  meters/sec.  >  that  is  used  to 
specify  the  lower  bound  of  tne  output  signal  if  some  of  the  signal 
is  considered  noise.  (  See  NOISE*  Entry  1*  also  see  CRED*  Entry  5 

UMIN  is  similar  in  function  to  UMAX,  but  used  to 
specify  the  upper  bound. 


H.  CRED 


Entry  5.  on*  lin* 

REFLTO*  ALPHAO,  ATTENO,  RHOO, 
(  6D10.  0  ) 


variable  names 
format 


REFLTO  is  the  RMS  surface  roughness  value*  in  meters. 
If  this  value  is  non-zero*  the  program  explicitly  set? 
the  next  three  values  to  zero. 


ALPHAO  is<  if  present*  the  sound  velocity  in  the 
fluid  upper  half  space*  in  meters  per  second.  This  is 
the  variable  that  determines  whether  the  present  case  has 
either  a  fluid  or  vacuum  upper  half  space. 

ATTENO  is*  if  present*  the  attenuation  in  the 
fluid  upper  half  space. 

RHOO  is*  if  present*  the  density  of  the  fluid 
upper  half  space*  in  grams  per  c. c. 

H  is*  if  present*  the  distance  between  the  source 
and  detector  in  the  fluid  upper  half  space.  To  use  this 
option: 

a.  The  first  two  sound  velocity  depths  must  be  zero. 

b.  There  must  be  some  variation  in  the  pressure 
sound  velocity  profile. 

c.  There  must  not  be  any  shear  sound  velocity. 

d.  ALPHAO  and  H  must  be  present. 

CRED  is  the  reducing  velocity.  For  each  output  signal* 
the  starting  time  is  offset  from  TO  at  rang*  0  <  point 
of  origin  )  bu  a  time  equal  to  that  rang*  divided  by  CRED. 
This  value  defaults  to  CO*  for  a  blank  or  zero  entry. 


IMPULS1 IB  USER'S  MANUAL 


Operating  Procedures 

The  primary  computer  at  Lamont-Doherty  Geological  Observatory 
is  a  VAX  11/780.  which  runs  under  the  VMS  operating  system. 

Listed  below  are  four  distinctly  different  operations  that 
are  used  to  process  IMPULSllB. 


1. 


Comp i 1 ina-Link ing 

a.  Compiling  (  FORT  command  > 

Compiling  the  IMPULS11S  program  on  the  VAX  11/780 
has  always  been  done  with  the  following  switches 
on: 

1.  /CHECK  Forbids  array  pointer  overflow. 

2.  /N0I4  Integer  variables  are  set  to  2  bytes. 

3.  /NQF77  Where  the  compiler  needs  to  choose 
between  Fortran  77  and  Fortran  4+.  the  latter 
is  chosen. 

4.  /LIST  A  compiler  listing  is  made  and 
stored  on  the  disk  as  IMPULSllB.  LIS. 


b.  Linking  (  LINK  command  > 

Linking  is  straight-forward.  No  extra  routines 
are  linked  with  iliPULSllB. 


2.  Running 

IMPULSllB  is  normally  run  through  a  batch  stream, 
in  particular  the  SLOW  batch  stream,  since  this  is  the 
cheapest  route,  especially  at  night.  However,  there  is 
another  batch  stream  with  higher  priority  (  and  higher 
cost  )  that  can  be  used  for  quicker  results,  the 
SYS*BATCH  batch  stream.  In  order  to  use  either  of  these 
two  batch  streams: 


a.  The  executable  file  (  IMPULS11B.EXE  )  must  be  ready 

to  use.  The  LINK  operation  satisfies  this  requirement. 

b.  The  data  file  must  also  be  ready  to  use.  This  is 
normally  prepared  by  the  user  using  the  editor. 

(  and  the  IMPULS1 IB  USER 'S  MANUAL,  if  necessary  > 
at  a  terminal. 


c.  Submit  the  command  file.  IMPULS11B.COM.  to  the 
batch  stream  <  or  batch  queue  )  using  a  command 
such  as  one  of  these: 

SUBMIT /NAM£*IES 1 50CT /QUEUE“SLOW/ AFTER*20: 00  IMPULSllB 
SLOW  queue.  Job  starts  at  8: 00  P.  M. 

SUBM I T /NAME* I EB 1 9QC T / QUEUE “SLOW  IMPULSllB 
SLOW  queue.  Job  starts  now 

SUBM I T /NAME» I EB 1 SOCT  IMPULSllB 

SYStBATCH  queue  (  the  'regular'  queue  ). 

Job  starts  now 

d.  The  following  command  file,  also  prepared  using 
the  editor,  is  a  typical  file  submitted,  using  one 
of  the  previous  submit  commands. 

•ASSIGN  IEB190CT  SYSBPRINF  (1) 

•ASSIGN  IEB150CT  SYS*OUTPUT  (2) 

•ASSIGN  DRO: IEB150CT.  PLT  F0R007  (3) 

•ASSIGN  DRO: IEB150CT.  OUT  F0R006  (4) 

•ASSIGN  DRO: STANDARD.  DAZ  F0R005  (5) 

•TYPE  IMPULSllB  (6) 

•COPY  DRO: STANDARD.  DAZ  DRO: I EB1 SOCT.  DAT  (7) 

•RUN  IMPULSllB  (8) 

This  run  is  called  the  IEB150CT  run.  'IEB'  stands  for 
IMPULSE  ELEVEN  B.  '150CT'  indicates  that  this  is  the 
19th  run  in  October. 


Ths  first  two  linss  in  this  command  fils  put  ths 
log  fils  on  ths  disk  instaad  of  ths  printsr. 

Ths  third  lins  puts  unit  7  output  into  fils  IEB190CT.  PLT. 
Ths  fourth  lins  puts  unit  6  output  into  fils  IEB150CT.  OU1 
Ths  fifth  lins  draws  unit  9  input  from  STANDARD.  DAZ. 

Ths  sixth  lins  copiss  IMPULS1 la.  LIS  to  IEB190CT.  LIS. 

Ths  savsnth  lins  copiss  STANDARD.  DAZ  to  IEB190CT.DAT. 

Ths  sighth  lins  runs  program  IMPULS11B. 

A  command  fils  such  as  this  yialds  a  fils  group; 
sach  fils  having  ths  sams  nams  but  a  diffsrsnt 
sxtsnsion.  Thsss  filss  ars  than  trsatsd 
collactivaly.  Thay  ars  printsd  togsthsr  and  storsd 
on  taps  togsthsr.  Such  a  group  is*  for  sxampls: 


IEB15QCT.DAT  a  copy  of  ths  data  fils 

XEB190CT.  LIS  ths  program  listing  fils 

1EB19QCT.LOG  ths  log  fils 

IEB190CT.  OUT  ths  output  fils 

1EB19QCT.  PLT  ths  numarical  output  fils 

Ths  ussr  normally  gats  a  complsts  printout  of  ths 
run  by  typing  ths  command: 

PRINT  IEB190CT.  *;# 

on  ths  tarminal. 

Copying  Runs  to  Taps  (  Archiving  > 

Pariodically*  it  is  nscsssary  to  clsar  ths  old  runs 
off  ths  disk  and  put  thsm  on  taps.  This  is  a  vary  simpls 
opsration  and  procssds  as  follows: 

a.  Initializs  a  taps.  Giva  it  a  nams  that  is  somshow 
indicativs  of  ths  prsssnt  month.  For  sxampls: 

INITIALIZE/DENSITY-l&OO  MT1 :  OCTNQV 

This  taps  would  havs  runs  mads  in  Octobar  and  Novsmbsr. 
This  stsp  is  not  nscsssary  if  ths  ussr  is  adding  runs 
to  a  taps  that  is  alrsady  partly  fillsd.  In  such  a 
cass>  bs  carsful  not  to  initializs! 

b.  Mount  ths  taps  with  ths  MOUNT  command. 

MOUNT  MT1 :  OCTNOV 

This  is  not  ths  sams  as  physically  mounting  ths  taps 
on  ths  taps  drivs*  which  is  assumad  to  havs  taksn 
placs  prsviously. 

c.  Copy  all  ths  run  filss  from  disk  to  taps  with  ths 
COPY  command.  For  sxampls  in  this  cass: 

COPY  *OCT.  *•  *  MTl: 

COPY  *NOV.  *j*  MTl: 

would  copy  all  ths  Octobar  and  all  ths  Novsmbsr  runs 
onto  ths  taps  mountsd  on  taps  drivs  1. 

d.  For  futurs  rsfsrsncs*  it  is  good  to  gst  a  list  of 
namss  of  ths  filss  on  ths  taps.  This  is  accomplishsd 
as  follows: 

DIRECTORY/PRINT  MTl: 

Ths  printout  rssulting  from  this  command  doss  not 
havs  ths  taps  nams  on  it*  so  ths  ussr  should  writs 
this  on  immsdiatsly. 

If*  for  soms  rsason*  this  printout  is  lost*  ths  taps 
can  bs  rsmountsd  and  ths  dirsctory  can  bs  taksn 
again.  But  this  is  no  rsason  to  6s  carslsss  with 
ths  original  dirsctory  printout. 
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•.  Dismount  th#  top*  with  th*  dismount  command: 

DISMOUNT  MT1: 

f.  D*l*t*  the  files  on  the  disk  that  have  just  been 
copied  to  th*  tap*.  In  this  case: 

DELETE  *OCT.  *;  * 

DELETE  *NOV.  *i  * 

4.  Distributing  th*  Program  by  Tap* 

Th*  program  and  test  cases  are  distributed  using 
tapes  with  the  same  format  as  th*  archival  tapes.  This 
is  th*  Files-11  format#  peculiar  to  Digital  Equipment 
Corporation.  The  tapes  from  this  system  have  a  density 
of  1600  bp i.  and  a  blocksize  of  2048  bytes. 

If  any  further  questions  arise#  feel  free  to  call: 

Dr.  H.  Kutschal*  (914)  399-2900  x382 

Mr.  L.  Carroll  (914)  359-2900  x221 


IEB040CT.  DAT;  2 


6- APR— 1983  15:29 
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TEST  CASE  DATA 


06  41  0 

111* 

0  0 

DATA  FOR  IMPULS1 10 

0.  00 

1500.  0 

00000150 

1.  00 

50.  00 

1500.  0 

00000150 

1.  00 

100.  0 

1500.  0 

00000150 

1.  00 

too.  0 

3000.  0 

1400.  0 

00000770 

.  00C016 

2.  50 

125.  0 

3000.  0 

1400.  0 

00000770 

. 000016 

2.  50 

125.  0 

3000.  0 

1400.  0 

00000770 

.  000016 

2.  50 

2  3 

00  2048 

00  001 

1  00.  50  10000. 

0.  5000000 

2048  60 

4.  1000000 

0.  010 

05.  90 

2.  0 

0.  0 

1500.0 

.  0000150 

1.  00 

0000.  0 

4000.  0 

1000.  0 

7500.  i 
00 


8000.  0 

1500.  0 
0  1 
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IMPULS22B  USER'S  MANUAL 
Introduction 

IMPULS22B.  one  of  our  streamlined  FFP  impulse  programs,  was  completed 
in  March  of  1983.  It  is  an  extension  of  IMPULS11B.  one  of  the  interim 
standard  impulse  programs. 

IMPULS22B  has  these  features  and  options  for  the  user: 

1. 

a.  The  complete,  unabridaed  version 

of 


The  FFP  method  of  computation 

a.  The  complete,  unabridged  version 

b.  An  abbreviated  version  to  calculate  the  effect 
one  particular  number  of  bounces 

A  comprehensive  selection  of  modelling  parameters. 
These  include: 

a.  A  detailed  sound  velocity  profile 
Choice  of  source  depth 
Choice  of  detector  depth 

Choice  of  frequency  band  and  source  frequency 
Choice  of  wave  number  integrand  limits 
Optional  beaming,  as  it  is  affected  by 
integrand  tapering 
Optional  surface  roughness 

Choice  of  either  fluid  or  vacuum  upper  half 
space 

Possibility  of  the  source  and  detector  in 
the  fluid  upper  half  space 
Optional  internal  roughness 

Choice  of  initial  range,  range  spacing,  and 
range  count 


b. 

c. 

d. 

e. 

f . 

X: 

i. 


3. 


t: 

Output. 


This  includes: 

a.  A  printout  of  all  input 

b.  The  pulse  signal  plots 

Additional  Optional  Tables  and  Plots. 

a 


These  include: 


The  wave  integrand  plots  of  the  FFP  for  certain 
frequencies 

Numerical  output  of  the  program  for  further  processing  or  plotting. 
Formerly,  this  was  called  punched  output,  since  it  went  out  to 
the  card  punch. 


>: 


CT 
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Input  Variables  and  Formats 

Program  input  for  IMPULS22S  includes  both  the  physical  model 

and  program  control.  This  section  of  the  manual  will  list 

and  explain  each  line  of  the  input.  - 

1.  Entry  1,  one  line  5 

MMAXF,  NOUTPT,  NREVRS,  NOISE.  NCARD,  NPLOT.  NREAD  variable  names 
(  713  )  format 

* 

MMAXF  is  a  control  variable  which  is  used  to  read 
in  the  correct  number  of  points  in  the  sound  velocity 
profile.  (  see  Entry  2.  -  next  Entry  ) 

NOUTPT  is  the  number  of  output  signal  plots. 

NREVRS  is  both  a  flag  and  a  pointer.  If  non-zero,  it 
signals  that  there  will  be  pulse  reversal,  and  it  points 
to  the  output  signal  which  is  reversed  and  used  as  a 
source  function. 

• 

NOISE  is  a  control  flag  that,  when  set  to  zero. 

(  be  careful!  )  zeroes  out  the  output  signal  which  is 
outside  the  UMIN  to  UMAX  interval.  This  eliminates  the  '' 

noise  outside  the  specified  time  interval,  which  can  £ 

be  severe  at  long  ranges.  * 


NCARD  is  a  control  variable  that,  when  set  non-zero, 
allows  the  program  to  put  out  signal  data  on  unit  7  in  a  form 
is  readable  as  input  by  another  program.  A  zero  value  will  p 
any  output  on  unit  7.  This  is  normally  used  to  make  graphic 
(  Originally,  on  the  IBM  360/91.  unit  7  was  the  punched  card 
hence  the  name  NCARD.  ) 


wh  ict 
reveni 
plots, 
unit. 


a 


j 


NPLOT  is  a  control  variable  that,  when  set  non-zero, 
allows  the  program  to  put  out  printer  plots  of  the 
output  signals. 


NREAD  is  a  control  variable  that,  when  set  non-zero, 
directs  the  program  to  interpret  Entry  2  data  (  see  below  ) 
as  thickness  input  instead  of  the  usual  depth  input. 


2.  Entry  2.  MMAXF  lines  (  max.  of  49  lines  )  •1 

((  V(M.  N) .  N-1.7  ),  M»l.  MMAXF  >  variable  name 

(  7D10.  0  )  format 


V  is  an  array  that  contains  the  sound  velocity  profile. 

For  each  point  in  the  profile,  and  on  each  line.  ? 

(  one  point  per  line  ).  the  program  expects:  * 

a.  a  depth  value,  in  meters,  except  when  NREAD  is  set  non-zero. 

in  which  case  the  program  expects  a  thickness  value,  in  meters  » 

b.  a  sound  (  pressure  )  velocity  value,  in  meters  per  second  » 

c.  a  sound  <  shear  )  velocity  value,  in  meters  per  second, 

if  applicable 

d.  a  corresponding  attenuation  value  for  pressure 

e.  a  corresponding  attenuation  value  for  shear.  n 

if  applicable 

f.  a  density  value,  in  grams  per  c.  c.  .  or  specific  density 

g.  a  roughness  value,  ir  applicable 

Note:  No  depth  value  can  be  less  than  a  previous 

depth  value. 


3.  Entry  3.  one  line  -4 

KSF,  LRF.  IRES.  LS.  KKK.  IBEAM.  NW,  RSPACE.  PO.  C.  CO 

variable  namesJ 

(  715,  4D10.  0  )  format 


KSF  is  a  pointer  to  the  source  depth  value. 

It  points  to  one  of  the  points  in  the  sound  velocity 
profile;  the  one  which  is  the  source  depth  value. 

LRF  is  a  pointer  to  the  detector  (  or  receiver  ) 
depth  value.  it's  function  is  similar  to  KSF. 

IRES  is  a  control  variable  that,  when  set  non-zero, 
requires  the  program  to  taper  the  response  function 
at  each  end. 


LS  is  the  number  of  points  used  in  the  wave  integration. 

KKK  is  a  control  variable  that  indicates: 

a.  If  0  or  2*  this  is  the  last  case. 

b.  If  -1  or  -2.  the  next  case  starts  with  Entry  1. 

(  new  profile  >. 

c.  If  +1  or  +3 <  the  next  case  starts  with  Entry  3. 

The  second  value  in  each  case  (  2 *  -2*  +3  >•  will  direct 
the  program  to  take  each  signal  and  compress  it  back  to 
its  source.  A  plot  of  this  compressed  pulse  follows 
each  output  signal. 

IBEAM  is  a  control  variable  that.  when  set  non¬ 
zero*  requires  the  program  to  taper  the  wave  integrand 
at  each  end. 

NW  is  a  flag  that*  when  set  non-zero*  directs  the 
program  to  print  out  selected  integrand  plots. 

RSPACE  is  the  range  spacing  between  output  signals* 
in  kilometers.  However*  the  program  can  only  approximate 
this  value. 

PO  is  the  source  level.  For  this  program  we  normally 
use  a  value  of  lOOOO. 

C  is  the  lower  bound  sound  velocity  in  the  wave 
integration*  in  maters  per  second. 

CO  is  the  upper  bound  sound  velocity  in  the  wave 
integration*  in  meters  per  second. 

Entry  4*  one  line 

FREQO.  NX*  INTF*  FREQM.  DT*  RO*  FREQS.  UMAX*  UMIN 

variable  names 

(  F10.  0.  219*  6F10.  0  )  format 

FREQO  is  the  lower  bound  of  the  frequency  integration, 
in  Hertz. 

NX  is  the  number  of  FFT  points  and  the  number  of 
points  in  each  output  signal. 

INTF  is  a  variable  that  specifies  how  often*  within  the 
frequency  integration,  the  wave  integrand  will  be  plotted  on  the 
printer.  (  every  INTF'th  frequency  ).  Enter  NX+1  for  no  wave 
integrand  plots. 

FREQM  is  the  upper  bound  of  the  frequency  integration, 
in  Hertz. 

DT  is  the  time  difference  between  the  successive  points  in 
the  output  signals*  in  secs. 

RO  is  the  initial  range  or  the  range  of  the  first 
output  signal*  in  kilometers. 

FREQS  is  the  sound  source  frequency*  in  Hertz. 

UMAX  is  a  sound  velocity  (  meters/sec.  )  that  is  used  to 
specify  the  lower  bound  of  the  output  signal  if  some  of  the  signal 
is  considered  noise.  <  See  NOISE*  Entry  I.  also  see  CRED*  Entry  9  ) 

UMIN  is  similar  in  function  to  UMAX*  but  used  to 
specify  the  upper  bound. 


Entry  5.  on*  1  in? 

REFLTO,  ALPHAO,  ATTENO,  RHOO,  H»  CRED.  BOUNCE 
(  6010  0.  15  > 


variable  nan* 
format 


REFLTO  is  th*  RMS  surfac*  roughness  value/  in  meters. 
If  this  value  is  non-zero#  the  program  explicitly  sets 
the  next  three  values  to  zero. 

ALPHAO  is#  if  present,  th*  sound  velocity  in  the 
fluid  upper  half  space#  in  meters  per  second.  This  is 
th*  variable  that  determines  whether  th*  present  case  has 
either  a  fluid  or  vacuum  upper  half  space. 

ATTENO  is.  if  present,  th*  attenuation  in  th* 
fluid  upper  half  space. 

RHOO  is.  if  present,  the  density  of  the  fluid 
upper  half  space,  in  grams  par  c. c. 

H  is.  if  present,  th*  distance  between  the  source 
and  detector  in  th*  fluid  upper  half  space.  To  us*  this 
option: 

a.  The  first  two  sound  velocity  depths  must  be  zero. 

b.  There  must  be  some  variation  in  the  pressure 
sound  velocity  profile. 

c.  There  must  not  be  any  shear  sound  velocity. 

d.  ALPHAO  and  H  must  be  present. 

CRED  is  th*  reducing  velocity.  For  each  output  signal, 
the  starting  time  is  offset  from  TO  at  rang*  O  (  point 
of  origin  )  bu  a  time  equal  to  that  rang*  divided  by  CRED. 
This  value  defaults  to  CO.  for  a  blank  o 


to  that  rang*  divided  by  CRED. 
for  a  blank  or  zero  entry. 


BOUNCE  is  a  control  variable  that,  when  set  non-zero, 
directs  the  program  to  calculate  only  th*  effects  of  a  certain 
number  of  bounces.  BOUNCE  being  that  number.  This  is  somewhat 
similar  to  th*  'H  variable'  output  and  should  not  be  used  in 
conjunction  with  that  feature. 
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C.  Operating  Procedure* 

The  primary  computer  at  Lamont-Doherty  Geological  Observatory 
is  a  VAX  11/780.  which  runs  under  the  VMS  operating  system. 

Listed  below  are  four  distinctly  different  operations  that 
are  used  to  process  IMPULS22B. 


Comp i 1 ing-Link ing 
a.  Compiling 


Compiling  (  FORT  command  ) 

Compiling  the  IMPULS22B  program  on  the  VAX  11/780 
has  always  been  done  with  the  following  switches 
on: 

1.  /CHECK  Forbids  array  pointer  overflow. 

2.  /NO I 4  Integer  variables  are  set  to  2  bytes. 

3.  /NQF77  Where  the  compiler  needs  to  choose 
between  Fortran  77  and  Fortran  4+»  the  latter 
is  chosen. 

4.  /LIST  A  compiler  listing  is  mad*  and 
stored  on  the  disk  as  IMPULS22B.  LIS. 

Linking  (  LINK  command  > 


b.  Linking  (  LIh 

Linking  is  straight-forward.  No  extra  routines 
are  linked  with  IMPULS22B. 


2.  Running 

IMPULS22B  is  normally  run  through  a  batch  stream, 
in  particular  the  SLOW  batch  stream,  since  this  is  the 
cheapest  rout*,  especially  at  night.  However,  there  is 
another  batch  stream  with  higher  priority  (  and  higher 
cost  >  that  can  be  used  for  quicker  results,  the 
SYS4BATCH  batch  stream.  In  order  to  use  either  of  these 
two  batch  streams: 

a.  The  executable  file  (  IMPULS22B.EXE  >  must  be  ready 
to  us*.  The  LINK  operation  satisfies  this  requirement. 

b.  The  data  file  must  also  be  ready  to  use.  This  is 
normally  prepared  by  the  user  using  the  editor. 

(  and  the  IMPULS22B  USER'S  MANUAL,  if  necessary  > 
at  a  terminal. 

c.  Submit  the  command  file.  IMPULS22B.COM.  to  the 
batch  stream  (  or  batch  queue  >  using  a  command 
such  as  one  of  these: 

SUBMIT/NAME-ITB17MAR/QUEUE-SL0W/AFTER*20: 00  IMPULS22B 
SLOW  queue.  Job  starts  at  8:00  PM. 

SUBM IT/NAME- I TB 1 7MAR / QUEUE aSLOW  I MPULS22B 
SLOW  queue.  Job  starts  now 

SUBM I T /NAME* I TB 1 7MAR  IMPULS22B 

SYS4BATCH  queue  (  the  'regular'  queue  >. 

Job  starts  now 

d.  The  following  command  file,  also  prepared  using 
the  editor,  is  a  typical  file  submitted,  using  one 
of  the  previous  submit  commands. 

♦ASSIGN  ITB17MAR  SYS*PRINT  (1) 

♦ASSIGN  ITB17MAR  SYSMUTPUT  <2> 

♦ASSIGN  DRO:  ITB17MAR.  PLT  F0R007  (3) 

♦ASSIGN  DRO:  ITB17MAR.  OUT  F0R006  (4) 

♦ASSIGN  DRO: STRATE. DAZ  F0R009  (3) 

♦TYPE  IMPULS22B  <6> 

♦COPY  DRO: STRATE.  DAZ  DRO:  ITB17MAR.  DAT  (7) 

♦RUN  IMPULS22B  (8) 

This  run  is  called  the  ITB17MAR  run.  'ITB'  stands  for 
IMPULSE  TWENTY  TWO  B.  ' 1 7MAR '  indicates  that  this  is  the 

17th  run  in  March. 


n.T 


POT 


The  first  two  lines  in  this  commend  file  put  the 
log  file  on  the  disk  instead  of  the  printer. 

The  third  line  puts  unit  7  output  into  file  ITB17MAR.  PLT. 
The  fourth  line  puts  unit  6  output  into  file  ITB17MAR.  OUT. 
The  fifth  line  draws  unit  5  input  from  STRATE.  DAZ. 

The  sixth  line  copies  IMPULS22B.  LIS  to  ITB17MAR.  LIS. 

The  seventh  line  copies  STRATE.  DAZ  to  ITB17MAR.DAT. 

The  eighth  line  runs  program  IMPULS22B. 

A  command  file  such  as  this  yields  a  file  group; 
each  file  having  the  same  name  but  a  different 
extension.  These  files  are  then  treated 
collectively.  They  are  printed  together  and  stored 
on  tape  together.  Such  a  group  is*  for  example: 

ITB17MAR.DAT  a  copy  of  the  data  file 

ITB17MAR.  LIS  the  program  listing  file 

ITB17MAR.LOG  the  log  file 

IT817MAR.  OUT  the  output  file 

ITB17MAR.  PLT  the  numerical  output  file 

The  user  normally  gets  a  complete  printout  of  the 
run  by  typing  the  command: 

PRINT  ITB17MAR.  *» * 

on  the  terminal. 

Copying  Runs  to  Tape  <  Archiving  ) 

Per iod ical ly.  it  is  necessary  to  clear  the  old  runs 
off  the  disk  and  put  them  on  tape.  This  is  a  very  simple 
operation  and  proceeds  as  follows: 

a.  Initialize  a  tape.  Give  it  a  name  that  is  somehow 
indicative  of  the  present  month.  For  example: 

INITIALIZE/DENSITY-1600  MT1 :  FEBMAR 

This  tape  would  have  runs  made  in  February  and  March. 

This  step  is  not  necessary  if  the  user  is  adding  runs 
to  a  tape  that  is  already  partly  filled.  In  such  a 
case,  be  careful  not  to  initialize! 

b.  Mount  the  tape  with  the  MOUNT  command. 

MOUNT  MT1 :  FEBMAR 

This  is  not  the  same  as  physically  mounting  the  tape 
on  the  tape  drive,  which  is  assumed  to  have  taken 
place  previously. 

c.  Copy  all  the  run  files  from  disk  to  tape  with  the 
COPY  command.  For  example  in  this  case: 

COPY  *FEB.  *i  *  MT1: 

COPY  *MAR.  *;  *  MT1: 

would  copy  all  the  February  and  all  the  March  runs 
onto  the  tape  mounted  on  tape  drive  1. 

d.  For  future  reference,  it  is  good  to  get  a  list  of 
names  of  the  files  on  the  tape.  This  is  accomplished 
as  follows: 

DIRECTORY/PRINT  MT1 : 


The  printout  resulting  from  this  command  does  not 
have  the  tape  name  on  it.  so  the  user  should  write 
this  on  immediately. 

If.  for  some  reason,  this  printout  is  lost,  the  tape 
can  be  remounted  and  the  directory  can  be  taken 
again.  But  this  is  no  reason  to  be  careless  with 
tne  original  directory  printout. 
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e.  Dismount  the  tape  with  the  dismount  command: 

DISMOUNT  MT1 : 

f.  Delete  the  files  on  the  disk  that  have  just  been 
copied  to  the  tape.  In  this  case: 

DELETE  *FEB.  *;  * 

DELETE  *MAR.  *;  * 

4.  Distributing  the  Program  by  Tape 

The  program  and  test  cases  are  distributed  using 
tapes  with  the  same  format  as  the  archival  tapes.  This 
is  the  Files-11  format*  peculiar  to  Digital  Equipment 
Corporation.  The  tapes  from  this  system  have  a  density 
of  1600  bpii  and  a  blocksize  of  2048  bytes. 


If  any  further  questions  arise*  feel  free  to  call: 

Dr.  H.  Kutschale  (914)  339-2900  x382 

Mr.  L.  Carroll  (914)  359-2900  *221 


ITB15MAR.  DATi  3 


6— APR- 1983  15:29 


TEST  CASE  DATA 


04  21 

0 

1  1  1 

0 

000.  0 

1490.  0 

.  000005 

1.  03 

050.  0 

1490.  0 

. 000005 

1.  03 

4200. 

1490.  0 

. 000005 

1.  03 

4200. 

1520.  0 

. 000010 

1.  03 

2 

02 

0  256 

0  00  001  05. 

10000.  0 

2.  441407 

512  40 

24.414063  0.0080 

50.  0 

8.  0 

0.  0 

0.  00 

00.  0  00.  00 

000.  0 

1519 

1488.  0 


1519.  999 


GE  K 


ITB18MAR.  DAT; 2 


6-APR-1983  15:29 


TEST  CASE  DATA 


04  21  0 

000.  0 
050.  0 
4200. 

4200. 

2  02 
2.  441407 
0.  0 


1  1  1  < 
1490.  0 
1490.  0 
1490.  0 
1520.  0 
0  256 


512 
0.  00 


600.  0 
0  00 


.  000005 
.  000005 
. 000005 
.  000010 
001  05. 


1.  03 
1.  03 
1.  03 
1.  03 
10000.  0 


1488.0 


40  24.  414063  0.  0080  50.  0  8.  0 

00.0  00.00  000.0  1519.999 


1519.  999 


-*.4 


ITB19MAR.  DAT;  2 


6-APR-1983  15:30 


TEST  CASE  DATA 
04 21  01110 

000. 0  1490.  0  000005 

000.  0  1490.  0  000005 

4200.  1490.  0  ‘  000005 

4200.  1520.  0  „  !  000010 

2  02  0  128  0  00  001  05 

2.441407  512  40  24.  414063  0  0080  °5' 

0.  0  1490.  0  .  000005  01.  03 


1.  03 
1.  03 
1.  03 
1.  03 

10000.0  1488.0 

50.  0  8.  0 

25.00  1519.999 


1519.  999 
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Introduction 


IMPULS11G.  our  more  complete  impulse  program/  was  completed  in  October 
of  1982.  It  is  an  extension  of  IMPULS10G/  the  farmer  interim  standard 
impulse  program. 

IMPULS110  has  these  features  and  options  for  the  user: 


1. 

A 

2. 

Op 

hor 

pre 

3. 

A 

The 

a. 

b. 


The  FFP  method 

The  Normal  Mode  method 


usual 


a. 

b. 

c. 

d. 

e. 

f. 


l. 


A  detailed  sound  velocity  profile 
Choice  of  source  depth 
Choice  of  detector  depth 

Choice  of  frequency  band  and  source  frequency 
Choice  of  wave  number  integrand  limits 
Optional  beaming/  as  it  is  affected  by 
integrand  tapering  (  in  the  FFP  mode  ) 

Optional  surface  roughness 

Choice  of  either  fluid  or  vacuum  upper  half 
space 

Possibility  of  the  source  and  detector 
the  fluid  upper  half  space 
Optional  internal  roughness 
Choice  of  initial  range/  range  spacing/ 
range  count 


in 


and 


Output. 

a. 

b. 

c. 


This  includes: 

A  printout  of  all  input 
The  pulse  signal  plots  (  usually  ) 

A  transmission  loss  plot  for  the  range 
ulse  signals  (  Normal  Mode  method  only 


pulse  signa 
Additional  Optional  Tables  and  Plots 


of 

> 


the 


These 


includi 


a. 


b. 


Dispersion  curves/  which  includes: 

Phase  velocities  and  group  velocities 
The  wave  integrand  plots  of  the  FFP  for  certain 
frequencies 


Note:  All  short  plots  except  the  wave 
are  repeated/  the  repeated  plot  being 
copy  of  the  first  and  having  the  same 
the  first. 


integrands 
a  smoothed 
bounds  as 


6.  Numerical  output  of  the  program  for  further  processing  or  plottini 
Formerly/  this  was  called  punched  output/  since  it  went  out 
to  the  card  punch. 
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B. 


Input  Variables  and  Formats 

Program  input  for  IMPULS11G  includes  both  the  physical  model 
and  program  control.  This  section  of  the  manual  mill  list 
and  explain  each  line  of  the  input. 

1.  Entry  1.  one  line 

MMAXFT  NOUTPT,  NREVRS,  NOISE.  NCAROi  NPLOT  variable  names 

(  613  )  format 

MMAXF  is  a  control  variable  which  is  used  to  read 
in  the  correct  number  of  points  in  the  sound  velocity 
profile.  (  see  Entry  2.  -  next  Entry  > 

NOUTPT  is  the  number  of  output  signal  plots 

NREVRS  is  both  a  flag  and  a  pointer.  If  non-zero,  it 
signals  that  there  will  be  pulse  reversal,  and  it  points 
to  the  output  signal  which  is  reversed  and  used  as  a 
source  function. 

NOISE  is  a  control  flag  that,  when  set  to  zero. 

(  be  careful!  )  zeroes  out  the  output  signal  which  is 
outside  the  UMIN  to  UMAX  interval.  This  eliminates  the 
noise  outside  the  specified  time  interval,  which  can 
be  severe  at  long  ranges  for  the  FFP  mode. 

NCARD  is  a  control  variable  that,  whan  set  non-zero, 
allows  the  program  to  put  out  signal  data  on  unit  7  in 
a  form  which  is  readable  as  input  bu  another  program. 

A  zero  value  will  prevent  any  output  on  unit  7.  This 
is  normally  used  to  make  graphic  plots.  (  Originally, 
on  the  IBM  360/91,  unit  7  was  the  punched  card  unit, 
hence  the  name  NCARD.  ) 

NPLOT  is  a  control  variable  that,  when  set  non-zero, 
allows  the  program  to  put  out  printer  plots  of  the 
output  signals. 

2.  Entry  2,  MMAXF  lines  <  max.  of  49  lines  > 

((  V(M.  N),  N-1,7  >,  M-l.  MMAXF  >  variable  name 

(  7D10.  0  )  format 


V  is  an  array  that  contains  the  sound  velocity  profile. 
For  each  point  in  the  profile,  and  on  each  line, 

(  one  point  per  line  ),  the  program  expects: 

a.  a  depth  value,  in  meters 

b.  a  sound  (  pressure  )  velocity  value,  in  meters  per  second 

c.  a  sound  (  shear  )  velocity  value,  in  meters  per  second, 

if  applicable 

d.  a  corresponding  attenuation  value  for  pressure 

e.  a  corresponding  attenuation  value  for  shear, 

if  applicable 

f.  a  density  value,  in  grams  per  c. c. ,  or  specific  density 

g.  a  roughness  value,  if  applicable 

Note:  No  depth  value  can  be  less  than  a  previous 

depth  value.  The  program  will  stop  if  this  is  so. 


3.  Entry  3,  one  line 

KSF,  LRF,  IRES,  LS,  KKK, 

(  713,  4D10.  0  ) 


IBEAM,  NM,  RSPACE,  PO,  C,  CO 
variable  names 
f ormat 


KSF  is  a  pointer  to  the  source  depth  value. 

It  points  to  one  of  the  points  in  the  sound  velocity 
profile;  the  one  which  is  the  source  depth  value. 

LRF  is  a  pointer  to  the  detector  (  or  receiver  > 
depth  value.  It's  function  is  similar  to  KSF. 

IRES  is  a  control  variable  that,  when  set  non-zero, 
requires  the  program  to  taper  the  response  function 
at  each  end. 


I 
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LS  is  the  number  of  points  used  in  the  wave  integration. 

KKK  is  a  control  variable  that  indicates: 

a.  If  0  or  2*  this  is  the  last  case. 

b.  If  -1  or  -2*  the  next  case  starts  with  Entry  1 . 

(  new  profile  ). 

c.  If  +1  or  +3*  the  next  case  starts  with  Entry  3. 

The  second  value  in  each  case  (  2.  -2.  +3  )•  will  direct 
the  program  to  take  each  signal  and  compress  it  back  to 
its  source.  A  plot  of  this  compressed  pulse  follows 
each  output  signal. 

IBEAM  is  a  control  variable  that.  when  set  non¬ 
zero*  requires  the  program  to  taper  the  wave  integrand 
at  each  end. 

NM  is  the  maximum  number  of  modes  needed  for  the 
Normal  Mode  case. 

RSPACE  is  the  range  spacing  between  output  signals> 
in  kilometers.  For  the  Normal  Mode  methods  this  number 
is  exact.  For  the  FFP  method#  the  program  approximates 
this  value. 


T 

'J 


PO  is  the  source  level, 
use  a  value  of  10000.  0 


For  this  program  we  normally 


4. 


C  is  the  lower  bound  sound  velocity  in  the  wave 
integration*  in  meters  per  second. 

CO  is  the  upper  bound  sound  velocity  in  the  wave 
integration*  in  meters  per  second. 

Entry  4*  one  line 
FREQO.  NX*  INTF.  FREQM. 

(  F10.  0*  213.  6F10.  0  ) 


DT,  RO,  FREQS.  UMAX.  UMIN 

variable  names 
format 


FREQO  is  the  lower  bound  of  the  frequency  integration* 
in  Hertz. 

NX  is  the  number  of  FFT  points  and  the  number  of 
points  in  each  output  signal. 

INTF  is  a  variable  that  specifies  how  often*  within  the 
frequency  integration,  the  wave  integrand  will  be  plotted  on  the 
printer.  <  every  INTF'th  frequency  ).  Enter  NX+l  for  no  wave 
integrand  plots. 

FREQM  is  the  upper  bound  of  the  frequency  integration, 
in  Hertz. 

DT  is  the  time  difference  between  the  successive  points  in 
the  output  signals,  in  secs. 

RO  is  the  initial  range  or  the  range  of  the  first 
output  signal,  in  kilometers. 

FREQS  is  the  sound  source  frequency*  in  Hertz. 

UMAX  is  a  sound  velocity  (  meters/sec.  >  that  is  used  to 
specify  the  lower  bound  of  the  output  signal  if  some  of  the  signal 
is  considered  noise.  (  See  NOISE.  Entry  1.  also  sc*  CRED*  Entry  3  ■ 


UMIN  is  similar  in  function  to  UMAX* 
specify  the  upper  bound. 


but  used  to 
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5.  Entry  5/  one  line 

REFLTO.  ALPHAO.  ATTENO.  RHOO,  H*  CRED.  KHPV,  IDISP.  KURVE, 

IVEL  variable  names 

(  6010.  0.  419  )  format 

REFLTO  is  the  RMS  surface  roughness  value*  in  meters. 

If  this  value  is  non-nro*  the  program  explicitly  sets 
the  next  three  values  to  zero. 

ALPHAO  is*  if  present*  the  souno  velocity  in  the 
fluid  upper  half  space*  in  meters  per  second.  This  is 
the  variable  that  determines  whether  the  present  case  has 
either  a  fluid  or  vacuum  upper  half  space. 

ATTENO  is*  if  present*  the  attenuation  in  the 
fluid  upper  half  space. 

RHOO  is*  if  present*  the  density  of  the  fluid 
upper  half  space*  in  grams  per  c. c. 

H  is*  if  present*  the  distance  between  the  source 
and  detector  in  the  fluid  upper  half  space.  To  use  this 
option: 

a.  The  first  two  sound  velocity  depths  must  be  zero. 

b.  There  must  be  some  variation  in  the  pressure 
sound  velocity  profile. 

c.  There  must  not  be  any  shear  sound  velocity. 

d.  ALPHAO  and  H  must  be  present. 

CRED  is  the  reducing  velocity.  For  each  output  signal* 
the  starting  time  is  offset  from  TO  at  range  0  (  point 
of  origin  )  bu  a  time  equal  to  that  range  divided  by  CRED. 

This  value  defaults  to  CO*  for  a  blank  or  zero  entry. 

KHPV  is  the  control  variable  that  determines 
whether  output  will  be: 

a.  pressure*  if  zero. 

b.  vertical  displacement*  if  positive. 

c.  horizontal  displacement*  if  negative. 

IVEL  modifies  this  if  it  is  set.  (  see  below  > 

FFP  only. 

IDISP  is  a  control  variable  that  directs  the  program 
to  forego  output  signals  in  favor  of  dispersion  curves. 

In  this  case,  the  dispersion  curve  plot  runs  from 
FREQO  to  FREGM  (  min  to  max  freq.  >  (  entry  4  >  and  from 
C  to  CO  (  min  to  max  sound  velocities  >  (  entry  3  >.  This 
option  will  also  print  phase  velocity  and  group  velocity 
plots  for  up  to  9  modes*  and  punch  out  these  phase  and 
group  values  for  further  processing. 

KURVE  is  a  control  variable  that*  when  set  non-zero, 
signals  the  program  to  contour  the  wave  integrands.  This 
produces  a  plot  similar  to  the  dispersion  curve  plot 
above  while  still  allowing  signal  output.  FFP  only. 

IVEL  is  a  control  variable  that*  when  set  non-zero* 
alters  the  output  signal  definition  from  displacement 
to  velocity.  KHPV  must  be  set*  or  this  has  no  effect. 

FFP  only. 

Note: 

The  program  determines  the  computational  method* 

(  whether  FFr  or  Normal  Mode  >*  by  examining  all  of 
the  attenuation  values  and  the  surface  roughness  value. 

If  they  are  all  zero,  the  Normal  Mode  method  is  selected. 
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IMPULSl 1G  USER'S  MANUAL 


Operating  Procedures 

Tha  primary  computer  at  Lamont-Doherty  Geological  Observatory 
is  a  VAX  11/780.  which  runs  under  the  VMS  operating  system. 

Listed  below  are  four  distinctly  different  operations  that 
are  used  to  process  XMPULSllG. 


Comp ilina -Linking 
a.  Compiling 


Comp i 1 ing 
Compiling  the  I 
has  always  been 


<  FORT  c 
IMPULSl 1G  progr 
n  done  with  tne 


command  ) 

ram  on  the  VAX  11/780 
e  following  switches 


/CHECK 
/NO  1 4 
/NQF77 


Forbids  array  pointer 
Integer  variables  are 


3.  /NQF77  Where  the  compiler  needs  t< 
between  Fortran  77  and  Fortran  4+.  the 
is  chosen. 

4.  /LIST  A  compiler  listing  is  made 
stored  on  the  disk  as  IMPULSl IG.  LIS. 


er  overflow, 
re  set  to  2  bytes, 
needs  to  choose 
4+,  the  latter 


Linking 
Linking  is 
are  linked 


(  LINK 
straight-f orward. 
with  IMPULSl 1G. 


command  > 

No  extra  routines 


Runnin 


mg 

IMpULSHG  is  normally  run  through  a  batch  stream. 
. —  — .  .  a|ch  s- 


in  particular  the  SLOW  b 


tream.  since  this  is  the 


cheapest  route,  especially  at  night.  However,  there  i« 
another  batch  stream  with  higher  priority  (  and  higher 
cost  )  that  can  be  used  for  guicker  results,  the 


SYS4BATCH  batch  stream, 
two  batch  streams: 


In  order  to  use  either  of  these 


a.  The  executable  file  (  IMPULS11G.EXE  >  must  be  ready 
to  use.  The  LINK  operation  satisfies  this  requirement. 

b.  The  data  file  must  also  be  ready  to  use.  This  is 
normally  prepared  by  the  user  using  the  editor. 

(  and  the  IMPULSHG  USER'S  MANUAL,  if  necessary  ) 
at  a  terminal. 

c.  Submit  the  command  file.  IMPULS11G.COM.  to  the 
batch  stream  (  or  batch  queue  >  using  a  command 
such  as  one  of  these: 

SUBM I T  /NAME*  I  EG  1 90C  T / QUEUE»SLQW / AFTER  »20 :  00  IMPULSHG 
SLOW  queue.  Job  starts  at  8:00  P.  M. 

SUBM IT/NAME*! EG 150CT/GUEUE*SL0W  IMPULSl 1G 
SLOW  queue.  Job  starts  now 

SUBMIT/NAME-IEG19QCT  IMPULSHG 

SYSdBATCH  queue  (  the  'regular'  queue  ). 

Job  starts  now 

d.  The  following  command  file,  also  prepared  using 
the  editor,  is  a  typical  file  submitted,  using  one 
of  the  previous  submit  commands. 

•ASSIGN  IEG150CT  SYS*PRINT  (1) 

•ASSIGN  IEG15GCT  SYS*OUTPUT  (2) 

•ASSIGN  DRO: IEG150CT. PLT  F0R007  <3> 

•ASSIGN  DRO:  IEG190CT.  OUT  F0R006  (4) 

•ASSIGN  DRO: STANDARD.  DAZ  F0R009  (9) 

•TYPE  IMPULSHG  (6> 

•COPY  DRO:  STANDARD.  DAZ  DR0-.IEG190CT.DAT  (7) 

•RUN  IMPULSHG  <8> 

This  run  is  called  the  IEG190CT  run.  'IEG'  stands  for 
IMPULSE  ELEVEN  G.  '190CT'  indicates  that  this  is  the 
19th  run  in  October. 


The  first  two  lines  in  this  command  file  put  the 
log  file  on  the  disk  instead  of  the  printer. 

The  third  line  puts  unit  7  output  into  file  IEG1SOGT. PLT. 
The  fourth  line  puts  unit  6  output  into  file  IEG150CT.  OUT. 
The  fifth  line  draws  unit  5  input  from  STANDARD.  DAZ. 

The  sixth  line  copies  IMPULS1 1G. LIS  to  IEG150CT. LIS. 

The  seventh  line  copies  STANDARD.  DAZ  to  IEG13QCT.DAT. 

The  eighth  line  runs  program  IMPULS11G. 


A  command  file  such  as  this  yields  a  file  group; 
each  file  having  the  same  name  but  a  different 
extension.  These  files  are  then  treated 
collectively.  They  are  printed  together  and  stored 
on  tape  together.  Such  a  group  is,  for  example: 


IEG150CT.  DAT 

a 

IEG150CT.  LIS 

th 

IEG150CT.  LOG 

th 

IEG150CT.  OUT 

th 

IEG1SOCT.  PLT 

th 

copy  of  the  data  file 
e  program  listing  file 
e  log  file 
e  output  file 
e  numerical  output  file 


The  user  normally  gets  a  complete  printout  of  the 
run  by  typing  the  command: 

PRINT  IEG150CT.  *;* 


on  the  terminal. 

Copying  Runs  to  Tape  (  Archiving  ) 

Periodically!  it  is  necessary  to  clear  the  old  runs 
off  the  disk  and  put  them  on  tape.  This  is  a  very  simple 
operation  and  proceeds  as  follows: 

a.  Initialize  a  tape.  Give  it  a  name  that  is  somehow 
indicative  of  the  present  month.  For  example: 

INITI ALIZE/DENSITY=1600  MTl :  OCTNOV 

This  tape  would  have  runs  made  in  October  and  November. 
This  step  is  not  necessary  if  the  user  is  adding  runs 
to  a  tape  that  is  already  partly  filled.  In  such  a 
case,  be  careful  not  to  initialize! 

b.  Mount  the  tape  with  the  MOUNT  command. 

MOUNT  MTl:  OCTNOV 

This  is  not  the  same  as  physically  mounting  the  tape 
on  the  tape  drive,  which  is  assumed  to  have  taken 
place  previously. 

c.  Copy  all  the  run  files  from  disk  to  tape  with  the 
COPY  command.  For  example  in  this  case: 

COPY  *OCT.  *;  *  MTl: 

COPY  *NOV.  *;  *  MTl: 


would  copy  all  the  October  and  all  the  November  runs 
onto  the  tape  mounted  on  tape  drive  1. 

d.  For  future  reference,  it  is  good  to  get  a  list  of 
names  of  the  files  on  the  tape.  This  is  accomplished 
as  follows: 


DIRECTORY/PRINT  MTl: 

The  printout  resulting  from  this  command  does  not 
have  the  tape  name  on  it,  so  the  user  should  write 
this  on  immediately. 

If,  for  some  reason,  this  printout  is  lost,  the  tape 
can  be  remounted  and  the  directory  can  be  taken 
again.  But  this  is  no  reason  to  be  careless  with 
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NL 
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e.  Dismount  the  tap*  with  the  dismount  command: 

DISMOUNT  MT1  : 

f.  Delete  the  files  on  ths  disk  that  have  just  bean 
copied  to  tha  tape.  In  this  case: 

DELETE  *OCT.  *i  * 

DELETE  *NOV.  *i  * 

4.  Distributing  tha  Program  by  Tape 

Tha  program  and  test  cases  are  distributed  using 
tapes  with  the  same  format  as  the  archival  tapes.  This 
is  the  Files-11  format.  peculiar  to  Digital  Equipment 
Corporation.  The  tapes  from  this  sustem  have  a  density 
of  1600  bpi.  and  a  blocksize  of  204b  bytes. 


If  any  further  questions  arise*  feel  free  to  call: 

Dr.  H.  Kutschale  (914)  359-2900  x382 

Mr.  L.  Carroll  (914)  359-2900  *221 
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PELOSSX  USER'S  MANUAL 
Introduction 

PELOSSX/  our  PE  transmission  loss  program#  was  completed  in 
November  of  1982. 

PELOSSX  has  these  features  and  options  for  the  user: 

1.  The  PE  method  of  computation 


2. 


3. 


A  comprehensive  selection  of  modelling  parameters. 

These  include: 

a.  A  detailed  sound  velocity  profile 

b.  Choice  of  source  depth 

c.  Choice  of  detector  depth 

d.  Choice  of  frequency 

e.  Optional  ability  to  vary  the  velocity  profile  with  range 

Output.  This  includes: 

a.  A  printout  of  all  input 

b.  Optional  contour  plot  of  the  transmission  loss  in  a  field 
of  depth  vs.  range 

c.  Optional  page  plot  of  transmission  loss  vs.  range  for  a 
given  detector  depth 

Numerical  output  of  the  program  for  further  processing  or  plotting. 
Formerly/  this  was  called  punched  output#  since  it  went  out 
to  the  card  punch. 
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PELOSSX  USER'S  MANUAL 


B.  Input  Variables  and  Format*  j 

Program  input  for  PELOSSX  includes  both  the  physical  modal 
and  program  control.  This  section  of  the  manual  will  list 
and  explain  each  line  of  the  input.  J 

1.  Entry  1.  one  line  ! 

MMAX#  KKK#  MRMAX,  NFL AO »  IDUMM,  LVX#  DEPLR  variable  names 

(  613#  2X#  FIO.  0  )  format 

MMAX  is  a  control  variable  which  is  used  to  read  in  the 
correct  number  of  points  in  the  sound  velocity  profile#  using 
input  type  A.  (  see  Entry  4.  >  If  depths  and  velocities  are  input# 

(  input  type  B  )»  a  non-zero  value  is  needed  here  also  because 
the  depths  and  velocities  arc  converted  to  type  A  input.  MMAX 
must  be  a  <  power  of  2  )  plus  1  because  of  its  close  association 
with  the  FFT.  Maximum  allowable  input  value  is  129.  (  See  also  LVX 

below#  the  next  to  last  item  in  Entry  1.  > 

KKK  is  a  flag  that  determines  whether  or  not  this  is  the  last 
case.  Enter  zero  for  the  last  case#  or  non-zero  to  go  back  to  the 
top  and  go  through  the  program  again. 

MRMAX  is  the  number  of  sound  velocity  profiles.  This  is  used 
for  sound  velocity  profile  range  dependency.  Maximum  value  is  16. 

NFLAG  is  a  flag  that  directs  the  program  to  plot  a  single  pag< 
plot  of  transmission  loss  and/or  a  contour  plot  of  transmission 
loss.  A  negative  value  will  yield  only  a  contour  plot,  a  zero  valui 
will  yield  both  types  of  plot#  and  a  positive  value  will  yield  a 
single  page  (  one  detector  depth  )  plot  only. 

IOUMM  is  a  dummy  variable  at  present. 

LVX  is  both  a  flag  and  a  counter  value.  If  zero#  the  sound 
velocity  profile.  (  Entry  4.  below  >  consists  of  eight  sound 
velocity  values  per  line;  MMAX  values  per  profile.  This  is  input 
type  A.  If  LVX  is  greater  than  zero,  the  sound  velocity  profile 
consists  of  'LVX'  lines  of  depth  and  velocity  values.  This  is 
input  type  B. 

DEPLR  is  the  receiver  or  detector  depth,  in  meters. 

2.  Entry  2.  one  line 

N.  NSTEP.  FREQ.  ZMAX.  SNORM.  FGW,  DRKM,  GD 

variable  name' 

(  0F1O.  0  )  format 

N  is  the  number  of  points  in  the  refined  velocity  profile. 
Velocity  values  in  Entry  5  are  interpolated  to  yield  N  velocity 
values.  Also#  N  must  be  a  power  of  2#  since  it  is  the  number  of 
FFT  points.  Maximum  value  is  912. 

NSTEP  is  the  number  of  computational  steps  taken  in  range. 
NSTEP  times  DRKM  <  see  below  )  gives  the  total  range  of  the  plot. 

FREQ  is  th»  input  frequency#  in  Hertz. 

ZMAX  is  the  maximum  depth#  in  meters. 

SNORM  is  the  reducing  velocity#  in  meters  per  sec. 

FGU  is  an  ini t ia 1 i za t ion  factor  which  is  normally  set  to  1.0 

DRKM  is  the  length  of  the  spacing  between  computational 
points,  in  kilometers. 

GD  is  the  source  depth#  in  meters. 


& 
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variable  name 
format 


Entry  3<  from  on*  to  thraa  line* 

(  RPTS(I),  I  »  1,  MRMAX  ) 

(  8F10.  2  ) 

RPTS  is  an  array  that  contains  rang*  points.  Th*s*  range 
points  specify  the  ranges  at  which  the  program  switches  from 
on*  velocity  profile  to  the  next.  The  last  rang*  point  of  every 
run  should  be  out  further  than  the  maximum  rang*  (i.e.  NSTEP*DRKM>. 
A  run  with  only  on*  velocity  profile  still  needs  a  range  point 
entered. 

Entry  4.  Input  tup*  A>  UMAX  entries  (  max.  17  lines*  129  values  > 

(  C(l!>  I  *  l.MMAX  )  variable  name 

(  SFtO.  0  >  format 

•  C  is  an  arrau  that  has  velocity  values  in  it.  These  velocity 
values  are  assumed  to  be  at  equally  spaced  intervals  along  the 
velocity  profile  starting  at  the  top  (  the  surface  of  the  water  ) 
and  ending  at  ZMAX. 

Input  type  B>  LVX  lines  (  max.  87  lines*  87  points  > 

((  V(M*  N)  *  N-1.2  ),  M-l.  LVX  )  variable  name 

(  2F10.  4  >  format 

V  is  an  array  that  contains  the  sound  velocity  profile. 

For  each  point  in  the  profile*  and  on  each  line* 

(  one  point  per  line  )>  the  program  expects: 

a.  a  depth  value*  in  meters 

b.  a  sound  (  pressure  )  velocity  value*  in  meters  per  second 

Input  type  B  is  flagged  by  a  positive  value  of  LVX. 

Entry  4  is  repeated  rlRMAX  number  of  times. 


•-  ^  V'  »  ^  v  -/  s"  v  *. 
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Operating  Procedure* 


The  primary  computer  at  Lamont-Doherty  Geological  Observatory 
is  a  VAX  11/780*  which  runs  under  the  VMS  operating  system. 

Listed  below  are  four  distinctly  different  operations  that  . 
are  used  to  process  PELOSSX. 


Comp i ling -Linking 
a.  Compiling 


Compiling  (  PORT  command  > 

Compiling  the  PELOSSX  program  on  the  VAX  11/780 
has  always  been  done  with  the  following  switches 


1.  /CHECK  Forbids  array  pointer  overflow. 

2.  /N0I4  Integer  variables  are  set  to  2  bytes. 

3.  /N0F77  Where  the  compiler  needs  to  choose 
between  Fortran  77  and  Fortran  4+»  the  latter 
is  chosen. 

4.  /LIST  A  compiler  listing  is  made  and 
stored  on  the  disk  a*  PELOSSX.  LIS. 

b.  Linking  (  LINK  command  ) 

Linking  is  straight-forward.  No  extra  routines 
are  linked  with  PELOSSX. 


Forbids  array  pointer  overflow. 


Running 

PELOSSX  is 


normally  run  through  a  batch  stream* 
e  SLOW  batch  stream*  since  this  is  1 


in  particular  the  SLOW  batch  stream*  since  this  is  the 
cheapest  route*  especially  at  night.  However*  there  is 
another  batch  stream  with  higher  priority  (  and  higher 
cost  )  that  can  be  used  for  quicker  results*  the 
SYSPBATCH  batch  stream.  In  order  to  use  either  of  these 
two  batch  streams: 

a.  The  executable  file  (  PELOSSX.EXE  )  must  be  ready 

to  use.  The  LINK  operation  satisfies  this  requirement. 

b.  The  data  file  must  also  be  ready  to  use.  This  is 
normally  prepared  bu  the  user  using  the  editor* 

(  and  the  PELOSSX  USER'S  MANUAL*  if  necessary  ) 
at  a  terminal. 

c.  Submit  the  command  file*  PELOSSX.COM*  to  the 
batch  stream  (  or  batch  queue  >  using  a  command 
such  as  one  of  these: 

SUBMIT /NAME -P EL07NOV / QUEUE "SLOW / AFTER -20 :  00  PELOSSX 
SLOW  queue*  Job  starts  at  8:00  P.  M. 

SUBM I T /N AME-PEL07NOV / QUEUE-SLOW  PELOSS X 
SLOW  queue*  Job  starts  now 

SUBMIT /NAME-PEL07N0V  PELOSSX 

SYSPBATCH  queue  (  the  'regular'  queue  )* 

Job  starts  now 

d.  The  following  command  file*  also  prepared  using 
the  editor*  is  a  typical  file  submitted*  using  one 
of  the  previous  submit  commands. 

PASSION  PEL07N0V  SYSPPRINT  (1) 

PASSION  PEL07NOV  SYSPOUTPUT  (2) 

PASSION  DRO: PEL07NOV.  PLT  F0R007  <3> 

PASSION  DRO: PEL07N0V.  OUT  FOROO&  (4) 

PASSION  DRO: STANDARD.  DPE  F0R005  (9) 

PTYPE  PELOSSX  (6) 

PCOPY  DRO: STANDARD. DPE  DRO: PEL07N0V.  DAT  (7) 

PRUN  PELOSSX  (8) 

This  run  is  called  the  PEL07N0V  run.  'PEL'  stands  for 
PARABOLIC  EQUATION  LOSS.  '07N0V'  indicates  that  this  is  the 
7th  run  in  November. 
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The  first  two  lints  in  this  command  filt  put  tht 
log  filt  on  tht  disk  insttad  of  tht  printtr. 

The  third  lint  puts  unit  7  output  into  filt  PEL07N0V.  PLT. 
Tht  fourth  lint  puts  unit  6  output  into  filt  PEL07N0V.  OUT. 
Tht  fifth  lint  draws  unit  9  input  from  STANDARD.  DPE. 

Tht  sixth  lint  copits  PELOSSX.  C.IS  to  PEL07N0V.  LIS. 

Tht  seventh  lint  copits  STANDARD.  DPE  to  PEL07N0V.DAT. 

Tht  tighth  line  runs  program  PELOSSX. 

A  command  filt  such  as  this  yields  a  filt  group; 
each  filt  having  the  same  name  but  a  different 
extension.  These  files  art  then  treated 
collectively.  They  art  printed  together  and  stored 
on  tape  together.  Such  a  group  is<  for  example: 


PEL07N0V.  DAT 
PEL07N0V.  LIS 
PEL07N0V.  LOG 
PEL07NOV.  OUT 
PEL07N0V.  PLT 


a  copy  of  the  data  filt 
the  program  listing  filt 
the  log  filt 
the  output  filt 
the  numerical  output  fill 


Tht  user  normally  gats  a  complete  printout  of  the 
run  by  typing  the  command: 

PRINT  PEL07NQV.  *»  * 

on  the  terminal. 

3.  Copying  Runs  to  Tape  (  Archiving  > 

Periodically*  it  is  necessary  to  clear  the  old  runs 
off  the  disk  and  put  them  on  tape.  This  is  a  very  simple 
operation  and  proceeds  as  follows: 

a.  Initialize  a  tape.  Give  it  a  name  that  is  somehow 
indicative  of  the  present  month.  For  example: 

INITIALIZE/DENSITY-1600  MTl:  OCTNOV 

This  tape  would  have  runs  made  in  October  and  November. 
This  step  is  not  necessary  if  the  user  is  adding  runs 
to  a  tape  that  is  already  partly  filled.  In  such  a 
case,  be  careful  not  to  initialize! 

b.  Mount  the  tape  with  the  MOUNT  command. 

MOUNT  MT1 :  OCTNOV 

This  is  not  the  same  as  physically  mounting  the  tape 
on  the  tape  drive<  which  is  assumed  to  have  taken 
place  previously. 

c.  Copy  all  the  run  files  from  disk  to  tape  with  the 
COPY  command.  For  example  in  this  case: 

COPY  *OCT.  *j  *  MT1 : 

COPY  *NOV.  *;*  MTl: 

would  copy  all  the  October  and  all  the  November  runs 
onto  the  tape  mounted  on  tape  drive  1. 


mounted  on  tape  drive  1. 


d.  For  future  reference*  it  is  good  to  get  a  list  of 
names  of  the  files  on  the  tape.  This  is  accomplished 
as  follows: 

DIRECTORY/PRINT  MTl: 

The  printout  resulting  from  this  command  does  not 
have  the  tape  name  on  it*  so  the  user  should  write 
this  on  immediately. 

If.  for  some  reason*  this  printout  is  lost*  the  tape 
can  be  remounted  and  the  directory  can  be  taken 
again.  But  this  is  no  reason  to  be  careless  with 
the  original  directory  printout. 


-\  *V'. 'WW 
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Dismount  tho  top*  with  ths  dismount  command: 
DISMOUNT  MT1: 
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f.  Doisto  tho  files  on  tho  disk  that  havo  just  boon 
copied  to  tho  tape.  In  this  case: 

DELETE  *OCT.  *;  * 

DELETE  *NOV.  *i  * 

4.  Distributing  tho  Program  by  Tape 

The  program  and  tost  cases  are  distributed  using 
tapes  with  tne  same  format  as  the  archival  tapes.  This 
is  the  Files-11  format#  peculiar  to  Digital  Equipment 
Corporation.  The  tapes  from  this  system  have  a  density 
of  1600  bpi.  and  a  blocksize  of  2043  bytes. 


If  any  further  questions  arise#  feel  free  to  call: 

Dr.  H.  Kutschale  (914)  399-2900  x382 

Mr.  L.  Carroll  (914)  359-2900  x221 
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Introduction 

PEPULSEX.  our  PE  impulse  program,  was  completed  in  October  of  1982 
PEPULSEX  has  these  features  and  options  for  the  user: 

1.  The  PE  method  of  computation 

2.  A  comprehensive  selection  of  modelling  parameters. 

These  include: 

a.  A  detailed  sound  velocity  profile 

b.  Choice  of  source  depth 

c.  Choice  of  detector  depth 

d.  Choice  of  frequency  band  and  source  frequency 

e.  Choice  of  initial  range,  range  spacing,  and 

range  count 

f.  Optional  ability  to  vary  the  velocity  profile  with  range 

3.  Output.  This  includes: 

a.  A  printout  of  all  input 

This  includes  a  printer  plot  of  all  the  velocity  profiles. 

b.  The  pulse  signal  plots 

4.  Numerical  output  of  the  program  for  further  processing  or  plotting 
Formerly,  this  was  called  punched  output,  since  it  went  out 

to  the  card  punch. 


s’ 
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Input  Variables  and  Formats 


Program  input  for  PEPULSEX  includes  both  the  physical  model 
and  program  control.  This  section  of  the  manual  will  list 
and  explain  each  line  of  the  input. 


1.  Entry  1> 
MMAX,  KKK* 


<  1113  ) 


one  line 

MR MAX i  N,  LDUMM*  NTIME,  NOUTPT,  NREVRS,  NCARD.  NPLOT .  LV] 

variable  namei 


format 


MMAX  is  a  control  variable  which  is  used  to  read  in  the 
correct  number  of  points  in  the  sound  velocity  profile*  using 
input  type  A.  (  see  Entry  9.  )  If  depths  and  velocities  are  input* 

(  input  type  B  )*  a  non-zero  value  is  needed  here  also  because 
the  depths  and  velocities  are  converted  to  type  A  input.  MMAX 
must  be  a  (  power  of  2  >  plus  1  because  of  its  close  association 
with  the  FFT.  Maximum  allowable  input  value  is  129.  (  See  also  LVX 

below*  the  last  item  in  Entry  1.  ) 


KKK  is  a  flag  that  determines  whether  or  not  this  is  the  last 
case.  Enter  zero  for  the  last  case*  or  non-zero  to  go  back  to  the 
top  and  go  through  the  program  again. 


MRMAX  is  the  number  of  sound  velocity  profiles.  This  is  used 
far  sound  velocity  profile  range  dependency.  Maximum  value  is  16. 


N  is  the  number  of  points  in  the  refined  velocity  profile. 
Velocity  values  in  Entry  9  are  interpolated  to  yield  N  velocity 
values.  Also*  N  must  be  a  power  of  2*  since  it  is  the  number  of 
FFT  points  at  one  location  in  the  program.  Maximum  value  is  912. 


LDUMM  is  a  dummy  variable  at  present. 


NTIME  is  a  flag  that,  when  set  non-zero*  moves  the  signal  tim 
window  up  by  the  length  of  the  window.  In  essence*  the  reducing 
velocity  marks  the  and  of  the  pulse  signal  plot  instead  of  the 
beginning,  which  is  the  normal  case. 


NOUTPT  is  the  number  of  output  signal  plots. 


NREVRS  is  both  a  flag  and  a  pointer.  If  non-zero*  it 
signals  that  there  will  be  pulse  reversal*  and  it  points 
to  the  output  signal  which  is  reversed  and  used  as  a 
source  function. 


NCARD  is  a  control  variable  that*  when  set  non-zero* 
allows  the  program  to  put  out  signal  data  on  unit  7  in 
a  form  which  is  readable  as  input  by  another  program. 


a  form  which  is  readable  as  input  by  another  program. 
A  zero  value  will  prevent  any  output  on  unit  7.  This 


is  normally  used  to  make  graphic  plots.  (  Originally, 
on  the  IBM  360/91.  unit  7  was  the  punched  card  unit* 
hence  the  name  NCARD.  > 


NPLOT  is  a  control  variable  that*  when  set  non-zero* 
allows  the  program  to  put  out  printer  plots  of  the 
output  signals. 


LVX  is  both  a  flag  and  a  counter  value.  If  zero*  the  sound 


velocity  profile*  (  Entry  9*  below  >  consists  of  eight  sound 
velocity  values  per  line*  MMAX  values  per  profile.  This  is  input 
type  A.  If  LVX  is  greater  than  zero*  the  sound  velocity  profile 
consists  of  'LVX'  lines  of  depth  and  velocity  values.  This  is 
input  type  B. 


Entry  z*  one  line 
RO,  RsPACE*  DEPLR.  ZMAX,  SNORM,  FGW*  DRKM.  GD 


one  line 


(  8F10.  0  ) 


variable  nami 
format 


RO  is  the  initial  range  or  the  range  of  the  first 
output  signal*  in  kilometers. 


RSPACE  is  the  range  spacing  between  output  signals* 
in  kilometers.  Unless  this  is  a  multiple  of  DRKM*  T  see  below  > 


the  program  can  only  approximate  this  value. 

DEPLR  is  the  receiver  or  detector  depth*  in  meters. 


218 


ZMAX  is  the  maximum  depth*  in  meters. 

SNORM  is  the  reducing  velocity,  in  meters  per  sec.  For  each 
output  signal,  the  starting  time  is  offset  from  TO  at  range  0 
(  point  of  origin  )  by  a  time  equal  to  that,  range  divided  by  SNORM. 

FGW  is  an  ini t ial i zat ion  factor  which  is  normally  set  to  1.0 

DRKM  is  the  length  of  the  spacing  between  computational 
points,  in  kilometers. 

GD  is  the  source  depth,  in  meters. 

Entry  3.  from  one  to  three  lines 

UMIN.  UMAX.  (  RPTS(I).  I  ■  1.  MR MAX  )  variable  names 

(  8F10.  0  )  format 


to  specify  the  lower  bound  or  the  output  signal. 

UMIN  is  similar  in  function  to  UMAX,  but  used  to 
specify  the  upper  bound. 

RPTS  is  an  array  that  contains  range  points.  These  range 
points  specify  the  ranges  at  which  the  program  switches  from 
one  velocity  profile  to  the  next.  The  last  range  point  of  every 
run  should  be  out  further  than  the  last  output  pulse  signal. 

A  run  with  only  one  velocity  profile  still  needs  a  range  point 
entered. 


Entry  4.  one  line 
FREQO.  NX.  DT.  FREQS.  PO.  FREQM 

(  F10.  0.  15.  5X,  4F10.  2  ) 


variable  names 
format 


FREQO  is  the  lower  bound  of  the  frequency  integration, 
in  Hertz. 

NX  is  the  number  of  FFT  points  and  the  number  of 
points  in  each  output  signal. 

DT  is  the  time  difference  between  the  successive  points  in 
the  output  signals,  in  sec. 

FREQS  is  the  sound  source  frequency,  in  Hertz. 

PO  is  the  source  level.  For  this  program  we  normally  use 
a  value  of  10000. 

FREQM  is  the  upper  bound  of  the  frequency  integration, 
in  Hertz. 

Entry  5.  Input  type  A.  MMAX  entries  (  max.  17  lines.  129  values  ) 
(  C<  15.  I  ■  1. MMAX  )  variable  name 

(  8F10.  0  ) 


variable  name 
format 


C  is  an  array  that  has  velocity  values  in  it.  These  velocity 
values  are  assumed  to  be  at  equally  spaced  intervals  along  the 
velocity  profile  starting  at  the  top  (  the  surface  of  the  wmter  ) 
and  ending  at  ZMAX. 

Input  type  B.  LVX  lines  (  max.  37  lines.  87  points  > 

(C  V<M,  N).  N*l.  2  7.  M«l.  LVX  >  variable  nam< 


(  2F10.  4  ) 


variable  name 
format 


V  is  an  array  that  contains  the  sound  velocity  profile. 
For  each  point  in  the  profile,  and  on  each  line. 

(  one  point  per  line  ).  the  program  expects: 

a.  a  depth  value,  in  meters 

b.  a  sound  (  pressure  )  velocity  value,  in  meters  per  second 

Input  type  B  is  flagged  by  a  positive  value  of  LVX. 

Entry  5  is  repeated  MRMAX  number  of  times. 


PEPULSEX  USER'S  MANUAL 


219 


Operating  Procedures 

The  primary  computer  at  Lamont-Doherty  Geological  Observatory 
is  a  VAX  11/780*  which  runs  under  the  VMS  operating  system. 

Listed  below  are  four  distinctly  different  operations  that 
are  used  to  process  PEPULSEX. 


1. 


Comp iling -Linking 

a.  Compiling  (  FORT  command  ) 

Compiling  the  PEPULSEX  program  on  the  VAX  11/780 
has  always  been  done  with  the  following  switches 
on: 


1.  /CHECK  Forbids  array  pointer  overflow. 

2.  /NO I 4  Integer  variables  are  set  to  2  bytes. 

3.  /NQF77  Where  the  compiler  needs  to  choose 
between  Fortran  77  and  Fortran  4+.  the  latter 
is  chosen. 

4.  /LIST  A  compiler  listing  is  made  and 
stored  on  the  disx  as  PEPULSEX.  LIS. 


b.  Linking  (  LINK  command  ) 

Linking  is  straight-forward.  No  extra  routines 
are  linked  with  PEPULSEX. 


2.  Running 

PEPULSEX  is  normally  run  through  a  batch  stream, 
in  particular  the  SLOW  batch  stream*  since  this  is  the 
cheapest  route*  especially  at  night.  However*  there  is 
another  batch  stream  with  higher  priority  (  and  higher 
cost  >  that  can  be  used  for  quicker  results,  the 
SYSPBATCH  batch  stream.  In  order  to  use  either  of  these 
two  batch  streams: 


a.  The  executable  file  (  PEPULSEX.EXE  >  must  be  ready 

to  use.  The  LINK  operation  satisfies  this  requirement. 

b.  The  data  file  must  also  be  ready  to  use.  This  is 
normally  prepared  by  the  user  using  the  editor* 

(  and  the  PEPULSEX  USER'S  MANUAL,  if  necessary  > 
at  a  terminal. 

c.  Submit  the  command  file.  PEPULSEX.COM*  to  the 
batch  stream  (  or  batch  queue  >  using  a  command 
such  as  one  of  these: 

SUBMIT/NAME-PEP09N0V/QUEUE-SL0W/AFTER-20: 00  PEPULSEX 
SLOW  queue*  Job  starts  at  8:00  P.  M. 

SUBM I T/NAME-PEP05NOV/GUEUE-SLOW  PEPULSEX 
SLOW  queue*  Job  starts  now 

SUBMIT /NAME*PEP05N0V  PEPULSEX 

SYSPBATCH  queue  (  the  'regular'  queue  )> 

Job  starts  now 


d.  The  following  command  file*  also  prepared  using 
the  editor,  is  a  typical  file  submitted*  using  one 
of  the  previous  submit  commands. 


PASSION  PEP09N0V  SYSPPRINT  (1) 
PASSION  PEP09N0V  SYSPOUTPUT  (2) 
PASSION  DRO: PEP09N0V.  PLT  F0R007  (3) 
PASSION  DRO: PEP09NOV.  OUT  F0R006  (4) 
PASSION  DRO: STANDARD.  DP2  F0R009  (9) 
PTYPE  PEPULSEX  <6> 
PCOPY  DRO: STANDARD.  DP2  DRO: PEP09N0V.  DAT  (7) 
PRUN  PEPULSEX  (8) 


This  run  is  called  the  PEP09N0V  run.  'PEP'  stands  for 
PARABOLIC  EQUATION  PULSE.  '09N0V '  indicates  that  this  is  the 
9th  run  in  November. 
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The  first  two  lints  in  this  command  filt  put  tht 
log  filt  on  tht  disk  insttad  of  tht  printer. 

Tht  third  lint  puts  unit  7  output  into  filt  PEPOSNOV.  PLT. 
Tht  fourth  lint  puts  unit  6  output  into  filt  PEPOSNOV.  OUT. 
Tht  fifth  lint  draws  unit  9  input  from  STANDARD.  DP2. 

Tht  sixth  lint  copits  PEPULSEX.  LIS  to  PEPOSNOV.  LIS. 

Tht  stvtnth  lint  copits  STANDARD.  DP2  to  PEPOSNOV.  DAT. 

Tht  tighth  lint  runs  program  PEPULSEX. 

A  command  file  such  as  this  yields  a  filt  group* 
each  filt  having  tht  samt  name  but  a  difftrtnt 
extension.  These  files  art  then  treated 
collectively.  They  art  printed  together  and  stored 
on  tape  together.  Such  a  group  is*  for  example: 

PEPOSNOV.  DAT  a  copy  of  the  data  filt 

PEPOSNOV.  LIS  the  program  listing  filt 

PEPOSNOV.  LOO  the  log  file 

PEPOSNOV.  OUT  the  output  filt 

PEPOSNOV.  PLT  the  numerical  output  filt 

Tht  user  normally  gets  a  complete  printout  of  the 
run  by  typing  the  command: 

PRINT  PEPOSNOV.  *j# 

on  the  terminal. 

3.  Copying  Runs  to  Tape  (  Archiving  ) 

Periodically*  it  is  necessary  to  clear  the  old  runs 
off  the  disk  and  put  them  on  tape.  This  is  a  very  simple 
operation  and  proceeds  as  follows: 

a.  Initialize  a  tape.  Give  it  a  name  that  is  somehow 
indicative  of  the  present  month.  For  example: 

INITIALIZE/DENSITY»1600  MTl :  OCTNOV 

This  tape  would  have  runs  made  in  October  and  November. 
This  step  is  not  necessary  if  the  user  is  adding  runs 
to  a  tape  that  is  already  partly  filled.  In  such  a 
case*  be  careful  not  to  initialize! 

b.  Mount  the  tape  with  the  MOUNT  command. 

MOUNT  MT1 :  OCTNOV 

This  is  not  the  same  as  physically  mounting  the  tape 
on  the  tape  drive*  which  is  assumed  to  have  taken 
place  previously. 

c.  Copy  all  the  run  files  from  disk  to  tape  with  the 
COPY  command.  For  example  in  this  case: 

COPY  #OCT.  *i*  MTl: 

COPY  #NQV.  *»  *  MTl: 

would  copy  all  the  October  and  all  the  November  runs 
onto  the  tape  mounted  on  tape  drive  1. 

d.  For  future  reference*  it  is  good  to  get  a  list  of 
names  of  the  files  on  the  tape.  This  is  accomplished 
as  follows: 

DIRECTORY/PRINT  MTl: 

The  printout  resulting  from  this  command  does  not 
have  the  tape  name  on  it*  so  the  user  should  write 
this  on  immediately. 

If*  for  some  reason*  this  printout  is  lost,  the  tape 
can  be  remounted  and  the  directory  can  be  taken 
again.  But  this  is  no  reason  to  be  careless  with 
the  original  directory  printout. 
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e.  Dismount  ths  taps  with  ths  dismount  command: 

DISMOUNT  MT1 : 

f.  Dslsts  ths  filss  on  ths  disk  that  havs  just  bssn 
copisd  to  ths  taps.  In  this  cass: 

DELETE  *OCT.  *;  * 

DELETE  *NOV.  *j  * 

4.  Distributing  ths  Program  by  Taps 

Ths  program  and  tost  casss  ars  distributsd  using 
tapss  with  ths  sams  format  as  ths  archival  tapss.  This 
is  ths  Filss-li  format*  psculiar  to  Digital  Equipmsnt 
Corporation.  Ths  tapss  from  this  sustsm  havs  a  dsnsity 
of  1-.600  bpi*  and  a  blocksizs  of  2048  bytss. 


If  any  furthsr  qusstions  ariss*  fssl  frss  to  call: 

Dr.  H.  Kutschala  (914)  339-2900  *382 

Mr.  L.  Carroll  (914)  339-2900  *221 
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129  0  1128  44  0  21  0  1  1  12  DATA  FOR  PEPULSEX 

50.0  005.0  61.0  3500.00  1499.0  1.0  1.0 

1490.0  1505.  900.0 

12.  50000  2048  0.  002  50.  10000.  0  102.  0000 

00.0  1494.0  0.000  0.00001  0.00000  1.03 

61.0  1494.7  0.000  0.00001  0.00000  1.03 

91.4  1495.2  0.000  0.00001  0.00000  1.03 

330.0  1496.8  0.000  0.00001  0.00000  1.03 

800.0  1487.9  0.000  0.00001  0.00000  1.03 

1950.0  1495.6  0.000  0.00001  0.00000  1.03 

3400.0  1511.7  0.000  0.00001  0.00000  1.03 

3400.0  1511.7  0.000  0.00010  0.00000  1.03 

3450.0  1511.7  0.000  0.00010  0.00000  1.03 

3450.0  1511.7  0.000  0.00010  0.00000  1.03 

3500.0  1511.7  0.000  0.00100  0.00000  1.03 

3500.0  1511.7  0.000  0.00100  0.00000  1.03 
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50.  0 
1490.  0 
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00.  0 
27.  344 
61.  0 
32.  031 
91.  4 
109.  375 
136.  719 
330.  0 
800.  0 
1950.  0 
3400.  0 
3500.  O 
00.  0 
27.  344 
61.  0 
82.  031 
91.  4 
109.  375 
136.  719 
330.  0 
800.  0 
1950.  0 
3400.  0 
3500.  0 
00.  O 
27.  344 
61.  O 
82.  031 
91.  4 
109.  375 
136.  719 
330.  0 
800.  0 
1950.  O 
3400.  0 
3500.  0 
00.  0 
27.  344 
61.  0 
82.  031 
91.  4 
109.  375 
136.  719 
330.  0 
800.  0 
1950.  0 
3400.  0 
3500.  0 
00.  0 
27.  344 
61.  0 
82.  031 
91.  4 
109.  375 
136.  719 
330.  O 
800.  0 
1950.  0 
3400.  0 
3500.  0 


5128  44 
005.  0 
1505. 
0  2048 

1494. 
1494. 

1494. 

1495. 

1495. 

>  1495. 

>  1495. 

1496. 
1487. 

1495. 
1511. 
1511. 
1494. 
1494. 
1494. 
1494. 

1494. 

\  1495. 

>  1495. 

1496. 
1487. 

1495. 
1511. 
1511. 
1494. 
1494. 
1494. 
1494. 

1494. 

>  1494. 

}  1495. 

1496. 
1487. 

1495. 
1511. 
1511. 
1494. 
1494. 

1494. 
1493. 

1493. 
5  1493. 

?  1494 

1496. 
1487. 

1495. 
1511. 
1511. 

1494. 
1494. 
1493. 
1493. 
1492. 

5  1492 

9  1494. 

1496 
1487 
1495 
1511 
1511 
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2 

32 
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8 
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6 
7 

7 
0 

314 

6 
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9 

02 
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8 
9 
6 
7 
7 

.  O 
314 
.  4 
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.  3 
.  52 
.  004 
.  8 
.  9 
.  6 
.  7 
.  7 
.  0 

.  314 
.  1 

.  846 
.  4 
.  82 
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.  8 
.  9 
.  6 
.  7 
.  7 
.  0 

.  314 
.  7 

i.  046 

:.  2 

!.  92 
.  104 

>.  8 
'.  9 

•.  6 

.  7 
.  7 


1  0  1 
61.  0 
050.  0 
0.  002 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
0.  000 
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3500.  00 
065.  0 
50. 

0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00100 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00100 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00100 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00100 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00001 
0.  00100 


1499.  0 
080.  0 
10000.  0 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
0.  00000 
o.  00000 
0. 00000 
0.  00000 
0.  00000 
o.  ooooo 
0.  ooooo 
o.  ooooo 
o.  ooooo 
o.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 
o.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 
0.  ooooo 


l.  0  l.o 

095.  0  900.  0 

102.  0000 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
1.  03 
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L29  0 

1128  12  0 

3  0  1 

1  25 

50.  0 

050.  0 

100. 

4200.  00 

1483.  0 

1.  0 

1434.  0 

1456. 

900.  0 

2.  441407  512 

0.  008 

8.  0 

10000.  0 

24.  41407 

000.  0 

1434.  0 

. 000005 

1.  0 

050.  0 

1434.  0 

. 000005 

1.  0 

100.  0 

1437.  0 

. 000005 

1.  0 

150.  0 

1440.  0 

. 000005 

1.  0 

200.  0 

1443.  0 

. 000005 

1.  0 

250.  0 

1446.  4 

. 000005 
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300.  0 
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. 000005 

1.  0 

350.  0 

1453.  2 
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1.  0 
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1456.  0 
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1.  0 

800.  0 

1460.  0 

. 000005 

1.  0 

1000. 

1463.  0 

. 000005 

1.  0 

1280. 

1467.  0 

. 000005 

1.  0 

1500. 

1472.  0 

.  000100 

1.  0 

1760. 

1476.  0 
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1.  0 

2000. 
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.  000100 
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2280. 

1484.  0 
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1.  0 
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1489.  0 
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1.  0 
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1494.  0 
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1.  0 
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.  001000 

1.  0 
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. 010000 

1.  0 

3520. 
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.  010000 

1.  0 

3760. 

1511.  0 

.  010000 

1.  0 
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1515.  0 

.  010000 

1.  0 

4000. 

1600.  0 

.  010000 

1.  0 

4200. 

1600.  0 

. 010000 

1.  0 

1.  0 
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trated  from  field  data  and  the  PFPP.  Hydroacoustic  signals  from  underwater 
explosions  that  have  propagated  over  the  Arctic  abyssal  plains  coononly 
display  marked  frequency  dispersion  in  pulses  that  are  bottom-interacting  and 
that  arrive  after  the  SOFAR  signal.  In  the  infrasonic  band  of  2  to  20  Hz,  the 
temporal  dispersion  for  each  pulse  that  has  interacted  with  the  flat  bottom  of 
the  plain  can  be  nearly  as  strong  as  that  observed  in  the  SOFAR  signal  for  the 
first  mode.  However,  the  bottom-interacting  pulses  correspond  to  a  coherent 
summation  of  many  higher-order  normal  modes  in  a  channel  bounded  above  by  the 
ocean  surface  and  below  by  the  zone  of  increasing  velocity  in  the  upper  400  m 
of  the  bottom  sediment Using  normal-mode  theory  and  the  Multiple  Scattering 
Pulse  Fast  Field  Program  (MSPFFP) ,  we  have  analyzed  the  dispersion  and  pulse 
shapes  and  have  derived  the  acoustic  properties  of  the  bottom  in  the  Pole, 
Barents,  and  Mendeleyev  Abyssal  Plains.  The  principal  properties  of  the 
bottom  controlling  the  propagation  are  congressional  velocity,  density,  and 
attentuation.  In  contrast,  the  ice  layer  has  a  negligible  effect  on  the 
dispersion  of  the  observed  waves.  The  effect  on  pulse  compression  of  this 
frequency  dispersion  of  the  bottom-interacting  signals  was  simulated 
numerically,  using  predistorted  waveforms  matched  to  the  dispersion  of  the 
SOFAR  channel  at  specified  ranges.  We  include  a  users  manual  (with  test 
computations)  for  our  computer  codes.  We  hope  that  this  manual  will  be  useful 
to  the  many  users  of  our  codes. 
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